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The quantum excitations in glasses have long presented a set of puzzles for condensed matter 
physicists. A common view is that they are largely disordered analogs of elementary excitations 
in crystals, supplemented by two level systems which are chemically local entities coming from 
disorder. A radical revision of this picture argues that the excitations in low temperature glasses 
are deeply connected to the energy landscape of the glass when it vitrifies: the excitations are 
not low excited states built on a single ground state but locally defined resonances, high in the 
energy spectrum of a solid. According to a semiclassical analysis, the two level systems involve 
resonant collective tunneling motions of around two hundred molecular units which are relics 
of the mosaic of cooperative motions at the glass transition temperature T g . The density of 
states of the TLS is determined by T g and the mosaic's length scale, which is a weak function 
of the cooling rate. The universality of phonon scattering in insulating glasses is explained. The 
Boson Peak and the plateau in thermal conductivity, observed at higher temperatures, are also 
quantitatively understood within the picture as arising from the same cooperative motions, but 
now accompanied by thermal activation of the mosaic's vibrational modes. The dynamics of 
some of the local structural transitions have significant quantum corrections to the semiclassical 
picture. These corrections lead to a deviation of the heat capacity and conductivity from the 
standard tunneling model results and explain the anomalous time dependence of the heat capacity. 
Interaction between tunneling centers contributes to the large and negative value of the Griineisen 
parameter often observed in glasses. 
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I. INTRODUCTION 

During the past several decades, it has been gradually 
recognized in the condensed matter and materials sci- 
ence community that amorphous materials, while sharing 
many characteristics with the more common crystalline 
solids, represent a distinct solid state of matter. On the 
one hand, glasses exhibit rigidity and elastic response on 
humanly relevant time scales, thus qualifying them as 
solids for many practical purposes. In fact, until the rel- 
atively recent advent of systematic studies of the materi- 
als' response to mechanical and electromagnetic pertur- 
bation, as well as of their detailed microscopic structure, 
the only commonly known distinct attributes of amor- 
phous substances had been their optical properties and 
the low magnitude and isotropic character of their ther- 
mal expansion. Those properties still undcrgird the main 
technological importance of amorphous materials. On 
the other hand, there are many ways in which glasses 
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are fundamentally different from crystals. This is most 
noticeable in their properties at cryogenic temperatures. 

We presently know very well that an amorphous solid 
is in reality a liquid caught locally in a small set o f 
metastable free energy minima (|Xia and Wolvnesl 12000). 
each of which are separated from the much lower free 
energy crystalline arrangement by high barriers. There- 
fore, the glass transition, as manifest in the laboratory, 
is not strictly speaking a phase transition in the regu- 
lar thermodynamic sense and is not accompanied by a 
symmetry change or appearance of a free energy singu- 
larity. In contrast, a liquid that was cooled below its 
melting point fast enough so as to avoid crystallization 
- i.e. has become supercooled - experiences a crossover 
to (highly viscous) activated transport. As the temper- 
ature is lowered further, the relaxation barriers grow in 
a very dramatic fashion thus confining the molecules in 
their metastable arrangements long enough to give the 
appearance of shear elasticity in the sample on the tech- 
nologically relevant frequency scales. A quantitative un- 
derstanding of the physics behind the glass transition has 
recently been achieved with the random first order tran- 
sition (RFOT) theory o f glasses (jKirkpatrick et aZJ . Il989t 
IXia and Wolvnesl 12000). This theory has provided a mi- 
croscopic picture of molecular motions in supercooled liq- 
uids, such as first principle predictions of the length scales 
of these motions and the cooperativity lengths and the 
barrier heights of the activated transport. At any given 
time, a supercooled liquid is a mosaic of cooperatively re- 
arranging regions, whose size becomes larger as the tem- 
perature is lowered. This article describes how the RFOT 
theory also provides the necessary microscopic input to 
understand the cryogenic anomalies observed in glasses. 

In spite of the absence of periodicity, glasses exhibit, 
among other things, a specific volume, interatomic dis- 
tances, coordination number and local elastic modulus 
comparable to those of crystals. Therefore it has been 
considered natural to consider amorphous lattices as 
nearly periodic with the disorder treated as a pertur- 
bation, often-times in form of defects, so such a study is 
not futile. This is indeed a sensible approach, as even 
the crystals themselves are rarely perfect, and many of 
their useful mechanical and other properties are deter- 
mined by the existence and mobility of some sort of de- 
fects as well by interaction between those defects. Nev- 
ertheless, a number of low temperature phenomena in 
glasses have persistently evaded a microscopic model- 
free description along those lines. A more radical re- 
vision of the concept of an elementary excitation on top 
of a unique ground state is necessary ifLubchenkd 2002: 
iLubchenko and Wolvnesl l200lLl2003a|) . 

Let us give a brief historical overview of some of the 
most outstanding issues in low temperature amorphous 
state physics. It was already noted in the 1960s that the 
thermal conductivity of amorphous solids is significantly 
lower than that of crystals. A low-temperature exper- 
imentalist using epoxy in his apparatus knew that its 
thermal conductivity at liquid helium temperatures went 



roughly as constant xT 2 , where the constant was prac- 
tically the same for other amorphous substances as well 
l)Andersorill999|) . Surprisingly, this had not particularly 
alarmed anyone, even though one would not a priori ex- 
pect low temperature properties of disordered solids to 
be different from crystals, as the appropriate thermal 
phonon length is much larger than the molecular scale 
which was presumed to characterize the relevant hetero- 
geneity scale. It w as not until Zeller and Pohl published 
their classic paper llZeller and Pohl!ll97lD that it became 
generally known that both the heat capacity and thermal 
conductivity of glasses were significantly different from 
those of crystals, and that these anomalies were corre- 
lated. The heat capacity turned out to be approximately 
linear in temperature and larger than the T 3 phononic 
contribution up to temperatures ~ 10 K. The challenge 
to the theorists was soon met by the so called Standard 
Tunn eling Model (STM) {Anderson et a^J . ll972t|Phillipsl 
1972), in which one assumes that due to a disordered 
pattern of molecular bonds in glasses, there are a num- 
ber of defects in the lattice (something like "loose" atoms 
or "dangling bonds" ) , which have two alternative posi- 
tions in space separated by a sufficiently low tunneling 
barrier. At low temperatures, the dynamics of such a 
system is described well by a two-level system (TLS) 
hamiltonian. If one assumes that the spectral density 
of these TLS's is flat, one recovers the linear heat capac- 
ity. One also finds that the inverse mean free path of a 
thermal phonon due to resonant scattering off the TLS's 
is equal to Z~ fp oc T, which implies thermal conductiv- 
ity k ~ I J2u Cp h M ZmfpM c s oc T 2 . Here, C P h(w) is 
the heat capacity of a phonon mode of frequency lu and 
c s is the speed of sound (one assumes here that heat is 
carried primarily by phonons, which was experimentally 
demonstra ted explicitly four years late r by Zaitlin and 
Anderson (|Zaitlin and AndersonL IT"975l) ). Note that the 
resonant character of phonon scattering implies that the 
scattering cross-section of low-frequency phonons would 
be independent of the scatterer size, but would scale with 
the phonon wavelength (squared) itself. Therefore no 
knowledge of scatterer 's microscopic details are needed. 
Rather, only a single coupling parameter is needed to es- 
timate the magnitude of scattering at lo w temperatures . 
The STM did prove to be very successful [|PhillipsllT98l|) . 
as it predicted, among other things, nonlinear sound ab- 
sorption due to the saturation of the resonant absorp- 
tion an d the phonon echo, both of which were later ob- 
serve d {Golding and Graebneit Il976t iHunklinger et all 
Il976l) . In spite of these successes, the microscopic na- 
ture of these defects had remained unknown, although 
there later appeared several indications in the literature 
that the tunneling centers are not single atom entities 
but rather involve motions within larger groups of atom s 
ijGuttman and Rahmarlll986l:lMon and Ashcro"ftlll978tK 
On the experimental front, there had been a growing 
amount of evidence that the number of these additional 
excitations and their coupling to the phonons are corre- 
lated and also depend on T g ijRavchaudhuri and Pohl 
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Il98lt iRevnolds. Jrl ll97St ll98Gf) . which culminated 
in the observation made by Freeman and Anderson 
l|Freeman and Anderson! Il986lb that the heat conductiv- 
ities of all studied insulating glasses, if scaled by elastic 
constants, fall onto the same line in two regions, con- 
nected by a non-universal flat piece corresponding to the 
so call ed "plateau". In Fig U we show a facsimile of Fig. 2 
from IjFreeman and Anderso n. 1986) that demonstrates 
this heat conductivity universality. The lower tempera- 




FIG. 1 Scaled thermal conductivity («) data for several 
amorphous materials is shown. The horisontal axis is tem- 
perature in units the Debye temperature To. The vertical 

axis scale K = B ^ p . The value of To is somewhat uncer- 
tain, but its choice made in (iFreeman and Anderson! Il986t) 
is strongly supported by that it yields universality in the 
phonon localization region. The solid lines are calculated 
using k ~ | J2uj Cph(o' )£ m fp(a>')cs w jth £ m f p /A = 15 and 
^mfp/A = 1 respectively ijFreeman and Anderson! Il98fit) . 

ture straight line corresponds to the value ~ 150 of the ra- 
tio of the thermal phonon mean free path l m f p to the ther- 
mal Debye wave-length A = hc s /ksT. This region spans 
roughly 1.5 decades in temperature between several mK 
(lowest T accessed so far for the heat conductivity mea- 
surements) and to 1-10 K, depending on the substance. 
The short linear region at higher temperatures (20-60 
K) corresponds to Z m f p /A ~ 1, which actuall y implies 
complete phonon localization ijGraebner et all Il986ji ac- 
cording to the heuristic loffe-Riegel criterion. This im- 
plies, among other things, that one can no longer use 
kinetic theory expressions for heat transfer at these tem- 
peratures, as a diffusive mechanism must prevail . The 



intermediate region ("plateau") is usually observed be- 
tween 1 and 30 K, and does not scale with the Debye 
temperature and speed of sound. The standard tun- 
neling model of non-interacting two-level systems men- 
tioned above is normally applied to the region where 
^mfp/A ~ 150, that is generically bel ow 1 K. The u ni- 
versality of Z m f p / A can be boiled down <|Phillipsl.lT98H) to 
the universality of the following combination of parame- 

— 2 — 

tcrs: P-^j, where P is the spectral and spatial density 

of the TLS's (empirically ~ 10 45±1 J _1 to -3 ), g is cou- 
pling to the elastic strain on the order of eV, p is the 

2 

mass density (for reference, ~ p( ?2 r 3 would be the inter- 
action strength between such TLS's at distance r from 
each other). Now, if the defects involved the motion of 
only a single atom, one would reasonably assume that the 
value of their spectral density and coupling to the lattice 
or their combination would be very strongly material de- 
pendent. Even though P and g 2 vary within almost two 
orders of magnitude (still surprisingly little), the combi- 

_ 2 

nation P-^-r is constant within 50% for different materials 

pc- 

(p and c 2 s vary considerably as well). It certainly takes 
a stretch of imagination to thi nk that this i s merely a 
coincidence, as pointed out in llLegg ett|. |199 1|). In 1988, 
Yu and Leggett proposed l|Yu and Leggettl I1988T1 that 
the density of states of the TLS might itself be a result 
of dipole-dipole interactions between some original non- 
renormalized excitations. In short, this idea is motivated 
by the observation that for TLS coupled to the phonons 
with strength g, the coefficient at the dipole-dipole inter- 
action term g 2 / pc 2 has dimensions energy times volume. 
Therefore the interaction induced renormalized density 
of states P has to be the inverse of g 2 / pc 2 with a coeffi- 
cient, hence the universality of Pg 2 / pc 2 for different ma- 
terials. However, it so far has not proved possible to use 
their approach to justify the value of that coefficient to 
yield the experimental Z m f p /A ~ 150. This is surprising, 
since one expects such a simple dimensional argument to 
be very robust. (Several other studies of the universal- 
ity (|Coppersmitbl Il99l! iMeissner and Spitzmanii Il98lj) 
were undertaken at the time, that followed the paper b y 
Freeman and Anderson ijFreeman and Andersoii H986Ll 
There has been subsequent work applying a renormal- 
izatio n group style calculation to a system of interacting 
TLS ijBurin and Kagarlll996jl . but it seems from the re- 
sults that renormalizations are rel evant only at ultr a-low 
temperatures (/iK's and below) ijNeu et all Il997(l . In 
spite of the di fficulties i n jus t ifying the strong inter ac- 
tion scenario aruzzol Il994t iLubchenko and Wolvnesl 
|200C|) . the works l|Leggettl Il99lt lYu and Leggettl Il988|) 
that first challenged the standard TLS paradigm remain 
the main conceptual motivation behind the present pa- 



1 This idea that the heat was transferee! b y a random walk was 
used early on by Einstein ( Eins teinl I1911D to calculate the ther- 
mal conductance of crystals but, of course, he obtained num- 
bers much lower than those measured in the experiment. As we 
now know, crystals at low enough T support well defined quasi- 
particles - the phonons - which happen to carry heat at these 



temperatures. Ironically, Einstein never tried his model on the 
amorphous solids, where it would be applicable in the £ m fp/A ~ 1 
regime. 
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per. One stresses however that the idea that the ob- 
served coupling constant, which is quite small, could 
be a result of some original "bare" strong interaction, 
is consistent with the microscopic theory if we argue 
it is the molecular interactions behind the glass tran- 
sition itself which become "renormalized" in a somewhat 
unexpected fashion. This mic roscopic theory suggests 
l)Lubchenko and Wolvne4l200l[) that the phenomenolog- 
ical two-level systems are discrete energy levels repre- 
senting resonantly accessible local degrees of freedom 
that exist in glasses due to the possibility of collective 
transitions between alternative structural configurations 
of compact regions encompassing roughly 200 molecular 
units. The theory of glassy ergodicity breaking shows 
the spectrum of these excitations is nearly flat and the 
density of states scales with the inverse glass transition 
temperature T s , echoing the excitation spectrum of a 
random energy model (REM) with that glass transition 
temperature. Furthermore, the transitions are an alter- 
native mode of motion that must be in equilibrium with 
phononic excitations at T g . This equilibrium requirement 
makes one realize that TLS-phonon coupling g, T g and 
the material's elastic constants are intrinsically related. 
The universality of the / m f p / A ratio is a consequence of 
this relationship reflecting the non-equilibrium charac- 
ter of the glassy state. The structural transitions, that 
become tunneling two-level systems at cryogenic temper- 
atures, exist because a glassy sample, when it falls out of 
equilibrium, resides in a metastable configuration chosen 
from a very high density of states. The sample is bro- 
ken up into a mosaic of dynamically cooperative regions. 
Alternatively speaking, the energy landscape is local in 
nature; that is rearrangements of compact regions will 
not change the structural state of the rest of the sam- 
ple, but only deform the surrounding regions weakly and 
purely elastically. A (small) fraction of these rearrange- 
ments requires overcoming only a very low barrier and 
can therefore occur even down to sub-Kelvin tempera- 
tures. The tunneling occurs by consecutive molecular 
displacement within the cooperativity length established 
at T g . The consecutive motion of atoms is conveniently 
visualized as a domain wall separating the two alternative 
local structural states, moving through the local region. 

The thermal conductivity plateau has traditionally 
been considered by most workers a separate issue from 
the TLS. In addition to the rapidly growing magnitude 
of phonon scattering at the plateau, an excess of den- 
sity of states is observed in the form of the so called 
"bump" in the heat capacity temperature dependence 
divided by T 3 . The plateau is interesting from several 
perspectives. For one thing, it is non-universal if scaled 
by the elastic constants (say u)t> and c s ). It is, how- 
ever, located between two universal regions and it is im- 
portant to understand which other scales in the problem 
determine its location and shape. The excitations that 
give rise to the dramatically increased phonon absorption 
at the corresponding frequencies have been circumstan- 
tially associated with the excitations observed as the so 



called Boson Peak (BP), directly seen in the inelastic X- 
ray and neutron scattering experiments, also observed 
in the optical Brillouin and Raman scattering measure- 
ments. These experimental developments date well into 
90-s and became possible, in the neutron spectroscopy 
case, due to the improved resolution in the neutrons' ve- 
locity detection, combined wit h the ability to g enerate 
higher energy incident beams ijForet et all ^96) . Sim- 
ilarly, meV resolution was needed to utilize the X-ray 
scattering technics to discern the small inelastic wings 
on th e sides of the strong elastic peak l)Benassi et all 
1996). The term "Boson Peak" comes from the fact 
that its intensity scales roughly according to the Bose- 
Einstein statistics. The extraction of the density of states 
from the spectra is unfortunately mod el dependent, and 
those models can be roug hly divided l|Pilla et all l2000|) 
into the ones where the Boson peak signifies the energy 
scale on the edge o f phonon localization, as promoted in 
ijForet et alJ . ll996|) . and those following the other school 
of thought which asserts that these modes are propagat- 
ing even well above the fre quency of the BP as supported 
by the interpretation in ijPilla et all l2000[) . As far as 
theoretical interpretation is concerned, it is our impres- 
sion that most of theories of the Boson Peak, existing 
until recently, have postulated a sort of spatial hetero- 
geneity in an otherwise per fectly elastic medium (see a 
partial list of references in (Gri eera et all l200l|) '). with 
the notable exception o f the soft-potential m odel (SPM) 
l)Buchenau et all Il992t iKarpov et aR Il983(l . It is, of 
course, always possible to recover the observed magni- 
tude of the heat capacity excess at the BP temperatures 
by a particular choice of parameters. While a contribu- 
tion of the lattice disorder to the density of states un- 
doubtedly exists and can be very significant (see, for ex- 
ample, sixrrrik^tiojisofsihca's heat capacity by Horbach 
at el. ijHorbach et all\l999h ). we must note that if amor- 
phous lattices were purely harmonic, the phonon absorp- 
tion at the BP frequencies would be of the Rayleigh type 
and should be significantly lower than ob served in the 
experiment llAndersonL Il98ll Uoshil fl979h . There must 
be internal resonances present in the bulk, that scatter 
phonons inelastically. Though phcnomcnologically intro- 
duced, this feature is present, for example, in the soft- 
potential model. An analysis of the higher temperature 
behavior of the tunneling transitions that give rise to the 
TLS at su bKelvin energies was provided in t he RFOT ap- 
proach in l)Lubchenko and Wolvneslf2003a(l . When these 
transitions occur at high enough temperature, the do- 
main wall separating the two alternative states can have 
its surface vibrations thermally excited. The large degen- 
eracy of these vibrational states, characteristic of a two 
dimensional membrane, that accompany the underlying 
structural transition, is sufficient to account for the en- 
hancement of phonon scattering at the plateau, as com- 
pared to the TLS regime. Finally, the superposition of 
the domain wall vibrations on the underlying tunneling 
transition leads to an excess of density of states that re- 
produces well the bump in the heat capacity (these com- 
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pound excitations we call "ripplons"). We therefore ar- 
rive at a unified physical picture that allows a unified 
quantitative explanation of previously seemingly unre- 
lated mysteries in the TLS regime and at the higher, 
plateau energies. 

The paper is organized as follows: the first section 
outlines the basics of the RFOT theory and then pro- 
ceeds in applying that theory to understanding the ori- 
gin of the tunneling centers in amorphous solids. The 
spectrum of the two-level systems, their coupling to the 
phonons and the origin of the universality of phonon scat- 
tering are then discussed. Additionally, we show how 
details of the derived TLS' tunneling amplitude distri- 
bution lead to a deviation of T dependence of the heat 
capacity from a strict linear form. The second section 
explains how the high energy vibrational excitations (rip- 
plons) of the tunneling interfaces gives rise to an excess 
of states which exhibits itself as the heat capacity bump 
and yields the rapidly rising phonon scattering at these 
higher energies. A short discussion of the relaxational 
absorption from these excitations is given and its fre- 
quency dependent part is derived. The contents of these 
first two sections are, for the most part, a detailed ac- 
coun t of the calculations underlying two ea rlier brief let- 
ters l|Lubchenko and Wolvnestl200lll2003a|) that have re- 
ported our explanation of the low temperature anomalies 
in glasses within a semiclassical approach. The third and 
fourth chapters are comprised of new results. There, we 
establish that, while not altering the main conclusions of 
the semiclassical picture, a purely quantum phenomenon 
of level mixing and repulsion has an observable effect on 
the density of states of the tunneling centers at low T. 
Finally, the interaction between tunneling centers, me- 
diated by phonons, is estimated and this is argued to 
make a significant contribution to the negative thermal 
expansivity (and thus a negative Griineisen parameter) 
observed in many amorphous materials. 



II. OVERVIEW OF THE CLASSICAL THEORY OF THE 
STRUCTURAL GLASS TRANSITION 

From a physicist's perspective, a theory of the glass 
transition describes what happens to a liquid when 
it is cooled down sufficiently but is not observed to 
crystallize. To a mathematician, this is a generalized 
problem of packing compact interacting objects of com- 
parable size given a specific constraint on the density 
distribution (it is not periodic) and total energy of the 
system. A nearly complete conceptual, microscopic 
picture of the amorphous state has emerged in the 
cours e of the two l a st dec ades (|Stoe^s^l_aiid_WolyjiesJ, 
19841) jSingh ef all Il985l) llKirkpatrick and Wolvnesi 
1987alblP llKirkpatrick and Thirumalai 1987alt : 
Kirknatrick et all Il989ft (IXia and Wolvnes I200C . 
2001alH) llLubchenko and Wolvnesi 120011 l2003alr . 
2004). This framework has lead to a unified, quan- 
titative understanding of many seemingly unrelated 



phenomena in supercooled liquids above and below 
the glass transition. The glasses we consider form 
at temperatures where quantum effects are small so 
classical statistical mechanics is used. We review such a 
classical glass transition in what follows. 

First, we make several comments on the phenomenol- 
ogy of supercooled liquids. Strictly speaking, these are 
nonequilibrium systems: When cooled sufficiently slowly, 
most simple liquids will crystallize at a temperature just 
below the melting temperature T m . Randomly atactic 
polymers become glassy but presumably never crystal- 
lize. The melting point is defined as the temperature at 
which the liquid and crystal free energies are equal. Cool- 
ing the liquid at least a bit below T m is necessary to create 
a free energy driving force so as to make the nucleation 
barrier finite and to allow the system to equilibrate. The 
crystal, once formed is different from the liquid in several 
ways, e.g. it scatters X-rays at precise angles and it is 
anisotropic. Crucially for us, a crystal supports trans- 
verse sound waves, at all frequencies (including u> = 0, 
thence the crystal retains its shape). In contrast, the 
supercooled liquid is a finite lifetime state since crystal- 
lization will eventually occur by nucleation. However, 
the growth of crystalline nuclei, inside the liquid, is sub- 
ject to the slowing of all motions in liquids. Owing to 
this dramatic slowing of liquid motions upon lowering 
the temperature, one can supercool the liquid substan- 
tially below its melting point, which is the key to forming 
glasses. The extra nucleation barrier ensures there is ad- 
equate time to study the properties of the supecooled 
noncrystalline state. Local structures in supecooled liq- 
uids persist for some time, call it l/w c . This time is 
longer than the time it takes to establish a Maxwell dis- 
tribution of velocity, which is at most a few vibrational 
periods. Such an amorphous system will support trans- 
verse waves at frequencies u > u> c , just as a crystal would, 
but will in contrast exhibit a liquid like, equilibrium re- 
sponse to time dependent perturbations at frequencies 
lu < lo c . As we have said, lu c drops rapidly upon cool- 
ing. If one is intent on observing equilibrium response 
at some frequency range, one must prepare the sample 
by cooling it more slowly than u) c . Conversely, for any 
given cooling rate, no matter how slow, the liquid will 
fall out of equilibrium on all time scales and the sample 
will appear to be mechanically solid. We say the liquid 
has undergone the glass transition. (The corresponding 
uo c usually ranges between 10 2 and 10 5 sec, depending on 
the experimenter's patience.) The liquid just below the 
glass transition temperature T g is only subtly different 
from the liquid just above T g . Structurally, first of all, 
the two are nearly identical. Even dynamically, both can 
flow, although the T-dependences of the corresponding 
transport coefficients are distinct in the two forms of the 
"equilibrium " supercooled liquid and the non equilibrium 
glassy state llLubchenko and Wolvnesi l2004|) . The resid- 
ual dynamics below T g is referred to as "aging" . Aging 
is at least as slow as the motions just above T g , but can 
be much slower when the sample is studied well below 
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T g . This requires a greater amount of the experimenter's 
patience in studying system properties than even needed 
for sample preparation. Finally, when the sample falls 
out of equilibrium at T g , a jump in the heat capacity 
is measured by differential calorimetry, thus resembling, 
crudely, a phase transition. 

The dramatic slowing down of molecular motions is 
explicitly seen in a vast area of different probes of liq- 
uid local structures. Slow motion is evident in viscosity, 
dielectric relaxation, frequency dependent ionic conduc- 
tance, as well as in the speed of crystallization itself. 
In all cases, the temperature dependence of the generic 
relaxation time obeys to a reasonable, but not perfect 
approximation the empirical Vogel-Fulcher law: 



Trlxn OC e 



DT /(T-T ) 



(1) 



For a review, see i|Angell et all . l2000t iBohmer et all 
1993). A specific example of a r(T) dependence is shown 
in the l.h.s. panel of Fig[7| In the expression above, 
To is a material dependent temperature at which the re- 
laxation times would presumably diverge, if the experi- 
menter had the patience to equilibrate the liquid at the 
corresponding temperatures. Needless to say, measure- 
ments of equilibrium dynamics near To are essentially 
nonexistent. The coefficient D is often called "fragility", 
with larger values of D corresponding to "stronger" sub- 
stances, while smaller values are associated with "frag- 
ile" liquids. This terminology apparently refer s to the 
degre e of covalent networking in the material ijAngell 
1985), a qualitative t rend later rational i zed by a density- 
functional study of l|Hall and Wolvnesl 12003). Fragility 
appears to correlate w ith the Poisson ratio, at lea st for 
non-polymeric glasses l|Novikov and Sokolovl l2004j) . At 
any rate, the value of coefficient D is directly related to 
what gla ssblowers refer t o as "short glasses" and "long 
glasses", ()P f aendeit 1 1 9 9 6j) : (molten) glass can be worked 
or shaped in the range of viscosities 10 4 — 10 9 Poise. If 
the corresponding temperature range is short, the glass 
is called "short", and vice versa for the "long" glass. The 
former and the latter obviously correspond to a small and 
large value of the parameter D respectively. 

The non-equilibrium character of a supercooled liquid 
is exhibited in the entropy of the liquid which is consider- 
ably larger at T g than that of the corresponding crystall 
at this temperature. This additional entropy corresponds 
to all the molecular translations, that would have other- 
wise frozen out at crystallization. In crystallization, this 
would appear as the latent heat of the liquid-to-crystal 
transition. In a supercooled liquid, the molecular struc- 
ture is dense enough to define a lattice locally. Vibra- 
tions around lattice sites are small. The excess entropy 
associated with the locations of these lattice sites has 
traditionally been designated as the "configurational" en- 
tropy. This excess entropy, s c , is temperature dependent. 
It refers to all possible liquid configurations that could be 
surveyed by the liquid if we wait long enough for molec- 
ular translations to occur. Experimentally, we determine 
the configurational entropy by relying on the third law of 



thermodynamics. Using the third law, we know the total 
entropy of the liquid at T m by integrating the crystal's 
heat capacity (over T) and adding the entropy of melting. 
Now for the supercooled liquid, we integrate the heat ca- 
pacity difference between the liquid and the crystal. To 
do this we, of course, assume the vibrational entropies 
of the ordered and aperiodic lattices are close. The heat 
capacity measured by differential calorimetry above the 
glass transition depends on the rate of the configurational 
and vibrational entropy decrease with temperature right 
above T g . Below T g the structure of the liquid remains 
the same as of the moment of vitrification, apart from 
some (normally insignificant) aging. The vibrational en- 
tropy decreases as it did above T g , but there is no com- 
ponent from configurational change. Thus one observes 
a non-zero heat capacity jump at T g . Above T g , the 
s c decreases and the density increases with lowering the 
temperature. This is expected because there are fewer 
ways to mutually arrange the molecules at higher densi- 
ties. When ex trapolated past T g , as was done by Simon 
and notably by Kauzmann in his review 



l KauzmannU l948'l. the configurational entropy vanishes 
at a temperature Tk, which is securely above the abso- 
lute zero. This suggests that only a non-extensive num- 
ber of low energy aperiodic, liquid arrangements could 
be found at Tk and the entropy of the liquid is thus 
equal to the corresponding crystal (correcting for differ- 
ences in their vibrational spectrum). This phenomenon 
is sometimes referred to as the "entropy crisis" , which, 
again, would presumably occur only under completely 
equilibrium cooling. Such an entropy crisis strictly oc- 
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configurational entropy at Tk 
l|Eastwood and Wolvnesl 12002 
crthcless, Tk is a useful fiducial point for the analysis. 
None of the results of the present theory in the experi- 
mentally accessible regime depend on the configurational 
entropy truly vanishing at any point. As we shall see, 
the configurational entropy is macroscopic but decreases 
with temperature. s c is typically ~ .8ks per movable 
unit at the conventional glass transition temperature cor- 
responding to cooling rate of inverse hour and decreases 
at a rate proportional to Ac p /T g . For simplicity, we will 
assume s c extrapolates so as to sc ale linearly with the 
proxi mity to the entropy crisis (see ijRichert and Angelll 
E223)): s c = Ac p {T-T K )/T K . 

Before our formal discussion, let us make several qual- 
itative statements about molecular transport above T g . 
The motions of a supercooled liquid are much slower than 
the local vibrations. The potential felt by an individual 
molecule comforms to a local "cage" . This local "cage" 
is formed by the neighboring molecules, of course. In 
order to translate irreversibly a given molecule, as op- 
posed to vibrating about the current position, the cage 
must be destroyed. In other words, a number of sur- 



7 



rounding molecues must be translated as well. Upon 
lowering the temperature, the density increases and s c 
decreases, therefore fewer alternative states are available 
to any given group of molecules. Thus it is clear that 
conforming the liquid to an arbitrary translation of a 
given molecular unit will require readjusting the posi- 
tions of more and more surrounding molecules at the 
same time. This leads to a larger cooperative region 
size, leading in turn to higher barriers for relaxation 
processes and higher viscosity. At a crude level, this 
pictu re underlies the arguments from ijAdam and Gibbsl 
119651) . but those arguments fail to relate the size of the 
moving regions to the energy landscape itself. In con- 
trast, the Random First Order Transition (RFOT ) theory 
l|Kirkpatrick et oi.lll989tlXia and Wolvnesl 120(1(1) explic- 
itly shows how these reconfigurational motions occur and 
thus establishes intrinsic connection between the kinetic 
properties and the thermodyn amics of supercooled liq- 
uids. Our account is based on ( Lubchcnk o and Wolvnesl 
l2004f) which also discusses the intrinsic connection be- 
tween cooperative, activated motions in the supercooled 
liquid both above and the classical aging dynamics below 
the glass transition. These arguments also pave the way 
for understanding the quantum dynamics at cryogenic 
temperatures. 

The main prerequisite of the RFOT theory is the exis- 
tence of time scale separation between vibrational ther- 
malization and equilibrating structural degrees of free- 
dom that involve crossing saddle points on the free en- 
ergy surface. This only occurs below a crossover tem- 
perature Ta which is predicted by the theory itself. The 
existence of local trapping in cages is well established by 
experiment: there is a long plateau in the time depen- 
dent struct ure factor as m easured by the inelastic neutron 
scattering l|Mezel [l991). In RFOT, such trapping was 
first established th eoretically using a density functional 
theory (DFT) in ijSingh et a/1 Il985|) : This paper shows 
there are aperiodic free energy minima by computing the 
free energy of an aperiodic variational density distribu- 
tion function: p(r) = p{r,{r t }) = ^ (f f' 2 e~ a ( r - r ^ 2 . 
The set of coordinates {r^} denotes a particular aperiodic 
lattice. The typical lattice spacing is a. A zero value for 
the parameter a would correspond to a completely de- 
localized, uniform liquid state, such as just below the 
liquid-vapor transition, a — > oo would imply freezing 
into an infinitely rigid lattice, a can also be interpreted 
as the spring constant of an equivalent Einstein oscillator 
forcing each molecule to remain near its proper location 
in the aperiodic lattice. F(a) develops a metastable mini- 
mum, at non-zero a — ao ^ 0, only below some tempera- 
ture Ta ■ This minimum has higher free energy than than 
the lowest minimum at a = (see Fig|2J). In the mean 
field limit, the appearance of such minimum would lead a 
lattice stiffness and would represent a state with a diver- 
gent viscosity. This localization transition and the viscos- 
ity catastrophe of mode-couplin g theories are essentially 
identic al as was established in ijKirkpatrick and Wolvnesl 
Il987a|h A single such high lying free energy minimum 



F(a) 




FIG. 2 This is a schematic of the free energy density of 
an aperiodic lattice as a function of the effective Einstein 
oscillator force constant a (a is also an inverse square of 
the localization length used as input in the density func- 
tional of the liquid. Specifically, the curves shown charac- 
terize the system near the dynamical transition at Ta, when 
a secondary, metastable minimum in f(a) begins to appear 
as the temperature is lowered. This figure is taken from 
llLubchenko and Wolvnesl l2003bft . 



would be thcrmodynamically irrelevant, but one must 
recall that this F(a) is computed for a single, particu- 
lar aperiodic lattice, which is actually only one of many 
possibilities. Taking into acount the thermodynamically 
large number of alternative aperiodic packings increases 
the entropy of the (set of) localized, aperiodic state(s) 
and thus lowers the metastable free energy minima just 
the right amount to make them competitive with the 
mean-field uniform, delocalized state. The correspon- 
dence between the free energy difference in mean field 
theory and the configuration al entropy was rigorously 
shown for the Potts Glass by ijKirkpatrick and Wolvnesl 
1987b) who argued such systems have similar symmetry 
properties to structural glasses. For structural glasses 
this correspondnce may also be shown more forma lly us- 
ing a replica formalism IjMezard and Par isi. 1999). The 
localization transition at Ta is accompanied by a discon- 
tinuous change in the order parameter a. This is why 
the transition is called "Random First Order" . Although 
there is a discontinuity in a, the actual structure in which 
the system freezes is chosen at random out of a multitude 
of possibilities (given by the configurational entropy) At 
the same time, such an ordered phase will persist only 
for finite times, therefore this is a true transition only 
for high-frequency motions, comparable at first to the 
vibrational time scale. This transition at Ta only signi- 
fies a soft cross-over, as far as the whole dynamical range 
is concerned. We emphasize, there are many different 
"phases" below Ta, all of which are random packings. 
The number of random packings, thermally available to 
a region of size TV, e Sc , decreases gradually with tem- 
perature. (This corresponds to gradual freezing out the 
translational degrees of freedom with lowering the tem- 
perature, as signified by the decreasing uj c .) Because the 
decrease is gradual, the random first order transition does 
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not exhibit a latent heat. In a finite range system, differ- 
ent minima can interconnect by barrier crossing. (Such 
barriers would be infinite in mean field.) Even though the 
transition at Ta is a crossover, the temperature Ta itself 
is a useful parameter characterizing material properties. 

The resulting time scale separation at and below Ta 
has two important consequences. First, one may per- 
form canonical averaging over the vibrations within a 
particular structural state. This gives a free energy of 
a particular structural state: $ = E — TS v ihr, where 
S'vibr is the vibrational entropy. Note the vibrations are 
not necessarily harmonic. To define <£, all that matters 
is that the local vibrations equilibrate much faster than 
the structural degrees of freedom. As a consequence, $ 
can be termed the bulk, microcanonical energy of a given 
structural state. To any value of this energy one may 
associate a bulk, microcanonical entropy S c ($) counting 
states with similar contributions from energy and vibra- 
trional entropy; both $ and 5 C ($) scale linearly with the 
size of the system. One may thus to work with mor- 
phologically distinct, globally defined aperiodic phases 
without actually specifying their precise molecular con- 
stitution, so long as we know their spectrum, i.e. their 
number as a function of the microcanonical free energy. 
These statistics are directly measurable by calorimetry 
just as in our discussion of the Kauzmann paradox. 

Having established the transitory existence of a global 
aperiodic structure, we may next enquire into how molec- 
ular motions allow the system to escape such a phase 2 . 
This occurs by replacing locally one part of the aperiodic 
packing by a different local packing. This will be an ac- 
tivated event. The RFOT theory allows one to compute 
the mean activation barrier and its distribution. Also, 
the theory determines critical region size and the spatial 
extent of the excitations corresponding to the cooperative 
rearrangement. The magnitude of an individual molecu- 
lar displacement during the transition is determined by 
a. To estimate the activation free energy, let us make 
the following construct. Considering a library of possi- 
ble local aperiodic arrangements at a particular location, 
as illustrated in Fig|3J This local library of states can be 
constructed based on the existence of the global library of 
states introduced earlier that we described by the energy 
variable <!> and the corresponding entropy S c (&) reflect- 
ing the spectrum. Clearly, the energy density exp[S' c ($)] 
is extremely high and grows rapidly with $. We might 
perform a full survey of local states by mentally carving 
out a small region of size N, while freezing in place the 
lattice sites surrounding the region. One can then heat 
the local region and then allow that region to equilibrate. 
Unless the new local arrangement is exactly the same 
as the original one, its energy will likely be significantly 



2 Of course, the issue of producing the aperiodic state in the 
laboratory would also involve estimating whether corresponding 
quenching rates can be experimentally achieved. 
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FIG. 3 This figure is taken from llLubchenko and Wolvnest 
120041) . In the upper panel on the left a global configuration 
is shown, chosen out of a global energy landscape. A region 
of N = 5 particles in this configuration is rearranged in the 
center illustration. The original particle positions are indi- 
cated with dashed lines. A larger rearranged region involving 
N = 7 particles is connected dynamically to these states and 
is shown on the right. In the lower panel, the left most fig- 
ure shows the huge density of states that is possible initially. 
The density of states found in the local library originating 
from a given initial state with 5 particles being allowed to 
move locally is shown in the second diagram. These ener- 
gies are generally higher than the original state owing to the 
mismatch at the borders. The larger density of states where 
7 particles are allowed to move is shown in the right most 
part of this panel. As the library grows in size, the states 
as a whole are still found at higher energies but the width of 
the distribution grows. Eventually with growing N, a state 
within thermal reach of the initial state will be found. At this 
value of N* we expect a region to be able to equilibrate. 



higher: A local substitution statistically must cost free 
energy, stemming from a structural mismatch between 
two randomly chosen aperiodic packings of a given en- 
ergy <F This mismatch energy corresponds to the usual 
surface energy, such as that between two different crystal 
forms or at a liquid-crystal interface. The free energy cost 
of locally replacing the initial phase (labelled as "in" ) by 
another phase, call it j, can therefore be written as 



R) 



Jib 



(R) = ^ j (N)-^ m (N)+T j 



(2) 



where T ' is the mismatch energy and R is the location 
of the local region. As before, the capital <£> denotes the 
bulk energy, corresponding to a distinct aperiodic pack- 
ing, with the vibrational entropy already included. To 
compute the likelihood of such a local rearrangement, 
substitute for the specific surface energy Tji n its average 
value which should scale with size: 7JV 1 . 7 depends on 
the material and on temperature. Naively, the usual sur- 
face energy scaling is N( d ~ 1 ^ d , expected in d dimensions. 
One can argue however that x will actually turn out to 
equal 1/2. Such a surfac e tension reno rmalization was 
first conceived by Villain ijVillainl Il985j) , in the context 
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of the random field Ising model (RFIM). In RFIM, the 
Ising spins, in addition to their coupling, are subjected 
to a random static magnetic field obeying certain fluctu- 
ation statistics. A flat interface, or domain wall, between 
spin-up and spin-down domains will distort so as to con- 
form to the local variation of the field. An RG argument 
incarnating this distortion on a hierarchy of length scales 
yields a scale dependent renormalization of the surface 
tension . givin g a surface free energy exponent x = 1/2 
ilVillainl ll985). The structure-structure interface in a su- 
percooled liquid resembles the RFIM, owing to the fluctu- 
ations of local energies of the various aperiodic packings. 
The statistics of thes e fl uctuating local energies require 
that 6&(8N) ~ QqVSN, where $ is <W-independent, 
echoing the fluctuation st atistics of the frozen ra ndom 
field of the RFIM. Thus, as ijKirkpatrick et al\E98^ sug- 
gest, the originally thin flat interface will become diffuse 
yielding x = 1/2. In the liquid vivid interpre- 

tation of the surface energy renormalization is possible: 
Since the interface is distorted down to the smallest scale 
(allowed by the material's discretness) , the region occu- 
pied by the now diffuse wall is neither of the two original 
structures it separates. Instead it may be interpreted 
as accommodating other structures. These intermediate 
packings interpolate structurally two randomly chosen, 
and thus a priori energetically disagreeable packings. In 
other words, the original thin interface separating two 
given packings, is "wetted" by other packings thus low- 
ering the overall interface energy. As we shall see, real 
liquids have only modest size regions of rearrangement, so 
it is hard to argue about the exact value of the exponent. 
Nevertheless, we note two felicitous observations: With 
x = 1/2, the usual scaling argument will give precisely a 
discontinuity in Ac p at any ideal transition, to be seen 
at Tk- Also, while the RFIM itself remains the subject 
of discussion, Villain's argument does give a length scale 
exponent agreeing with the majority of experiments a nd 
numerical studies ( Be langeit Il998t iNattermannl Il998'l . 

The role of the interface mismatch energy in the recon- 
figuration process can be beneficially understood from a 
statistical point of view, as illustrated in FigOU It costs 
free energy to reconfigure a small number N of molecules 
because considering a small region severely limits the 
number of available liquid configurations. The interface 
energy grows with AT, however the available density of 
states, too, grows with N, both in terms of its abso- 
lute value and the distribution's width. At some size N* , 
that will be computed shortly, all relevant liquid states 
become available. The rate of escape of a group of N 
molecules to another structural state can be determined 
by a canonical type sum accounting for the multiplicity 
of the final states at energy cf>j : 

- r mi 1 cl . o e Sc( *^ )/fcB e" ( ^"~ ' ! " b:)/ ' £BT . (3) 

In the second step, a steepest descent evaluation is made 
where <f> c(l maximizes the integrand, is some constant 



of units energy that reflects the local curvatures of the 
energy landscape. The quantities (f) l - ,b and <&j are related 
through Eq. (0) . The time scale r m i cr o is the time scale of 
a molecular scale non-activated process, typically of the 
order a picosecond. The value </> e q that maximizes the 
integrand above will be the internal (equilibrium) free 
energy characteristic of the system at the ambient (i.e. 
vibrational) temperature T . In other words, the greatest 
kinetically accessibility of a state, as embodied in the op- 
timization in Eq.®, implies that the state will be most 
frequently visited by the system, therefore it must be the 
equilibrium state. The integration in Eq.JJJl is similar 
to a canonical sum; yet it is different in an important 
way: The summation in Eq.® is far more general than 
the usual expression for the partition function because 
when relaxation times are continuously distributed, one 
must explicitly weigh the contribution of a state (to the 
canonical sum) by its kinetic accessibility. The latter, 
in general, will depend on the spatial extent of the exci- 
tation corresponding to a transition between two states; 
in this regard, the integration variable <f>j is, in a sense, 
a local microcanonical energy. Consequently, the energy 
4> cq corresponds to a canonical energy. Yet, (frj and 4> eq 
would strictly become a microcanonical and canonical 
energy, in their conventional sense, only in the large N 
limit, when the boundary effects are small. In contrast, 
the very thermodynamic relevance of the glassy state is 
due to the locality of the landscape and non-smallness of 
the surface term. Finally, since the bulk entropy S c (& cq ) 
corresponds to the equilibrium energy, it will be given by 
the equilibrium configurational entropy S^T), measured 
by calorimetry. Thus given oq , one can compute the 
value of the typical escape rate to a structure where N 
particles have moved. This gives: 

HN) = r m i ro exp {S C (AT,T) - ^-^} ■ (4) 

The number of particles that must be moved for com- 
plete equilibration is determined by the minimum of this 
expression over TV. We thus determine an activation free 
energy profile 

F* (AT) = -TS C (N,T) 

= ^ eq (N)-^ in (N)+ 1 VN-TS c (N,T),(5) 

where we used Eq.(j3J| in the second equality. The maxi- 
mum of the F(N) curve defines the bottleneck location. 
This equation is suitable for finding the rate of structural 
rearrangement both in the equilibrated supercooled liq- 
uid (before it crystallizes!) and in the nonequilibrium 
glass, which ages below T g . 

Let us first consider equilibrium liquid rearrangements. 
In this case typically <f> cq = <f>i n , apart from fluctuations. 
Thus one arrives at the following simple expression, 

F(N) — jVN — Ts c N, (6) 

where we have used the thermodynamic scaling of the 
configurational entropy, S C (N) = s c N. In the super- 
cooled equilibrated liquid, molecular transport is driven 
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by only the multiplicity of mutual molecular arrange- 
ments. For this reason, the reconfigurations following 
the activation profile from Eq.® have been called "en- 
tropic droplets" . The graph of the function in Eq. JBJ is 
shown in Fig0] The transition state configuration will 



F(N) 




FIG. 4 The droplet growth free energy profile from Eq.© is 
shown. 

satisfy dF/dN = 0, corresponding to an unstable sad- 
dle point of this free energy. This gives for fixed 7 a 
rearranging region size TV* that grows as s c diminishes: 
N* = ( r -f/2s c T) 2 . The resulting barrier also scales in- 
versely proportionally to s c : 



F* = 



7 



4s,T 



(7) 



An inverse scaling of the barrier with the configura- 
tional entropy was arri ved at by Adam and Gibbs 
(lAdam and CxibbsL I1965D in a different (and inequiva- 
lent) way. Notice if 7, as function of temperature, is 
smooth around Tk and s c is described by the linear law 
s c oc (T — Tk), the resulting activation barrier is ex- 
actly of the Vogel-Fulcher law form J[J, which, as we 
have said, fits data well. Many arguments can lead to in- 
creasing relaxation times at low temperatures and with 
enough adjustable parameters, can fit data. What is dif- 
ferent about the RFOT theory is that it establishes an 
intrinsic link between the rate law and the entropy cri- 
sis. In addition, if the entropy of the equilibrated fluid 
can be estimated, the density functional theory allows 
the vibrational entropy and thus, by substraction, the 
configurational entropy to be determined. Therefore Tk 
can be estimated from the microscopic force laws. This 
has be en done for simple soft s pheres by Mezard and 
Parisi ijMezard and PariS . Il999l). giving reasonab le re- 
sults. Hall and Wolvnes piall and Wolvnel F2003) have 
also calculated To and Ta for a simplified model of a net- 
work fluid. Their study is consistent with known chemi- 
cal trends for Ta and Tk as the network becomes more 
thoroughly crosslinked. 

The idea of the configurational entropy itself driving 
liquid rearrangements still appears to generate some con- 
fusion. One possible reason for this is that s c is totally 
unambiguously defined only in the mean field limit. In 
the latter limit, rearrangements are infinite so dynamics 
driven by s c do not arise. This is a good place to empha- 
size that the RFOT theory is not mean-field! Only the 
local landscape, within an entropic droplet, is actually 



well described by a mean-field, Random Energy Model 
like approximation. We took advantage of this in ex- 
tracting the energy spectrum of low energy structural 
excit ations in a frozen glass ijLubchenko and Wolvne4 
l200lj) . as explained in detail in the following Section. We 
wis h to point the reader to the recent elegant treatment 
of l|Bouchaud and Biroli l2004[) re-analyzing the RFOT 
conclusions for rearrangements in an equilibrated fluid 
from th e viewpo i nt of Derrida's Random Energy Model 
(REM) llDerridalll98l . 

Now, calculations of Ta and Tk are plagued by the 
usual difficulties of liquid state structure theory and the 
accuracy of approximations some of which are hard to 
control. Still, even in the face of such approximations, 
such microscopic considerations lead us to expect a uni- 
versal value of j/Tg at T g as we shall discuss below. 

The RFOT theory allows the coefficient 7 in the mis- 
match energy to be estimated from a microscopic ar- 
gument. It turns out to be proportional to Tk and 
to depend logarithmically on the inverse square of the 
so called Lindemann ratio. Early in the 20th century, 
Lindcmann argued that the thermal fluctuations of an 
atom's position could only be a fraction of the lattice 
spacing a in a solid, if the pac king is to be mechan- 
ically stable l)Lindemannl Il910[) . Since the threshold 
value of the vibrational amplitude of an atom in the 
lattice is finite, the transition in which the lattice dis- 
integrates must be first order. For crystals, the Lin- 
demann ratio of this threshold displacement to the 
lattice spacing is about 1/10. For amorphous materi- 
als, the d^/a ratio can be obtained from the plateau 
in the self correlation functions me asured by neutron 
scattering experiments (Mezei, 1991). Again, this ratio 
turns out to be approximately one-tenth (universally!). 
This number is reproduced in several microscopic calcu- 
lations consistent with the RFOT theory, such as the self- 
consistent phonon theory and density functiona l theories 
l|Singh et all Il985t IStoessel and Wolvneslll984h. and dy- 



namical mode coupling theory ( Bcngtzelius et 



„... , TS 

iKirkpatrick and Wolvnelll987albl:IWolvneslll992l) . with 
modest quantitative variations. The meaning of a ~ 0.1 
as a mechanical stability criterion h as been also corrob- 
orate d within the replica formalism ijMezard and Parisil 

1999) . In terms of the DFT calculation dicussed earlier, 
corresponds with the metastable minimum that the 

free energy F(a) develops below the dynamical transi- 
tion temperature Ta (see Fig|2} ■ It has a relatively weak 
temperature dependence. The logarithmic scaling of the 
surface tension coefficient with the Lin demann length fol- 
lows from a detailed calculation by ijXia and Wolvnesl 

2000) , but can be rationalized in a simple way: Below 
Ta, motions span only the length d^, while in the liquids, 
they can move a distance a before losing their identity 
with a neighboring molecule. The entropy of the "caged" 
fluid is less and thus the free energy cost of confining a 
molecule within length di,, as opposed to a, can be as- 
sessed by recalling the free energy expression for an ideal 

monatomic gas: — / — ^/c^Tm 



' eV 
v N 



i 2 / 3 mT 



written 



11 



deliberately here so as to have a length scale squared in 
the logarithm. 

7 is proportional to and only logarithmically de- 
pends on a nearly universal quantity, the Lindemann ra- 
tio. If T g is near Tk, i.e. for slow quenches, j/Tg is 
thus nearly material independent and calculable: 7 = 
I y/STrksTg \a.(a.La 2 /ire). Quantifying the mismatch en- 
ergy this specifically leads to many predictions about the 
dynamics near T g , for a range of substances. First, the 
coefficient in the Vogel-Fulcher law D is predicted to fol- 
low from the measured thermodynamics. Using the 7 
value above, we find not only the VF dependence of the 
relaxation times on the temperature, oc e DT °^ T ~ T °\ but 
also a remarkably simple formul a relating D and the hea t 
capacity jump: D = 32.R/Ac p ijXia and Wolvnesl. [2000 b 
The coefficient 32. is nearly universal and, as we see, fol- 
lows numerically from the microscopic theory since the 
universal value of the Lindemann ratio enters only loga- 
rithmically in the localization entropy cost. The numeri- 
cal relation between D and Ac p from this simple explicit 
calculation is in rather remarkable agreement with exper- 
iment. In FigEI we plot the so called fragility index m, 
as computed from calorimetry and extracted from direct 
relaxation measurements, m is proportional to the slope 
of the logr vs. 1/T relation at T g and thus is directly 
related to D if the VF law is valid. (D values in the lit- 
erature are obtained from global fits of logr vs. 1/T and 
depend somewhat on the fitting procedure.) 

Two other remarkable universalities emerge from the 
value of 7. First, at a reference laboratory time scale of 1 
hr ~ 10 17 To we have a universal value of s c ~ 0.8fcs. This 
implies s c (T g ) / s c {T m ) ~ 0.7, where s c (T m ) is, of course, 
also the fusion entropy. This relation is independent of 
question of what is the moving subunit. The relation 
holds very well. A second important universal feature 
emerges from the universal value of 7 /T g : the cooperative 
size at T g is nearly universal. 

Let us now consider in greater detail the pattern of 
cooperative structural rearrangements in a supercooled 
liquid. These turn out to presage the existence of the 
residual degrees of freedom in a glass below T g . Within 
a period of time shorter than the typical relaxation time 
t, the molecular motions within regions of size £ 3 will be 
highly correlated and, at the same time, approximately 
decoupled from the surrounding. That is, the liquid is 
broken up in to a (flickering) mosaic pattern of coop- 
erative regions. This mosaic structure is directly mani- 
fested in the dynamical heterogeneity recently observed 
in supercooled liquids usin g single molecule experiments 
IjRusse l and IsraelofJ l2000f) , nonlinear relaxation experi- 
ments (ISilesc ul. Il999ft and non-linear NMR experiments 
ijTracht et all Il998j) . (These experimental tools became 
available only a decade after the RFOT theory was first 
formulated.) The size of a typical mosaic cell is found 
from the thermodynamic condition F(N*) = 0. Unlike 
the regular nucleation of one distinct phase within an- 
other (as in crystal growth in the liquid), by crossing 
the barrier from Eq.© the local region arrives at a sta- 
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FIG. 5 The horisontal axis shows the value of fragility as 
computed from the thermodynamics by the RFOT theory, 
and the vertical axis contains the fragility directly measured 
in kinetics exoeriments. Here m is the so called fragility 
index, defined according to m = T[dlog 10 t(T) / d(l/T)]. 
m is somewhat more useful than the fragility D, be- 
cause deviations from the strict Vogel-Fulcher law, r = 
Toe DTK K t ~ Tk ' , are often observed, see text, m essentially 
gives the apparent activation energy of relaxations at T g , 
in units of T g , it is roughly an inverse of D. This fig- 
ure is taken from iStevenson and Wolvnest 12004) . In eval- 
uating m theoretically, one needs to know the size of the 
moving unit, or "bead", in each particular liquid. The lat- 
ter c an be estimated using the entr opy loss at crystalliza- 
tion iLubchenko and Wolvnesll2003b() . resulting in m^or oc 
Afl^^ffi) oc Ac i^Tg)Tg ^ j n v - iew Q £ -j-jjg near universality of 

s c (T g ) (see text). 



tistically similar but an alternative solution of the free 
energy functional, thus that solution still represents a 
typical liquid state! An informal analogy here is that 
distinct low energy dense local liquid packings are like 
the fingerprints of different individuals - different in de- 
tail, yet generically equivalent liquid states. Since we 
have agreed that F = is the liquid equilibrium free 
energy at this temperature (the crystalline state is as- 
sumed to be hidden behind a high enough crystal nu- 
cleation barrier), the condition F(N*) = specifies the 
size of region to which an arbitrary liquid configuration 
is available. Therefore, a region of size N* is able to 
reconfigure on the experimental time scale characterized 
by F*. In terms of physical length, F(N*) = implies 

2/3 
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£ = 

(N* ~ 190). The critical radius r* at T g is a multiple 
of £. Droplets of size N > N* are thermodynamically 
unstable and will break up into smaller droplets, in con- 
trast to what prescribed by F(N), if used naively beyond 
size N* . This is because V = and V = V* repre- 
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sent thermodynamically equivalent states of the liquid 
in which every packing typical of the temperature T is 
accessible to the liquid on the experimental time scale, 
as already mentioned. In view of this "symmetry" be- 
tween points N — and N* , it may seem somewhat odd 
that F(N) profile is not symmetric about N*. Droplet 
size N, as a one dimensional order parameter, is not a 
complete description. The profile F(N) is a projection 
onto a single coordinate of a transition that must be de- 
scribed by e s ° order parameters - the effective number 
of distinct aperiodic packings explored by the liquid. At 
the point N* , the free energy functional actually has a 
minimum as a function of the (multi-component) order 
parameter. A more detailed discussion of t his can be 
found in Ref. ljLubchenko and Wolvnesl l20"0 3bh where we 
compute the softening of the barrier F$ near Ta due to 
order parameter magnitude fluctuations that are impor- 
tant near Ta- 

We thus see that the length scale of the mosaic and 
number density of the mosaic domain walls is determined 
by the competition between the energy cost of a domain 
wall and the entropic advantage of using the large num- 
ber of configurations. We emphasize again, the relative 
domain size £/a depends only on the logarithms of the 
relaxation rate and the Lindemann ratio, nearly univer- 
sal parameters themselves, and is therefore the same for 
different substances. This high temperature phenomenon 
of universality at T g has direct consequences for the uni- 
versality of the ultra-low temperature glassy anomalies. 

We have seen that the cooperative region, which 
represents a nominal dynamical unit of liquid, is of 
rather modest size, resulting in observable fluctuation ef- 
fects. Xia and Wolynes computed the rel axation bar- 
rier distribution ijXia and Wolvnesl l2001a|) . The con- 
figurational entropy must fluctuate, with the variance 
given by the usual e x pressi on: ((5S C ) 2 ) = C p oc 1/D 
IlLandau and Lifshitd. [1980). The barrier height for a 
particular region is directly related to the local density of 
states, and hence to the configurational entropy itself by 
Eq.@), F* oc l/s c . As a result, the barrier distribution 
width must correlate with the fragility. A gaussian ap- 
proximation leads to a sim ple formula 5F$ /F$ = l/2\fD 
(|Xia and Wolvnesl 12001 a)) . There are also calculable de- 
viations from gaussianity. The barrier distribution gives 
rise to non-exponentiality of relaxations. These are well 
fitted by a stretched exponential e - ^/' '' 3 . The measured 
f3 correlates with the fragility, in good agreement with 
the theory, see FigEl 

We have so far presented a simplified picture of acti- 
vated relaxation in liquids, which is more accurate at 
temperatures close to Tk, and thus sufficiently lower 
than Ta - the temperature at which activated processes 
emerge. The transition at Ta where metastable min- 
ima emerge, along with a mosaic structure with inter- 
mediate tense regions, i.e. domain walls, is in many 
respects similar to a spinodal for an ordinary first or- 
der transition, except that the number of alternative 
phases is very large (e ScN for a region of size N). The 
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FIG. 6 Shown is the correlation between the liquid's fragility 
and the exponent /3 of the stretched exponential relaxations, 
as predicted by the RFOT theory, superimposed on the mea- 
sured values in many liquids taken from the compilation of 
(iBohmer et all 1 1993) . Fhe dashed line assumed a simple 
gaussian distribution with the width mentioned in the text. 
The solid line takes into account the existence of the highest 
barrier by replacing the barrier distribution to the right of 
the most probable value by a narrow peak of the same area; 
the peak is lo cated at that most proba ble value. Fhis figure 
is taken from (IXia and Wolvnesl l2001al) . 



proper treatment of this transition must include fluctu- 
ations of the order parameter and consequent softening 
of the droplet surface tension at temperatures close to 
Ta- As a result of this, closer to Ta the structural re- 
laxation barriers are lowered from what would be ex- 
pected extrapolating from near Tk - this gives devia- 
tions from the VF law. The corresponding length scales 
r* and £ also should be smaller than would be predicted 
by the "vanilla" , T^-asymptotic version of the RFOT 
theory. These barrie r "softening" effects were quanti - 
tatively estimated in l)Lubchenko and Wolvnesl l2003b(l . 
They demonstrated that softening effects do vary be- 
tween different substances and are more pronounced for 
fragile liquids. As a result, the value of the configura- 
tional entropy at T g , as predicted by the RFOT theory 
with softening varies somewhat, within a factor of two or 
so among different substances. This is in contrast to the 
universal s c (T g ) = .82 of the vanilla RFOT. Nonetheless, 
the value of £ at T g is much less sensitive and seems to 
be always within 5% of the simple estimate above. This 
is shown in the r.h.s. panel of Fig[7| Understanding of 
the softening effect has allowed us to compute the acti- 
vation barrier for liquid rearrangements in the full tem- 
perature range, including the high T part near Ta, where 
the barriers become low, and the transport is dominated 
by activationless, collisional phenomena. Consistent with 
this predicted softening, the T dependence of relaxation 
times, r = r micro e' ) ' 2 / 4T2s <=(' r \ as predicted by the RFOT 
(see Eq.©), fits well the experimental dependences in 
the low frequency range, but underestimates the viscosity 
near boiling. After softening is included, one can com- 
pute the activation component of the molecular trans- 
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FIG. 7 Experimental d ata (symbols) for TNB's viscosity 
jPlazek and MariH [l968), superi mposed on the results of 
the fit ting procedure (line) from llLubchenko and Wolvnesl 
l2003bfl are shown. Ta is shown by a tickmark. The temper- 
ature T C r signifies a cross-over from activated to collisional 
viscosity, dominant at the lower and higher temperatures re- 
spectively (see text). The temperature is varied between the 
boiling point and the glass transition. The r.h.s. pane de- 
picts the temperature dependence of the length scales of co- 
operative motions in the liquid. The thick solid and dashed 
lines are and £ respectively. This figure is taken from 
llLubchenko and Wolvnesl l2003bft ■ 



port, with the temperature Ta as a fitting parameter of 
the theory ijLubchenko and Wolvnest l2003hj) . Fitting the 
viscosity was performed using the following obvious con- 
straints: (a) at low temperatures, the order parameter 
a fluctuations are negligible, the barriers are fully estab- 
lished and high, and the transport is thus fully activated; 
(b) near boiling, the barrier vanishes, and the viscosity 
(known to be around a centipoise for all liquids) gives the 
value of Tmicro- The fit, shown in Fig0 demonstrates that 
of the 16-17 orders of the total dynamical range, about 
three orders, on the low viscosity side, are dominated by 
collisions. The experimental and activation-only theo- 
retical curve differ from each other above a temperature 
T cr . The three order of magnitude time scale separation, 
arising internally in the theory, is indeed consistent with 
the prerequisite of the transport being fully activated 
at T cr and below. The discussion above indicates that 
samples quenched (sufficiently fast) from a temperature 
T > T cr may exhibit somewhat distinct detailed molec- 
ular motions, also implying quantitative deviations form 
the RFOT predictions. At any rate, these sample, being 
caught in a very high energy state, are expected to have 
small cooperative regions, and also be very brittle and 
in general mechanically unstable. Such rapid quenches 
would be extremely difficult to produce in a lab, because 
T cr corresponds to relaxation times of the order 10~ 8 -- 9 
sec. On the other hand, it is these ultra fast quenches, 
that must be currently employed by simulations owing 
to the limitations of computer power. We speculate that 
the thin "amorphous" films made by vapor deposition on 
a cold substrate also may sometimes correspond to such 
ultra-quenches. While one may expect a number of be- 
haviors in the bulk that are qualitatively distinct from 



what we have discussed here, various surface effects are 
likely to be important too: For one thing, such films are 
thin, have a large free surface, and strongly interact with 
the substrate. Further, there is a good reason to believe 
these films undergo loca l cracking, and spontaneous crys- 
tallization l|PerrvM200l . 

The present article deals with phenomena in glasses at 
temperatures much much lower than the temperatures 
at which the samples form. If a sample, upon vitrifi- 
cation, is cooled significantly below T g , its lattice re- 
mains practically the same as of the moment of freez- 
ing. Indeed, the typical reconfiguration barrier is at least 
ln(10 15 ) ~ 35fcsT g , as already mentioned. If, on the 
other hand, the sample in maintained at some temper- 
ature T close enough to T gi exceedingly slow structural 
relaxations take place. These attempts of the sample to 
equilibrate to a structure characteristic of temperature 
T can be detected. Achieving quantitative accuracy in 
such experiments is difficult. Consistent with the notion 
that the lattice, and the barrier distribution, freeze in 
at the glass transition, the relaxation below T g , obeys 
approximately the following temperature dependence: 



feo exp 
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AE* 



(1 - inmt) 



AE* 



(8) 

where E* is the equilibrated apparent activation en- 
ergy at T g and xnmt lies between and 1. This 
equation is part of the Nayaranaswan y-Moynihan-Too l 
(NMT) empirical descr i ption of aging l|Movnihan et all 
119761: iNaravanaswamvt Il971i iTool Il946h . The differ- 
ence in the apparent activation energy above and below 
T g , as expressed by the parameter xnmt, will depend 
on how fast the barrier itself was changing, with cool- 
ing, above T g , under "equilibrium" cooling conditions. 
Since the rate of that change depends on the fragility, 

log 10 e, one expects that 
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£nmt an d m are correlated. The RFOT based theory of 
aging in l)Lubchenko and Wolvnesl l2004|) analyzes struc- 
tural rearrangements in a non-equibrium glassy sample 
by means of Eq. JSJ , where the initial state is not equilib- 
rium, but instead corresponds to the structure frozen-in 
at T g . The predicted correlation between xnmt and m 
is very simple: m ~ 19/xnmt, and is consistent with 
experiment, see FiglHl 

For some of the comparison of theory and experiment 
it is necessary to be specific about the molecular length 
scale a (a very detailed discussion of this q uantity can be 
found in l|Lubchenko and Wolvnesl l2fl03bJ )V The molec- 
ular scale denotes the lattice spacing between molecu- 
lar units (or "beads" ) that act as idealized spherical ob- 
jects at the ideal glass transition at Tk- The determi- 
nation of a, though approximate, is rather unambigu- 
ous and can be done using the knowledge of chemistry 
to give values accurate within 15%. For example, the 
number of beads in a chain molecule, that interacts with 
the surrounding only weakly, is always close to the num- 
ber of monomers. Highly networked substances, such as 
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FIG. 8 The fragility parameter m is plotted as a function 
of the NMT nonlinearity parameter 2?nmt ■ The curve is pre- 
dicted by the RFOT theory when the temperatu re variation of 
7o is neglected. The data are taken from Ref. llAngell et all 
2000). The disagreement may reflect a breakdown of phe- 
nomenology for the history dependence of sample preparation. 
The more fragile substances consistently lie above the predic- 
tion, which has no adjustable parameters. This discrepancy 
may be due to softening effects. 



amorphous silica, present a more difficult case, because 
it is not clear how covalent the intcrmolecular bonds in 
these substances are. Since melting also involves free- 
ing up molecules, with encreased entropy, an indepen- 
dent check on the soundness of a particular bead num- 
ber assignment can be done by comparing the fusion en- 
tropy of the substance (if it exists in crystalline form) 
with the known entropy of fusion of a hard sphere liq- 
uid or Lenna rd- Jones liquid, equal to 1.16A :b and 1.68fcs 
respectively l)Hansen and McDonaldL LI976|) . Note, how- 
ever, the knowledge of the absolute value of a is not re- 
quired for most of the numerical predictions the theory 
will make in the quantum regime. 

We thus see that the RFOT theory provides a rather 
complete picture of vitrification and the microscopies of 
the molecular motions in glasses. The possibility of hav- 
ing a complete chart of allowed degrees of freedom is very 
important, because it puts strict limitations on the range 
of a priori scenarios of structural excitations that can 
take place in amorphous lattices. This will be of great 
help in the assessment of the family of strong interaction 
hypotheses mentioned in the introduction. To summa- 
rize, the present theory should apply to all amorphous 
materials produced by routine quenching, with quanti- 
tative deviations expected when the sample is partially 
crystalline. The presence and amount of crystallinity can 
be checked independently by X-ray. It is also likely that 
other classes of disordered materials, such as disordered 
crystals, will exhibit many similar traits, but of less uni- 
versal character. 



III. THE INTRINSIC EXCITATIONS OF AMORPHOUS 
SOLIDS 



A. The Origin of the Two Level Systems 



In this section we discuss how phenomena near the 
glass transition temperature, described in the previous 
subsection, dictate the existence and character of the 
quantum excitations in glasses at liquid helium tempera- 
tures and below. As mentioned earlier, a dynamical pat- 
tern of cooperative regions forms in a supercooled liquid 
below Ta- Each cooperative region is defined by the ex- 
istence of at least two distinct configurations mutually 
accessible within the time scale r, which chatacterizes 
the life-time of the local mosaic pattern. Conversely, a 
molecular transport event is made possible by rearrang- 
ing molecules within the cooperative length scale. The 
mosaic pattern "flickers" on the time scale r; this pro- 
cess slows down dramatically upon vitrification and, be- 
low T g is referred to as "aging", as it corresponds to 
(very slow) structural changes. At T g , the existent pat- 
tern of transitions (with distributed energy changes and 
reconfiguration barriers) freezes in because each cell is 
now surrounded by a rigid lattice (this is because the 
rearrangements of the neighboring domains were uncor- 
rected at T g ). Each region of the material can now ex- 
plore the phase space as prescribed by the environment 
at the time of freezing. Below T g , the mosaic is spatially 
defined by the molecular motions that were not arrested 
at T g , and is thus strictly speaking only dynamically de- 
tectable. It is true that the weaker walls will probably 
be the site of (unstable) instantaneous normal modes in 
the fluid state with imaginary frequencies. This dynam- 
ical correlation pattern does not necessarily imply any 
easily discernible spatial heterogeneity in the atomic lo- 
cations. In fact, there has been no direct evidence for 
any static type of heterogeneity of the appropriate scale 
in glasses so far, which definitely contributes to the (un- 
derappreciated) mystery of glasses 3 . But can the dy- 
namical heterogeneity be seen directly? We will claim 
later that this is done for us by thermal phonons: the 
magnitude of scattering at the plateau can only be ex- 
plained by presence of dynamical heterogeneities. The 
latter are signified by structural transitions that scatter 
the phonons inelastically. Apart from aging (which we 
will ignore in the rest of the work), a particular pattern 
of flipping regions, as frozen-in at T g , will persist down to 
the lowest temperature. The apparent size of each cell in 
this mosaic of flippable regions will depend on the obser- 
vation time. The longer this time is, the more structural 
relaxation degrees of freedom (from the high barrier tail 
of the barrier distribution) one should observe. Even- 



3 We note, however, that there have been instances of mistaking 
polycrystalline samples for truly amorphous ones. 
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tually, in fact, the glass should crystallize 4 . In order to 
estimate the number of tunneling centers that are ther- 
mally active at low temperatures we will have to find the 
size of the regions that allow for a rearrangement accom- 
panied by a small energy change and, at the same time 
with a low barrier. It may reasonably seem that typically 
such barriers for multiparticle events would be very high. 
Nevertheless, the lattice is arrested in a high energy state. 
We can thus foresee the possibility of stagewise barrier 
crossing (or tunneling) events, when the width of the bar- 
rier for each consecutive atomic movement is only a small 
fraction of a typical interatomic distance, thus rendering 
individual atomic movements nearly barrierless. This is 
as if one could define an instantaneous mode of nearly 
zero frequency, at each point along the tunneling trajec- 
tory. (Yet at no point is the motion harmonic per se!) 
The presence of such low frequency modes should be ex- 
pected given the high number of configurational states 
available to the sample as the moment of freezing, as re- 
flected in the high value of the configurational entropy 
at T g . After all, the material is unstable, both glob- 
ally and locally! (Note, the extent of bond deformation 
during an individual atomic movement is low - within 
the Lindemann length - actually affording a few "hard" 
places along the tunneling trajectory, where the "instan- 
taneous" frequencies are not necessarily low.) One may 
contrast the situation above with, say tunneling of a 
substitutional impurity in a crystal, a system which is 
indeed near its true ground state. Such tunneling would 
not contribute to the very low T thermal properties ow- 
ing to a large barrier. Also, we note that multiparticle 
barrier crossing events have bee n seen in computer sim- 
ulatio ns of amorphous systems jG^taa^md Rahman, 
1986), anticipated theore tically ijHeuer and Silbe\ill993t 



Mon and AshcroftLll97l . and recently infer red from sim 



ulatio ns of dislocation motions in copper l)Vegge et all 
120011) . 

We summarize the discussion so far by noting that 
the preceding Section has demonstrated that the possi- 
ble atomic motions in a supercooled liquids are either 
purely vibrational excitations or structural rearrange- 
ments. Any possible motions below T g , in terms of 
the classical basis set must be a subset of the motions 
above T g , although the dynamics of these events become 
quantum-mechanical at low enough temperatures. Even 
after the system is cooled to an arbitrarily low tempera- 
ture, it remains essentially in the configuration in which 
it got stuck at the glass transition. The density of di- 
rectly accessible states at that high energy configuration 
is rather high; the total density of states is, of course, ex- 
ponentially larger, but inaccessible on realistic time scale 
without other regions of the glass rearranging. Since the 



typical rearrangements near T g span about a length £ 
across, we may make the following, preliminarily conclu- 
sion: The non-equilibrium character of the glass transi- 
tion necessarily dictates the existence of intrinsic addi- 
tional non-elastic degrees of freedom in a glass, tenta- 
tively one per region of roughly size £, in addition to the 
usual vibrations of a stable lattice. The universality of £, 
in a sense, is the main clue to the cryogenic universality 
that is observed. A schematic of a cooperative region is 
shown in FigOU 
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4 Note that there are, in principle, other ways to move molecules 
in a glass, in addition to the cooperative rearrangements: for 
example by creating defects such as vacancies (the corresponding 
barriers are prohibitively high, of course) 
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FIG. 9 A schematic of a tunneling center is shown. £ is 
its typical size, dh is a typical displacement of the order 
of the Lindemann distance. The doubled circles symbolize 
the atomic positions corresponding to the alternative inter- 
nal states. The internal contour, encompassing iV* beads, 
illustrates a transition state size, to be explained later in the 
text. 

Note that showing the existence of low energy tunnel- 
ing paths is really a mathematically problem of finding 
hyper-lines, connecting two points of particular latitude 
on a high-dimensional surface, that meander within a 
certain latitude range. Visualizing high-dimensional sur- 
faces is prohibitively difficult, while the field of topology, 
at its present stage of development, is of little help. Yet, 
a completely general argument is not required here: We 
only need to consider a very small subset of all surfaces, 
such that they satisfy the (very severe) constraint on the 
liquid density distribution above T g , namely such that 
conforms to an "equilibrated" liquid at T g . Because (and 
with the help) of this constraint, it is possible to put forth 
a formal argument showing that there are indeed enough 
low energy structural transitions in a frozen glass: This 
argument will follow (albeit in the reverse order, in a 
sense) the argument from the preceding Section, where 
we found the typical trajectory for rearrangement. The 
key point of the microcanonical-like library construction 
from Section^ is that the distribution width of energies 
of a region increases with region size. A region is guaran- 
teed to have a state at some low energy, call it Eqs(N), 
as found by integration in Eq.©. Past a certain critical 
size, this energy decreases as N grows larger, giving rise 
to the existence of a resonant state at a large enough N. 
One must bear in mind, however, that Eqs(N) reflects 
the typical freezing energy. It really gives an upper bound 
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on the lowest energy level. The actual lowest energy state 
fluctuates and always lies below Eqs(N), although most 
likely not much below. Here we will look in detail the 
statistics of these energy states below the typical recon- 
figuration profile, with the aim to find the probability of 
a low energy trajectory for reconfiguring a region size N. 

We will make several preliminary, quite general notions 
that will guide us in constructing an adequate approxi- 
mation for the local statistics of the energy landscape of 
a frozen lattice. First, we give a general argument of 
the density of frozen-in excitations, valid, as we will see 
shortly, in the limit of infinitely slow aging: Since the 
atomic arrangement does not change upon freezing, the 
classical density of states of a frozen glass is that of the 
supercooled liquid at T g . Those states correspond to con- 
figurations in which the system could have frozen at T g 
and in principle can explore, provided they are thermally 
accessible and have a sufficiently low barrier separating 
them from a given configuration. Take a generic liquid 
state at T g as the reference state. Then the Boltzmann 
probability to switch to a conformation higher in energy 
by amount e is cx e~ t l T< > . That a configuration with that 
energy was one of the allowed configurations upon freez- 
ing means there must have been n(e) = Jj-e e / Tfl of them. 
The factor 1/T g arises because the energy spectrum by 
construction is continuous, while the actual local spec- 
trum is discrete and e = gives the upper bound on the 
location of the actual ground state. The latter must be 
somewhere between and — oo: J_ oQ den(e) = 1. This 
argument, however, is silent as to what the spatial char- 
acteristics of such excitations or their time scale are. 

This inverse "Boltzmann" density of states has been 
computed expl i citly i n frustrated mean-field spin systems 
l|Mezard et aZllT985j) . but is of more general nature. In- 
deed, such distributions arise universally when describing 
the statistics of the lo west energy state of a wide cla ss of 
energy distributions ijBouchaud and Mezardi Il997ft . in- 
cluding the random energy model l|Derridal |1981'). that 
will be used later on. Kinetic considerations did not 
explicitly enter our heuri stic derivation above (or, the 
mean-field es t imate s in l|Bouchaud and Mezardi Il997t 
iMezard et all IT985(l '). This is directly seen by differen- 
tiating d\ogn(e)/de = 1/Tg. Clearly, n(e) is the micro- 
canonical density of states corresponding to the trans- 
lational (liquid-like) degrees of freedom, and the system 
is assumed to be completely ergodic within that set of 
states. This corresponds to an approximation where we 
consider all degrees of freedom which are faster than a 
given time scale as very fast, and everything slower than 
that chosen time scale is regarded to be much slower than 
can be detected in the experiment. By using this same 
density n(e), as it was at T g , also at T < T g , we for- 
mally express the fact that this subset of the total density 
of states no longer thermally equilibrates but stays put 
where it was at T g - the subsystem of the translational de- 
grees of freedom has undergone an entropy crisis, a glass 
transition. Everywhere in the discussion above, we have 
been ignoring the contribution of the purely vibrational 



excitations to the total free energy. We thus assume that 
the spectrum of those elastic excitations is independent 
of precisely where on the glassy landscape the liquid is. 

We now give the argume nt, first laid out in 
ijLubchenko and Wolvne4 l200l|) , that allows one to es- 
timate the classical density of states and will also simul- 
taneously yield the size of the region where the excita- 
tion takes place. First we address the question of how 
many structural states are available to a compact frag- 
ment of lattice of size N, regardless of the barrier that 
separates those alternative states from the initial ones. 
This corresponds to the assumption of time scale sepa- 
ration mentioned just above. Within this assumption, 
the low energy limit of the spectrum must obey e E ^ Tg so 
as to give a glass transition at T g . Next, the spectrum, 
when integrated, must give e ScN for the total number of 
states available to the region. Notice further that we ex- 
pect the reconfiguring regions to be relatively small. The 
atomic motions within these small regions are directly 
coupled and so a mean-field, gaussian density of states, 
that only describes lowest order fluctuations around the 
mean, should be accurate enough. An energy density sat- 
isfying the requirements above actually corres ponds to 
the w ell known Random Energy Model (REM) l)Derrida| . 
11981ft . which also describes the pure state free energy in 
mean field frustrated spin models: 



Q, N (E) ~ exp { s c N 



[E - (NAe + -yVN)} 
2SE 2 N 



(9) 



where 5E 2 is the variance to be determined shortly. Here, 
the factor e ScJV gives the correct total number of states, 
the term j-t/N takes into account the interface energy 
cost of considering distinct atomic arrangements with the 
region. Note the fluctuations in the surface term are ex- 
pected but are automatically included in the fluctuation 
of the microcanonical energy E itself.The term NAe is a 
bulk energy necessary to account for the observed excess 
energy of the frozen structure relative to the energy of the 
ideal structure at Tk- It is easy to relate to measured 
quantities: To do this, recall that the system freezes in 
its "ground state", with energy Eqs, when its entropy 
becomes non-extensive: 



ft N (E GS ) = 1. 



(10) 



We take the energy E g of the liquid state at T g as 
the reference energy. Next note that in the absence of 
the surface energy term js/N, the lowest available en- 
ergy state is that of the liquid at [Ek ~~ E g )/N = 
- ff B K dTAc p (T) ~ ~Ac p {T g - T K ) =i ~T g s c . (The two 
latter equalities are accurate fo r T„ close to 7V. The cor- 
rections wou ld be observable l)Lubchenko and Wolvnes( . 
2003b, 2004), but small.) One immediately gets from 
Eqs.© and fifty that Ae = \f25E 2 s c -T g s c (7 = must 
be used in this estimate, but nowhere else!). Further, 
using the microcanonical d\n£lN(E)/dE\E=E GS = 
fixes the value of the variance SE 2 = 2T 2 s c . The resul- 
tant density of states is proportional to e^ E ~ E Gs)/ T 9 at 
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T g , as already shown above by a general argument. Now 
that we have determined the thermodynamical quantities 
entering Eq.©, we can find how the excess energy of an 
alternate ground state depends on the size N: 

E GS (N)=-f^N-T gSc N, (11) 

where Eqs is defined by Eq. (|10|) . Only low energy excita- 
tions will be thermally active at the lowest temperatures. 
Therefore, we are looking for excitations that are nearly 
isoenergetic with the reference state. This imposes an 
additional condition Eqs(N) = thus prescribing the 
minimal size No of a region such that has fijv(O) > 1 for 
N > No- A region of this size has at least one alter- 
native structural state at the same energy. One obtains 
from Eq.lJTTJ that N a = {^/T g s c ) 2 = N* , consistent with 
our previous argument that any region of size N* has a 
spectral density of states equal to -k^e E l T <>. Note that 

Eci. (|llfl echoes the free energy profile of droplet growth 
from Eq.©, but unlike Eq.©, it can be used for N > N* 
as well. Ea. l|ll|) explicitly shows that a droplet of size 
larger than N* has an exponentially increasing number 
of available configurations corresponding to lattices typi- 
cal of T g , consistent with the instability of droplets larger 
than N* at temeperatures above T g mentioned earlier. 

The microcanonical argument above is basically a 
gedanken experiment in which we had the demon-like 
ability to browse through all possible atomic arrange- 
ments, given the total number of allowed states equal to 
e Ns c _ rp^g £ f. a j sam pi e i s thus comprised of regions of 
the type considered in the argument (the interface en- 
ergy has been taken into account by the term jy/N, this 
energy may be viewed as the penalty for considering the 
states of a given compact region as if this region were 
totally independent from the rest of the sample; c.f. our 
earlier comments on the locality of the liquid's energy 
landscape). Therefore, if the rate of conversion between 
the alternative glassy states can be ignored, the argu- 
ment immediately yields the density of residual excita- 
tions in a frozen glass: ~ N'a»T e ^ Ta (N*a 3 = £ 3 , of 
course). However, even though each of these imaginary 
regions has an alternative resonant state, there is so far 
only an undetermined chance to reach it within any par- 
ticular time. In fact, the typical classical barrier for the 
excitations available to the regions of size < N* ~ 190 
is F* ~ 39ksT g . Such a barrier would seem to exclude 
the possibility for tunneling for a typical domain of size 
N*. But to account for kinetics issues, we should re- 
peat the argument for the critical size, but also simul- 
taneously include the life-time of each considered con- 
figuration as a selection criterion. In other words, one 
should compute the combined distribution of the excita- 
tion energies, their spatial extent and the corresponding 
tunneling amplitudes. Later in the paper, we will discuss 
one source of correlation between the excitation energy 
and the tunneling amplitude owing to level repulsion ef- 
fects. No netheless here, we present a sim pler argument, 
given in ijLubchenko and Wolvnesl l200lf> . in which the 



tunneling rate distribution is assumed to be independent 
from that of e. Simpler yet, we will look for the density 
of regions that allow for a rearrangement with a zero- 
height barrier. As vindicated post factum, all of these 
simplifications can lead only to at most a 10% error in 
the resulting density of states. 

Imagine the process of conversion to another state as 
a step-wise process where the "nucleus" of this new state 
is increased by adding one atom at a time, as signified by 
the horisontal axis in Fig llUI Such addition involves mov- 
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FIG. 10 The black solid line shows the barrier along the 
most probable path. Thick horizontal lines at low energies 
and the shaded area at energies above the threshold represent 
energy levels available at size N. The red and purple line 
demonstrate generic paths, green line shows the actual (lowest 
barrier) path, which whould be followed if hu>^ < fesT '/2ty. 

ing the atom a distance of the order of the Lindemann 
distance dr,. It follows then that the path connecting 
the two states is likely to encounter a high barrier of the 
order F*, which effectively disconnects those two states. 
However, the possible configurations through which one 
can pass and therefore the barrier heights are distributed 
and there is a chance even for a region of size N* to have 
an arbitrarily low barrier. What would such a distribu- 
tion be? We first have to decide whether the tunneling 
probability is a sum of contributions of many (interfer- 
ing) paths or, whether it is dominated by a single path, 
which has the lowest barrier. The first scenario would 
be realized in a highly quantum glass, where Debye tem- 
peratur e rivals or exceeds the glass t ransition temper- 
atures l)Schmalian and Wolvnesl feOOOlb Such a highly 
quantum glass could in fact melt due to quantum fluc- 
tuations. In our case, freezing is a completely classical 
process, which is signified by the fact that the barriers 
are proportional to a classical energy scale T g . We now 
assume more specifically that the contribution of a tun- 
neling path is proportional to e — ' rV '*/ fiw * , where is a 
quantum frequency scale, a multiple of ujb, and barrier 
V* scales with T g , as mentioned earlier. This would be 
an accurate assignment in the case of a parabolic bar- 
rier. The form of the tunneling amplitude e~^ v con- 
forms to our expectation that the tunneling trajectory is 
dominated by a single path with the lowest barrier, as 
V* and hui$ are taken from distributions characterized 
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by scales fcsT s and const xHlod respectively, the former 
one generally much larger than the latter (the const self- 
consistently will turn out to be less than one in subsection 
lIV.Bf) . Since the energy profile along the tunneling tra- 
jectory has a complicated shape formed by many inter- 
mediate states separated by small intermediate barriers 
(see below), it is fair to say that the state of the system 
at the highest barrier corresponds to the highest energy 
intermediate state (the "transition" state). The statis- 
tics of energy states have already been found earlier. We 
therefore use distribution @ with only one difference: 
we must double the variance because the barrier height 
is actually the difference between two fluctuating quan- 
tities: the energy of the final (or initial) energy and the 
highest energy along the path. As a result, 



fijv(V) ~exp^ s c N ■ 



[V -{T a s c N + ~fVN)} 
A5E 2 N 



(12) 



Distribution (|12fl thus gives the typical value of the bar- 
rier for the (quantum) growth of a droplet. It is easy 
to see from (|12[1 that the highest barrier corresponds to 
rearranging a region of size N$ ~ 14 and is equal to 



1^=^/(2^2-1) 



26T g s c 



(13) 



Since this is the hardest place to get through, we must 
take it as the transition state. Hence, the final distribu- 
tion of (transition state) barriers is the density of pure 
states corresponding to Eq. H12f> with N = N$ (similar to 
the e e / Tg , obtained above). Thus, 



exp 



V t ~ Vrr, 

V2T g 



exp ■ 



-18- s c + 



V2T g 



(14) 

As one can see, the probability to have a small barrier 
path is exponentially suppressed. Nevertheless, owing to 
the large value of the energy parameter in distribution 
(|12fl the fraction of zero barrier paths per mosaic cell 
~ e -i8s c ^ 3 . i[p 7 is actually not prohibitively small. 
A region larger by only 18 molecules (less than a sin- 
gle layer) will have e 18Sc more final states (and therefore 
paths) to go to. We therefore conclude, any region of 
size ~ 200 molecules will have an accessible alternative 
state with spectral density 1/T g . Finally, we stress a re- 
markable feature of the tunneling paths statistics in glass. 
Mark the very rapid - exponential - scaling of the number 
of paths leading out of a particular local structural state 
on the size of the respective region. This means that 
the final estimate of the density of structural transitions 
that have low enough barriers to be thermally relevant 
is rather insensitive to the details of correlation between 
the energy of the transition and its tunneling amplitude. 
Consequently, even a very simple estimate of this density, 
such as the one above, is very robust. Finally, note that 
the tunneling argument above is, again, a microcanoni- 
cal argument, such as the one leading to Eq.©, that also 
takes into account (in a rather crude manner) the mutual 
accessibility between alternative energy states. 



As we will see later, the tunneling barriers, and hence 
the relaxation times of the tunneling centers, are dis- 
tributed. This would lead to a time dependent heat ca- 
pacity. Ignoring this complication for now, the classical, 
long time heat capacity is easy to estimate already (as- 
suming it exists): Since our degrees of freedom span a 
volume £ 3 and their spectral density is 1/T g at low en- 
ergy, one obtains for the low T heat capacity per unit 
volume: T/T g £ 3 , up to an insignificant coefficient. The 
coefficient at the linear heat capacity dependence is often 
denoted P. For silica, T g =1500 K and realistic £ = 20 A 
yields P ~ 6 • 10 45 m _3 J _1 in agreement with the exper- 
iment (we took a — 3. 5 A - a length scale appropriate 
for a tetrahedron formed by four oxygens around a sili- 
con atom. These tetrahedra appe ar to be moving units 
in a-Silica ijTrachenko et all l2000(l '). The assumption of 
the existence of the long-time heat capacity is empirically 
justified (within logarithmic accuracy), but is also con- 
sistent with the present theory, see Subsection IIII. CI and 
Section El 

In conclusion, the main result of this Subsection is 
that the non-equilibrium nature of the glass transition 
results in the existence of residual motional degrees of 
freedom, a significant fraction of which remain thermally 
active down to the lowest temperatures. These degrees of 
freedom are collective highly anharmonic atomic motions 
within compact regions of size £ 3 , determined mainly by 
the length scale of the entropic droplet mosaic deter- 
mined at T g . The energy scale in the spectrum of these 
excitations is set by the glass temperature T g itself. We 
now turn to the question of what determines the strength 
with which these entropic droplet excitations couple to 
the phonons. This will explain the universality in the 
heat conductivity at temperatures below ~1 K. 



B. The Universality of Phonon Scattering 

First of all, what do we mean by "phonons" in amor- 
phous materials? There is no periodicity, therefore one 
can only strictly speak of elastic strain, even if the struc- 
ture is completely stable. In the latter case, low gradient 
strains V</> are still described by a simple bi-linear form: 



Hph = J d 3 r P£| (w) i 



(15) 



The V<^> field is defined on a isotropic, translation- 
ally invariant flat metric, as in continuum mechanics 
ijLandau and Lifshitzl . fl986). and so a wave- vector k is 
an operational concept. It is easily seen, by dimensional 
analysis, that strains arising specifically due to disorder 
will be of higher order in k than the term in Eq. i|15[l : The 
corresponding energy terms should scale with some pos- 
itive power of the lattice inhomogeneity lengthscale(s), 
^nhom (C > 0) , so as to vanish upon "zooming out" . The 
terms will subsequently go as fc 2 (Zj n hom&) < ° ■ But, as we 
have already seen, the amorphous lattice is not stable, 
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that is there are anharmonic transitions with arbitrarily 
small energy and barrier. Still, the regions encompassing 
the transitions are quite small, at most 6 lattice spacings 
across, which is much less than the thermal de Broglie 
wave-length at 1 K (about 10 3 a in a-silica). The unstable 
regions interact with the strains of the otherwise stable 
lattice. We conceptualize this interaction by approxi- 
mating the strain with pure phonons and computing the 
phonon mean-free path. The latter will turn out to be 
about 150 times longer than the phonon wave-length, so 
that the phonon approximation is internally consistent 
in that the phonons are indeed reasonably good quasi- 
particles, at the plateau temperatures and below. Fi- 
nally, note in Eq. Ijl5(l , we have used only one phonon po- 
larization for simplicity (it will be usually obvious how 
to account approximately for all three acoustic phononic 
branches at the end of a calculation; using this "scalar" 
version of the lattice dynamics, for the purposes of this 
paper, boils down to neglecting the difference between 
the longitudinal and transverse sound, except in the later 
discussion of the Griineisen parameter). 

The structural transitions interact with the phonons 
because the energy of the transition changes in the pres- 
ence of a strain: To see this exlicitly, consider the elastic 
energy within a droplet-sized region capable of undergo- 
ing a low energy transition, as relevant at low T . For 
the sake of argument, assume there is no sheer deforma- 
tion; a similar argument applies to the transverse phonon 
branches. The stress energy is then j v d 3 rKu 2 i /2, where 

K is a, compressibility constant on the order of pc 2 , which 
we are allowed to treat as a constant, with the error con- 
tributing in a higher order in fc, as a lready mentioned, uu 
is th e trace of the elastic tensor ijLandau and Li fshitz. 

which has the same meaning as A<f>. Since the 
transition energy is low, the lowest order, quadratic ex- 
pansion suffices. This implies that all individual bonds 
within the region distort by a very small amount (al- 
ready shown not to exceed dj,, even at the transition's 
bottle-neck). We have demonstrated that such regions 
do indeed exist and found their density in the previous 
section. We separate the total elastic tensor ity into a 
contribution <pij due to the elastic stress and dy due to 
the tunneling displacement. The du<t>u cross term repre- 
sents the coupling between the transition and the strain. 
If the phonon wave-length is larger than £, 4>u is constant 
within the integration boundaries and can be taken out 
of the integral. One consequently arrives at the follow- 
ing energy difference for the defect states in the presence 
of a phonon: pc 2 s (j)ii j v d 3 rdu. Here, du corresponds to 
the transition induced displacement between two given 
structural states. Each tunneling center is a multilevel 
system. However, since we are presently interested in 
the coupling of the lowest energy transition to the strain, 
we consider here the set of du corresponding to the two 
lowest energy states of the region. (Higher energy transi- 
tions turn out to be intimately related to the lowest en- 
ergy transition, and are discussed later, in Section IrVT ) 
We therefore conclude that a tunnelling transition, ac- 



tive at low temperatures, is linearly coupled to a lattice 
strain with the strength defined as g — pc 2 J v d 3 rdu/2; 

the corresponding term in the Hamiltonian reads: 

H int =gV^. (16) 

We present next two independent ways to estimate the 
coupling constant g. The first one is based on the realiza- 
tion that at the glass transition, purely phononic excita- 
tions and a frozen-in structural transition must coexist, 
that is they are are of marginal stability with respect 
to each other. On the one hand, a local region posed 
to harbor a structural transition below T g , must not be 
"crumpled" by a passing phonon. On the other hand the 
energy of the transition can only be so high, as to be 
sustainable by the lattice stiffness. In other words, an 
atom will be part a of frozen-in transition, if that atom 
is roughly in equilibrium between the transition driving 
forces and the ambient lattice strain. This stability con- 
dition gives at the molecular scale a, by Eas. i|15fl and 
(|16fl : ger 2 = —pc 2 a?V(f). The lattice strain will be dis- 
tributed throughout the material in the usual manner, 
subject, of course, to the equipartition that fixes the vari- 
ance of the strain so as to conform to the thermally avail- 
abe energy. We take advantage of this by multiplying the 
equilibrium condition by V0 and noting that thermal av- 
eraging is also ensemble averaging. This yields that for 
an atom posed to be part of low energy structural tran- 
sitions below T g , it is generically true that 

| ( g V0 a z ) | ~ pc 2 s ((W) 2 ) a 3 ~ k B T g . (17) 

Noting that (gV</>cr z ) ~ (|gV</>|), one arrives at a simple 
relation: 

g = ^pc 2 s a 3 k B T g , (18) 

which is the main result of this Subsection. It is under- 
stood that this estimate is accurate up to a number of 
order one, and t/'s are likely distributed. At any rate, 
we observe that the TLS-phonon coupling is the geomet- 
ric average between the glass transition temperature and 
an energy parameter pa?c 2 s (~ 10 1 eV) related to the co- 
hesion energy of the lattice (note the quantity pc 2 is a 
multiple of the Young's modulus). We point out the es- 
timate above applies quantitatively to low barrier tran- 
sitions only. The mechanical stability criterion is a zero 
frequency, static condition. However, it takes a finite 
time for a structural transition to respond to an external 
stress. In other words, a region harboring a slow transi- 
tion will likely appear perfectly elastic to a high frequency 
phonon. We thus arrive at the conclusion that TLS- 
phonon coupling must be frequency dependent, however 
deviations from the result obtained above will enter in a 
higher order in u>, k. 

Alternatively, one may attempt to estimate the inte- 
gral over the derivative of the displacement field that 
entered in the expression for the coupling constant g = 
pc 2 J v d 3 vdu/2. Since du is the divergence of a vector, 
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the integral is reduced to that over a surface within the 
droplet's boundary: f s dsd(r), where d(r) is the tun- 
neling displacement itself, near the boundary. How near? 
The Gauss theorem applies so long as the field d is con- 
tinuous. This field is roughly 4 in magnitude close to 
the droplet's center and is zero outside of the region. A 
function defined on a discrete lattice is expressly discon- 
tinuous. Requiring that one be able to cast a continuous 
tunneling displacement field on a discrete manifold, so 
that the field interpolates smoothly between in the 
middle and zero outside, imposes constraints on the d 
values at the droplet's boundary. We argue this value 
should generically go as (a/£)dL, up to a constant, in 
order to realize the interpolation and spread evenly the 
tensile field of the domain wall throughout the droplet. 
In other words, ~ (a/£)dj, is quite obviously the lower 
bound, while higher values statistically imply a higher 
stress, and higher transition energy. Now, since d(r) is 
randomly oriented, the integral over the droplet's bor- 
der is of the order a 2 V N* 2/3 (a/£,)d L: where N* 2/3 is the 
number of molecules at the boundary. The Lindemann 
distance at T g is equal to the magnitude of a thermal 
fluctuation, hence dhja— |V^>| 5 . As a result, the cou- 
pling to the extended defects is still about g ~ \/ pc 2 .a 3 T g , 
again within a factor of two or so. 

Note that, when considering a particular value of 
lattice distortion Vcj) in the discussion above, we did 
not specify the wave-length(s) of the phonons that con- 
tributed to this distortion, therefore the estimate of g in 
Eq. (|18|l is correct as long as the form of the interaction 
term (Eq. jJSt ) is adequate. This surely holds for long- 
wave phonons relevant at the TLS temperatures. 

With the knowledge of g, we can estimate the inverse 
mean free path of a phonon with frequency u. As done 
originally within the TLS model, the quantum dynamics 
of the two lowest energies of each tunneling center are 
described by the Hamiltonian Htls = ecr z/2 + Aa x /2. 
This expression, together with Eas. (|15|l and form a 
complete (approximate) Hamiltonian of the TLS plus the 
lattice vibrations. The phonon inverse m ean fre e path 
is then calcul ated in a standard fashion ij.IackleL Il972t 
PhilliDl ll98l : 



w P H=^tanh(^). (19) 



5 For reference, |V</>| ~ (T g / pc^a 3 ) 1 / 2 at T g is about 0.05 for Si0 2 , 
0.06 for B2O3 (oxide glasses), 0.03 for PS and PMMA (polymer 
glasses), in agreement with the Lindemann criterion. We stress 
the sensitivity to the value of the molecular size a, which is some- 
what arbitrary. Here, we have not calculated the bead size based 
on chemistry, but instead used the values of the speed of sound, 
as employed in the scaling pro cedure of Freeman and Anderson 
(Freeman and Anderson, 1986). According to the definition of 
the Debye frequency, a = (c a /uijy)(6ir 2 ) 1 / 3 . 



This yields 

A rfB /U P = ^-tanh(l/2) (J J , (20) 

where factor 1/3 comes from the averaging with re- 
spect to different orientations of the defects and we used 
P ~ l/T g C 3 . It follows that l mip /\ dB ~ (C/a) 3 ^ 200 
up to a constant of order one. Hence, the analysis 
above explains the universality of the combination of 
parameters Pg 2 /pc 2 , and relates it to the geometrical 
factor (a/£) 3 ~ 10~ 2 , which is the relative concentra- 
tion of cooperative regions in a supercooled liquid, an 
almost universal number within the random firs t or- 
der glass transition theory ijXia and WolvnesLl2000|) . de- 
pending only logarithmically on the speed of quenching. 
Strictly speaking, our argument predicts the universality 
in l m fp/X only within 10%-20% or so. This is a conse- 
quence of indeterminacy of e vs. A correlation that may 
be system specific, or could be due to deviations of the 
£/a ratio from the universal 5.8 at T g . Since the latter 
ratio depends on the glass preparation time, the corre- 
sponding experimental study seems worthwhile. 

Numerically, Eq. l(2TJ|) yields Z mfo/A ~ 70, a factor of 
2 les s than the empirical 150 I Freeman and AndersonL 
1986). We could not have expected much better accu- 
racy from our estimates, that used no adjustable param- 
eters. Although it may seem that we have slightly over- 
estimated the number of scatterers, the size of the error 
is too small to reliably support this suggestion. However, 
this is a good place to make a few comments on the dis- 
tribution of the tunneling matrix elements A, which will 
also prove useful for the discussion of the phonon scat- 
tering at higher frequencies in Chapter II VI The estimate 
for the phonon mean free path in Eci. (|19f) is not terribly 
sensitive to the form of tunneling amplitude distribution 
(Phillips, 1981) (within reasonable limits). This is be- 
cause the contribution of an individual TLS to the total 
scattering cross-section is proportional to A 2 /E 2 , where 
E = VA 2 + e 2 . Two common distributions have been 
used in the literature. One distribution simply assumed 
(A 2 /_E 2 ) ~ 1 in the absence of a more specific knowledge 
and a flat distribution of the total energy splitting E (this 
is actually the original T LS model) . The earlier Standard 
Tunn eling Model (STM) {Anderson et adll972t|phillipsl 
Il972j) . on the other hand, postulates P(A) oc ^ (approx- 
imately supported by our own conclusions too), which 
predicts a nearly flat E distribution as well. In the end, 
both models differ only in that the TLS model has to pos- 
tulate the average (A 2 /E 2 ) value when calculating the 
scattering cross-section (this number is absorbed into the 
TLS-phonon coupling constant g). On the other hand, 
the STM allows for somewhat more closed-form deriva- 
tions, however it still has to introduce the cut-off values 
A m i n and A max as parameters (fortunately, many mea- 
sured quantities depend on these parameters only loga- 
rithmically). (The distinction between the two models 
is described in detail in the front article by Phillips in 
iPhillips].IT98llU One point in favor of the STM is that it 



necessarily predicts time dependence of the specific heat. 
While a time dependence has been observed, its specific 
functional form has not been unambiguously established 
in the experimen t (see (Hunkli nger and Ravchaudhuril 
Il986t IPohi Il98j) '). We also mention, for completeness, 
there is a different way to parametrize the two-level 
system motions within the more general, soft-pote ntial 
model IjBuchenau et all . Il992t iKarpov et a/.Ul983(l . At 
any rate, we see that while a two-level system's contribu- 
tion to the total phonon scattering depends on the value 
of its A ~ E ratio, the precise form of the A distribution 
will change the answer quantitatively, but not qualita- 
tively. We note, that in the context of the present calcu- 
lation it is preferable to consider the simpler, TLS setup 
that does not specify the A distribution, because our ar- 
gument has so far been only semi-classical. Indeed, so far 
the tunneling amplitudes have only interested us from the 
perspective of the volume density of allowed transitions. 
We saw that an indeterminacy of the density could not 
exceed a 10% or so due to a weak (logarithmic) depen- 
dence of that density on a specific A distribution. There- 
fore, we are more confident in the numerical estimates us- 
ing the TLS-model setup that does not require introduc- 
ing additional parameters (such as A m j„) explicitly. Nev- 
ertheless, the special role of A ~ E defects in scattering 
is worth noting. These defects have low classical energy 
splittings e < A and their dynamics are mostly deter- 
mined by the quantum energy scale. These are the "fast" , 
or "zero -barrier" excitations discu s sed earlier in the lit- 
erature ijBlack and Halperinlll977tlGesztill982() . whose 
tunneling matrix element probably can not be directly 
estimated by WKB, but we can still guess that it scales 
with lu£). This suggests that using the same distribution 
function P(A) for all thermally defects may not be jus- 
tified, as circ umstantially supported by results of Black 
and Halperin l)Black and HalperinLll977|) who found that 
the density of TLS derived from the heat capacity and 
conductivity measurements respectively are not exactly 
equal to each other. While this indeterminacy in the ex- 
act barrier distribution intro duces only an error of order 
one in quantitative estimates IjBlack and H alperin. 1977) 
of the density of states and is not of special concern here, 
we note that the present theory, upon inclusion of the ef- 
fects beyond the strict semi-classical picture, does in fact 
provide a mechanism for the excess of the "fast" two-level 
systems, as will be explained in Section Ivl 

Strong Interaction Scenarios. By deriving the density 
of states of structural transitions, and their coupling to 
the phonons based on the known properties of the amor- 
phous lattice, we have constructively established the mi- 
croscopies of glassy excitations in excess to the purely 
elastic excitations. It follows from the discussion that no 
other excitations are present in glasses at 1 K and below 
(see also the discussion on the exhaustive classification 
of excitations in amorphous lattices in Subsection liV.A]l . 
Importantly, the density of states (DOS) in excess of the 
phonons, is due to local motions. This is in contrast with 
Strong Interaction Scenarios (SIS) ijBurin and Kaganl 
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| l996[ICoppersmithl.[T99lllLeggettlll99ltlYu and Leggettl 
1988) that posit that any local excitations (other than 
pure strain) would give rise to a "universal" density of 
states. Such universal density of states arises in SIS as 
a consequence of long range, 1/r 3 , interaction mediated 
by the phonons, so that the actual observed DOS is a 
highly renormalized entity. The corresponding excita- 
tions are expressly non-local, possibly infinite in extent. 
The idea is very attractive because of its generality but 
remains a pure abstraction, until those bare excitations 
are constructively shown to exist in the first place. Ad- 
ditionally, even upon assuming some bare excitations are 
present, demonstrating the quantitative relationship be- 
tween the effective density of states and the phonon cou- 
pling that conforms to the experimental Pq 2 I pc 2 ~ 10~ 2 
has so far proven elusive llCaruzzol Il994l ILeggettl Il999l 
iLubchenko and Wolvnesl l2000|) . On the other hand, we 
have shown that local structural transitions, that inter- 
act with phonons with a particular strength, must indeed 
take place in amorphous solids. In order to determine 
where the current theory stands in relation to the SIS, 
one may inquire whether the phonon-mediated interac- 
tion leads to the emergence of some collective density 
of states. It should be immediately clear that no such 
additional, collective DOS appears at T g , because the ar- 
gument in the previous subsection has already included 
all the effects of the surrounding of a structural tran- 
sition, which simply amounted to the thermal noise at 
T g delivered to the transition by the elastic waves. It, 
of course, does not matter what the phonon source is. 
What happens at low temperatures should be considered 
separately. The effects of interaction turn out to be small 
in the TLS regime and are discussed in detail in the fi- 
nal Section of this paper. Here, we give several qualita- 
tive estimates for the sake of completeness, both at high 
and low temperatures. The phonon-mediated interaction 
goes roughly as ^ 73 , with a numerical factor less than 
one (see Section IVl|) . Ignoring the factor, the interac- 
tion is expressed, with the help of Ea. (|18|) . via the glass 
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transition temperature according to ksT g ^j. Two neigh- 
bouring domains, a distance £ apart, would thus couple 
with strength J ne i g h = T g (a/^) 3 . In order to assess the 
effects of interaction on the effective energies of individ- 
ual transitions, or whether it even makes sense to talk 
about on-site energies after the interaction is turned on, 
one must compare the interaction strength to the width 
of the distribution of the on-site energies as derived in 
the absence of interaction, exactly the same way it is 
done in the context of Anderson localization. According 
to the previous Subsection the latter width, call it AE, 
is of the order T g . The ratio J nc i g h/A-E ~ (a/£) 3 is a 
small number, as expected. There will be no long range 
effects, due to resonant interactions, at high tempera- 
tures near T g . At very low temperatures, only tunneling 
centers (TC) with energy splitting ~ ksT or less are 
thermally active. While the relevant spread of the on- 
site energies Et ~ fcsT is now down by a factor T/T g 
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compared to the glass transition temperature, the con- 
centration of active TC is also down by the same fac- 
tor, namely (T/T g £ 3 ), thus increasing the mutual separa- 
tion between the regions of mobile transitions. The total 
dipole-dipole induced static field due to all those ther- 
mally active two-level systems at a given spot is simply 
^(T/Tg£ 3 ) ~ k B T(a/C) 3 , again much smaller than the 
relevant on-site energy range Et ~ ksT. The motions 
within the tunneling centers are quantum-mechanical at 
these low temperatures, and so one may consider possible 
effects of resonant interactions between distinct TC's, a s 
in the Burin-Kagan scenario ijBurin and Kaeanl Il996fl . 
These effects have been shown to become imp ortant only 
at ul tra-cold temperatures of /iK and below l)Neu et all 
Il997|) . as already mentioned in the introduction. 

Besides having explained the origin of the universality 
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of combination P-^j ubiquitous in the STM, we have also 
seen why the value of this parameter is different from 1, 
suggested by the strong defect-defect interaction univer- 
sality scenario. This value can be traced to the relative 
concentration of the domains (a/£) 3 < 0.01, as just men- 
tioned, ft is curious that the defect-phonon interaction in 
the long wave-length limit can be expressed as a surface 
integral. Besides supporting our picture of the residual 
excitations as motions of domain walls, it points at a con- 
nection with string theories, where the elementary par- 
ticles exhibit internal structure at high enough energies, 
which is also true in our case. In fact, this internal struc- 
ture is ultimately the cause of the phenomena observed 
in glasses at higher temperatures, namely the so called 
bump in the specific heat and the thermal conductivity 
plateau, which are dealt with in Chapter IIVI 

C. Distribution of Barriers and the Time Dependence of 
the Heat Capacity 

The STM postulated tunneling matrix element distri- 
bution P(A) tx 1/A implies a weakly (logarithmically) 
time d ependent heat capacity . This was pointed out early 
on by ijAnderson et all Il972|) , while the first spec ific es- 
timate appeared soon afterwards in (Jackie, 1972). The 
heat capacity did indeed turn out time-dependent, how- 
ever its experimental measures are indirect, and so a de- 
tailed comparison with theor y is difficult. Reviews on 
the s ubject can be found in llNittke et all l!99St iPohll 
119811) . Here, we discuss the A distribution dictated by 
the present theory, in the semi-classical limit, and eval- 
uate the resulting time dependence of the specific heat. 
While this limit is adequate at long times, quantum ef- 
fects are important at short times (this concerns the heat 
condictivity as well). The latter are discussed in Subsec- 
tion EH 

In the tunneling argument from Section flll.AI we have 
suggested a WKB type expression for the tunneling am- 
plitude: 



which would be completely correct in the case of a 
parabolic barrier with frequency ufi and height V$ and 
was motivated by the necessity to maintain the proper 
scaling with H in the denominator of the exponent, given 
that the typical barrier height is determined by the clas- 
sical landscape characteristics and should scale with T g . 

According to Ea. lfblj) . P(V^) cx e^ T i . It follows then 
that 

P(A)dA = A(^J^, (22) 

where c = tuv$ /V%T g should be less then 0.1 according to 
our estimates of or (see section HV.B|I . A is a constant, 
to be commented on very shortly. The distribution in 
Ea. (|2*2*|) becomes Pstm(A)gJA cx dA/A postulated in the 
standard tunneling model, if c — > 0. As shown next, the 
non-zero c gives rise to an anomalous exponent a = c/2 
in the heat capacity C cx T 1+a and a power law t c ^ 2 for 
the specific heat time dependence at long times, as op- 
posed to a logarithmic one, predicted by the STM. While 
c ~ 0.1 implies a ~ 0.05, experimentally, a seems to 
vary from 0.1 to 0.5. This larger value is consistent with 
quantum mixing effects that go beyond the scmiclassical 
analysis as we will discuss later. 

Scaling Tvur in the denominator of the tunneling ex- 
ponent implies that must be a quantum energy scale 
and it is indeed shown in Section fl V .Bl that ui* is propor- 
tional to the Debye frequency u>d- While the tunneling 
argument from Section IIII.AI only explicitly considered 
the statistics of the highest energy state along the tun- 
neling trajectory, the expression in Ea. (|2T]> actually does 
not use such a simplified picture but considers a finite 
vicinity of the barrier top. The conclusion of Section 
IIV.BI that the barrier heights are distributed exponen- 
tially, such as in in Ea. H14J) . remains true in either case. 
The leads to a non-zero value of c, and here we explore 
what consequences this has for the low temperature heat 
capacity and conductivity. As follows from the discus- 
sions in Section Till. Al constant A in Ea. l)22|l is of order 
one. 

Since the temperatures in question here are so low (IK 
and below), we will ignore the energy dependence of n(e) 
in this section and take n(e) = P. In order to see the 
time dependence of the heat capacity we obtain the com- 
bined distribution of the TLS energy splittings E and 
relaxation rat es r _1 - P(E,t^ 1 ) 1 much as was done in 
lUackl4lT97l . - and then count in only those TLS whose 
relaxation time r is shorter than a particular experimen- 
tal observation time t. 

The ( phonon irradiation induced) relaxation rate of a 
TLS is ll,IackleUl972t) : 

T -i^Zp^E CQtHf3E/2) (23) 



A = A e wf, 



(21) 



It follows from Eqs. and that 
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P{E,t~ 1 ) = PA 
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Ao 
E 



2vT 



(24) 



r 



where 



~iin( E ) 



^L C oth(f3E/2) (25) 



is the fastest relaxation rate of a TLS with energy split- 
ting E, achieved at A = E, of course. As follows from 
Ea. (|25ll . the rate scales roughly as the cube of temper- 
ature and is empirically of order an inverse millisecond 



at 10 2 K. The resultant sample's heat capacity per unit 
volume is then: 



dr- 1 PiE,^ 1 ), 
(26) 

where \f)E/2 cosh(/3£:/2)] 2 is the TLS heat capacity. 
With a change of variables, Eq. lffi)|l reads: 
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(27) 
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At long times the expression Ea. H27l) yields a power 
law for both time and temperature dependence of the 
specific heat: 



lim C(t) cx t°/ 2 T 1+c / 2 , 



(28) 



where, note, the temperature dependence also comes 
from the energy dependence of T m i n (E) oc E~ 3 in 
Ea. (l27|l . The value c ~ 0.1 implies the long time 
heat capacity should obey C oc T 1+a at low T, where 
a ~ 0.05, a smaller number than observed in amor- 
phous materials. We must bear in mind that the issue 
of the exact form of the time dependence in the labora- 
tory still appears to be unresolved, as there is no definite 
agre ement between different experiments; for references , 
see Hunkling er and Ravchaud huril 1 19861 iNittke et all 
119981 iPohli 119811 ISahling et all 1200 jT while there is 
no doubt that the specific heat is time-dependent, some 
experiments agree with the logarithmic time profile, as 
predicted by the STM, others give a lower or higher speed 
of variation with time. The present semiclassical predic- 
tion with c = 0.1 would be hard to distinguish from a 
logarithmic law. Finally, even though a correction to 
the linear temperature heat capacity dependence with 
c = 0.1 is most likely smaller than what is experimentally 
observed, the value of this correction is non-universal, 
consistent with empirical data. 

The expression in Ea. (|27jl can be evaluated numeri- 
cally for all values of t, and the results for three different 
waiting times are shown in Fig llll for c = 0.1. The value 
of T m i n = 2.0/isecat E/Td = 5.7- 10~ 4 , derived from the 
prese nt theory (also consistent with ijGoubau and Taut 
was used. The results for t = lOfisec demonstrate 
that due to a lack of fast relaxing systems at low ener- 



gies, short time specific heat measurements can exhibit 
an apparent gap in the TLS spectrum. Otherwise, it 
is evident that the power-law asymptotics from Ea. H28(l 
describes well Ea. (|2*7jl at the temperatures of a typical 
experiment. 




T/T n 



FIG. 11 Displayed are the TLS heat capacities as computed 
from Ea. H27^ appropriate to the experiment time scales of or- 
der a few microseconds, seconds and hours. Value of c = 0.1 
was used here. If one makes an assumption on the specific 
value of Ao, it is possible to superimpose the Debye contribu- 
tion on this graph, which would serve as the lowest bound on 
the total heat capacity. As checked for Ao = u)d, the phonon 
contribution is negligible at these temperatures. 

As clear from the discussion above, the long-time 
power law behavior of the heat capacity is determined by 
the "slow" two-level systems corresponding to the higher 
barrier end of the tunneling amplitude distribution, ar- 
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gued to be of the form shown in Eq. i|22[l . If one assumes 
that this distribution is valid for the zero-barrier tail of 
the log(A) distribution as well, one would expect that 
the heat conductivity would scale as T 2+c at the TLS 
temperatures, in contrast to an observed experimental 
sub-quadratic dependence T 2 ~ a . As we shall see in Sec- 
tion other quantum effects are indeed present in the 
theory and we will discuss then how these contribute 
both to the deviation of the conductivity from the T 2 
law and the way the heat capacity differs from the strict 
linear dependence, both contributions being in the di- 
rection observed in experiment. Finally, when there is 
significant time dependence of cy, the kinematics of the 
thermal conductivity experiments are more complex and 
in need of attention. When the time-dependent effects 
are included, both phonons and two-level systems should 
ideally be treated by coupled kinetic equations. Such 
kinetic analysis, in the context of the time dependent 
hea t capacity, has been conducted before by other work- 
ers l|Strehlow and Meissneii Il99^) . 

IV. THE PLATEAU IN THERMAL CONDUCTIVITY 
AND THE BOSON PEAK 

In this section we continue to explore the consequences 
of the existence of the low temperature excitations in 
amorphous substances, which, as argued in Chapter IIHI 
are really resonances that arise from residual molecular 
motions otherwise representative of the molecular rear- 
rangements in the material at the temperature of vitrifi- 
cation. We were able to see why these degrees of freedom 
should exist in glasses and explain their number density 
and the nearly flat energy spectrum, as well as the uni- 
versal nature of phonon scattering off these excitations 
at low T (< 1 K). 

At higher temperatures {K B T ~10~ 2 ' to lO -0 ' 5 ?^), 
an apparently different kind of excitations begins to ap- 
pear, leading to the so called bump in the heat capac- 
ity and plateau in the thermal conductivity, as was dis- 
cussed in the Introduction. We argue in this chapter that 
the transitions between the mutually accessible frozen-in 
minima in the amorphous lattice that give rise to the two 
level system behavior at the lowest T also explain the ex- 
istence of the modes responsible for this "Boson Peak" 
and the intense phonon scattering at the corresponding 
frequencies. This thus removes the need to invoke theo- 
retically any additional mechanisms, although other con- 
tributions may well be present to some extent; we will 
try to assess this possibility in the following Section. 

A. Introduction: Classification of Excitations in Glasses 

While we believe to have mostly achieved a microscopic 
understanding of the excitations that are specific to the 
amorphous solids and are not present in other types of 
materials, this description is rather new and, naturally, 
there is certain lack of established language that could be 



efficiently used to characterize these excitations. In this 
section, we will introduce some terminology that will be 
used in the rest of the article. At the same time, we will 
provide a brief general analysis of what possible quali- 
tatively distinct types of molecular motions can exist in 
glasses. 

Any atomic motions that take place in a frozen glass, 
obviously may also be present in the liquid above T g . 
For example, those high T motions that correspond 
to shear attain stiffness (on realistic time scales) be- 
low freezing. The motions in the liquid, apart from 
pure volume change, corresponding to the longitudinal 
sound, are molecular translations, or, informally speak- 
ing, jumps. Above Ta, such jumps are not accompanied 
by a noticeable volume change and bond stretching, as no 
metastable structures form in the liquid at these temper- 
atures. The barriers are therefore largely entropic. (It is 
nice to compare Feynman's disc ussion of the abs ence of 
energy barriers in superfluid He ijFevnmanl . Il954j) in this 
regard.) Below Ta, such hopping already involves mov- 
ing a number of molecules from one local minimum of 
the free energy functional to another such minimum and 
thus requires structural rearrangement within a certain 
cooperative length owing to the formation of metastable 
local arrangements. Molecular translations do not con- 
serve momentum, which subsequently must be provided 
by the rest of the bulk. We thus call these degrees of 
freedom, which are relics of the translational motions in 
the liquid above T g , inelastic degrees of freedom. They 
are truly inelastic also in the macroscopic sense of the 
word, because the existence of alternative configurations 
in the solid bulk, which are also coupled to the phonons, 
ultimately leads to irreversible relaxation, if the sample 
is subjected to mechanical stress thus causing a shift in 
the thermal population of the alternative internal, struc- 
tural states. T his is the mechanism behind the so called 
bulk viscosity l)Landau and Lifshitzl . 11986^ (incidentally, 
it also contributes to the so called relaxational phonon 
absorption, which we discuss in subsection lIV.Gfl . The 
switching from one energy minimum to another is ac- 
complished by moving the domain wall - the interface 
between the two alternative configurations - through the 
local region. As mentioned earlier, this domain wall is 
something of an abstract entity, really a quasiparticlc of 
a sort. Yet it has many ponderable attributes. For one 
thing, it has a mass (per unit area), which will be ob- 
tained in this section. It also has surface tension, there- 
fore it can support surface vibrations, again, of a sort. Al- 
though these vibrations are realized through real atomic 
motions, it is more beneficial to think of them as vibra- 
tional modes of an imaginary membrane. In fact, as will 
be argued later, the oscillations of this membrane cor- 
respond to the indeterminacy in the exact boundary of 
the frozen-in domain that has more than one kinetically 
accessible internal state. Therefore, highly anharmonic 
atomic motions in the real space correspond to harmonic 
motions in the space where the domain walls are defined. 
This mental construction does the trick of enabling us to 
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calculate the ripplon spectrum, as de monstrated in sec- 
tion llV.Cl Now, since it was shown in ijXia and Wolvnesl 
2000) , that the liquid degrees of freedom below Ta con- 
sist of switching to alternative local energy minima; we 
can claim our assignment of different inelastic modes is 
exhaustive (but not unique, of course!). These are, again, 
translations and vibrations of the domain walls. 

On the other hand, any purely elastic motions in the 
glassy lattice can be thought of as a sum of ordinary, 
affine, displacements a nd no n-affine displacements (see 
e.g. (|Wittmer et all l200l[l '). The affine component 
would be the only one present in a perfectly isotropic 
medium and would follow the stress pattern according 
to a Poisson equation (the situation with a non-isotropic 
crystal is conceptually the same). The non-affine dis- 
placements are a consequence of the absence of period- 
icity. They involve a small number of molecules and are 
characterized by a non-zero circulation of the displace- 
ment field, ft is not clear at present whether the size 
of these non-affine "islands" could be inferred in present 
day computer simulations, since the amorphous struc- 
tures that can be currently generated on modern com- 
puters still correspond to unrealistically rapid quenching 
rates. The resulting structure corresponds to a sample 
caught in a very high energy state with extraordinar- 
ily low barriers. As is clear by now from the random 
first order transition theory, such structures correspond 
to temperatures close to Ta and will have very small co- 
operative regions approaching the molecular scale a. 

We conclude this subsection by repeating ourselves 
that one important difference between the elastic and the 
inelastic modes is in how they absorb the phonons. While 
any static disorder can only provide Rayleigh scattering 
with a characteristic length scale equal to the size of the 
heterogeneity, the inelastic (resonant) absorption's cross- 
section scales as the square of the phonon wave-length, it 
thus will considerably dominate the Rayleigh mechanism 
for the longer wave-length phonons (absorption satura- 
tion in the TLS's does not occur at the sound intensities 
typical of heat transfer). 



B. The Multilevel Character of the Entropic Droplet 
Excitations 

We hope to have convinced the reader by now that 
the tunneling centers in glasses are complicated ob- 
jects that would have to be described using an enor- 
mously big Hilbert space, currently beyond our compu- 
tational capacity. This multilevel character can be antic- 
ipated coming from the low tem perature perspective in 
l)Lubchenko and Wolvnesl l200ll) . Indeed, if a defect has 
at least two alternative states between which it can tun- 
nel, this system is at least as complex as a double well po- 
tential - clearly a multilevel system, reducing to a TLS at 
the lowest temperatures. Deviations from a simple two- 
level behavior have been seen dire ctly in single-molecule 
experiments ijBoiron et all Il999|) . In order to predict 



the energies at which this multilevel behavior would be 
exhibited we first estimate the domain wall mass. Ob- 
viously, the total mass of all the atoms in the droplet is 
so large that the possibility of simultaneous tunneling of 
all atoms is completely excluded. The tunneling, we ar- 
gue, occurs stagewise; each individual motion encounters 
a nearly flat potential, implying low frequency instanta- 
neous modes. 

In addition, the effective mass of the domain wall turns 
out to be low, also owing to the collective, barrierless 
character of the tunneling events. This is because moving 
a domain wall over a molecular distance a involves dis- 
placing, at any one (imaginary!) time, individual atoms 
only a Lindemann length d^. Suppose this occurs on 
the (imaginary) time scale r. The resulting kinetic en- 
ergy is M w (a/r) 2 = N w m(d L /T) 2 , where N w ~ (£/<z) 2 
is the number of molecules in the wall and m is the 
molecular mass. Thus the mass of the wall M w is only 
m(£/a) 2 (d L /a) 2 . Using (£/a) ~ 5.8 and (d L /a) 2 ~ 0.01 
gives M w ~ tto/3. This implies the mass of the wall per 
atom is very small - about a hundredth of a molecular 
mass, consistent with the simulatio ns of certain barrier - 
less dislocation motions in copper ijVegge et all l200l[) . 
Using (dL/a) 2 ~ ksTg/pc 2 ., derived earlier, one can ex- 
press the wall's mass through the material constants as 
M w ~ (^/a^ksTg/cj.. The wall mass estimate above, in- 
spired by the Fey nman's argumen t on the effective mass 
in liquid helium ijFevnmanl Il953|) . is entirely analogous 
to the well known estimate of the sol iton mass in poly- 
acetylene, see e.g. (|Heeger et all . 1988). In the latter, the 
soliton moves a large distance, while individual atoms un- 
dergo only small displacements leading to a low soliton 
mass. 

We can now use the typical value of the barrier cur- 
vature from our tunneling argument in section lill.AI fsee 
Fig llOf) to estimate the typical frequency of motion 
at the tunneling barrier top. We now express the bar- 
rier profile V(N) as a function of the droplet 's radius 
r = a(3N/4:iT) 1 / 3 and obtain 

ut = -d 2 V/dr 2 /M w ~ I.6(o/£Vd- (29) 

According to the quantum transition state theory 
iWolvnes). fl98l . and ignoring damping, at a temper- 
ature T ~ fu J j t /2Trk B ^ [a/£)T D /2ir, the wall mo- 
tion will typically be classically activated. This tem- 
perature lies within the plateau in thermal conductivity 
ijFreeman and Andersonl Il986j) . This estimate will be 
lowered if damping, which becomes considerable also at 
these temperatures, is included in the treatment. Indeed, 
as shown later in this section, interaction with phonons 
results in the usual phenomena of frequency shift and 
level broadening in an internal resonance. Also, activated 
motion necessarily implies that the system is multi-level. 
While a complete characterization of all the states does 
not seem realistic at present, we can extract at least the 
spectrum of their important subset, namely those that 
correspond to the vibrational excitations of the mosaic, 
whose spectral/spatial density will turn out to be suf- 
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ficiently high to account for the existence of the Boson 
Peak. 



C. The Vibrational Spectrum of the Domain Wall Mosaic 
and the Boson Peak 

At low temperatures the two-level system excitations 
involve tunneling of the mosaic cells typically containing 
N* ~ 200 atoms. The tunneling path involves stage- 
wise motion of the wall separating the distinct alterna- 
tive configurations through the cell untill a near reso- 
nant state is found. At higher temperatures, other final 
states are possible since the exact number and identity 
of the atoms that tunnel can vary. These new config- 
urations typically will be like the near resonant level 
but will also move a few atoms at the boundary, i.e. 
at the interface to another domain. This is schemati- 
cally shown in Fig ll2l Alternatively, due to the quan- 




FIG. 12 Tunneling to the alternative state at energy e can 
be accompanied by a distortion of the domain boundary and 
thus populating the ripplon states. The doubled circles denote 
atomic tunneling displacements. The blue line signifies, say, 
the lowest energy state of the wall, and the blue circles corre- 
spond to the respective atomic displacements. An alternative 
wall's state is shown in red, the corresponding alternative sets 
of atomic motions are coded using different colors. The do- 
main boundary distortion is shown in an exagerated fashion. 
The boundary does not have to lie in between atoms and is 
drawn this way for the sake of argument; its position in fact is 
not tied to the atomic locations in an d priori obvious fashion. 

turn mechanical uncertainty of the exact location of the 
domain wall, its shape is intrinsically subject to fluc- 
tuations (these are zero-point vibrations of the domain 
wall). It is thus not surprising that the ripplon's fre- 
quencies turn out to be proportional to u>d, the basic 
quantum energy scale in the system. These fluctuations 
of the domain boundary shape can be visualized as do- 
main wall surface modes ("ripplons"). A detailed cal- 
culation of the ripplon spectrum would require a con- 
siderable knowledge of the mosaic's geometry. At each 
temperature below Ta the domain wall foam is an equi- 
librium structure made up of flat patches of no ten- 



sion (remember the renormalized <j(r) cx r -1 / 2 ; how- 
ever fluctuations will give rise to finite curvature and 
tension). To approximate the spectrum we notice that 
the ripples of wave-length larger than the size of a patch 
will typically sense a roughly spherical surface of radius 
R = ^(3/47r) 1 / 3 . The surface tension of the mosaic has 
been calculated from the classical microscopic theory and 
is given by a(R) = MknT ja 2 ) log((a/d L ) 2 / 7 re)(a/i?) 1 /2 
ijXia and Wolv nes. 200(J, where (Il/o, is the universal 
Lindcmann ratio. It could appear that such tension 
could collapse an individual fragment of the mosaic but 
this tension is, of course, compensated by stretching the 
frozen-in outside walls. We approximate the effect of this 
compensation by an isotropic positive pressure of a ghost 
(i.e. vanishing density) gas on the inside. 

The eigen-frequency spectrum of the surface modes of 
a hollow sphere with gas in side is well known (see e.g. 
ijMorse and Feshbachlfl9 53?l. as well as our appendix |A"|> . 
If we pretend for a moment that the surface tension coef- 
ficient a is curvature independent, the possible values of 
the eigen-frequency lu are found by solving the following 
equation: 



col[n/U'./?,/V,,)] = ( ) - 
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(30) 

where pw is the membrane's mass per unit area, p g 
and c g are gas' mass density and sound speed respec- 
tively. As stated earlier, Ea. H30|l only applies for I > 2. 

Finally, function cot[a;(z)] = (-1 + z lj ±f^ 

ji (z) is the spherical Bessel function of l-th. order, does ex- 
hibit behavior similar to that of the regular trigonometric 
cotangent for arguments of the order unity and larger, 
going however to — l/l as z — > 0. Its graph for I = 2 
is shown in Fig ll3l An inspection shows that for each 



where 




FIG. 13 The functions entering Ea. pijl 

are shown for some arbitrary parameter values. Here, / = 2. 

I, the smallest solution of Ea. H30|) gives the frequency of 
the proper eigen-mode of the shell itself (shifted due to 
the presence of the gas inside) , whereas the rest of the so- 
lutions represent the standing acoustic waves in the gas. 
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This is especially clear in the p g — > limit, when the 
lowest frequency does not even depend on the gas' sound 
speed, whereas the rest of the solutions are obviously de- 
termined by the inverse time it takes the sound in the 
gas to traverse the sphere. 

Since we are interested only in the wall's proper modes 
in the limit p g — > 0, we get unambiguously for the fre- 
quency of an Z-th harmonic: 

^ = (l-l)(l + 2)(-^\- (Z>2). (31) 

Accounting for the unusual r dependence of the sur- 
face tension o~(r) oc r^ 1 / 2 modifies the standard result 
from Ea. (|31fl by a factor of 9/8. The reason is, the pe- 
culiar surface energy dependence F sur f(R) = A-kR 2 <j = 
4-7rcro-R 3 ' /2 a 1 / 2 calls for the following dependence of pres- 
sure on the curvature: p = gja^Cg) _ 3 cr ^ ag com _ 
pared to the regular p — 2-^). The eigen- frequencies, in 
their turn, are determined by calculating the (frequency 
dependent) excess pressure due to a variation in curva- 
ture. Since now p oc i?~ 3 / 2 , varying p with respect to 
R brings down another factor of 3/2, thus giving 9/4 in- 
stead of the 2 of the curvature independent case. Hence 
the (barely s ignificant, but curious) correctio n factor of 
9/8 used in flLubchenko and Wolvnel l2003a|) . bmce we 
have been assuming that the amplitude is infinitesimally 
small, this factor is the only consequence of having a cur- 
vature dependent a, which should have made the mem- 
brane oscillations even more non-linear (as compared to 
a = const case) in the case of finite displacements. Pin- 
pointing this effect, however, is clearly beyond the accu- 
racy attempted by the present model. Finally, one finds 
a spectrum with 

^ = lp~^ {l ~ 1){l + 2)] (Z - 2) ' (32) 

where each l-th mode of a sphere is (21 + l)-fold degen- 
erate. Using pw = (rfi/o) 2 pa, obtained earlier in the 
chapter and T g ~ pc 2 a 3 (d,L/a) 2 fsection flII.B|) . one finds 

-1.34 UB ( fl /0 5/ V((-l)(' + 2)/4 

~ 0.15 l)(/ + 2)/4. (33) 

Because of the unive r sality of the (a/£) ratio 
l|Lubchenko and Wolvnel [2001), u>i is a multiple of the 
Debye frequency. Apart from the barely significant 
(a/^) 1 / 4 factor, again, due to the R dependent a, the 
ubiquitous scaling u>i ~ (a/£)cj£> stresses yet another 
time the significance of the scale £. Such a scale has been 
previously empirically deduced by interpreting inelastic 
scattering experiments but has been usually ascribed to 
the static heterogeneity length scale, in contrast with the 
dynamical nature of the mosaic in the present theory. We 
note, again, that this "static heterogeneity" has never 
been unambiguously seen in X-ray diffraction. Owing to 
the material's discreteness, one does not expect harmon- 
ics of higher than tt (^N*) 173 [(R - a/2)/R] ~ 9. .10th 



order, a relatively large number, which justifies the tac- 
itly assumed continuum approximation. The lowest al- 
lowed ripplon mode is I = 2 (corresponding frequency 
is ~ ITHz for silica, in remarka ble agreement with 
the in elastic neutron scatering data ijWischnewski et all . 

mn)). 

The requirement I > 2 can be understood from the 
symmetry considerations. 1 = 1, the case of no restoring 
force, corresponds to a domain translation. Within our 
picture, this mode corresponds to the tunneling transi- 
tion itself. The "translation" of the defects center of 
mass violates momentum conservation and must be thus 
accompanied by absorbing a phonon. Such resonant pro- 
cesses couple linearly to the lattice strain and contribute 
the most to the phonon absorption at the low temper- 
atures, dominated by one-phonon processes. I = 0, on 
the other hand, corresponds to a uniform dilation of 
the shell. This mode is formally related to the domain 
growth at T > T„, and is described by the theory in 
l)Xia and Wolvnesl l2000|) . It is thus possible, in princi- 
ple, to interpret our formalism as a multi-pole expansion 
of the interaction of the domain with the rest of the sam- 
ple. Harmonics with I > 2 correspond to pure shape 
modulations of the membrane. 

The existence of the domain wall vibrations explains 
and allows us to visualize, at least in part, the multi- 
level character of the tunneling centers as exhibited at 
temperatures above the TLS regime. Curiously, the ex- 
istence of TLS's, even though displayed at the lowest T, 
is basically of classical origin due to the non-equilibrium 
nature of the glassy state. Yet the ripplons, even though 
seen at higher T, are mostly due to quantum effects and 
would not be predicted by a strictly semi-classical the- 
ory, in which h — > 0. A schematic of the resultant droplet 
energy levels is shown in Fig ll4l The arrangement of the 




FIG. 14 Tunneling to the alternative state at energy e can 
be accompanied by a distortion of the domain boundary and 
thus populatiting the ripplon states. All transitions exem- 
plified by solid lines involve tunneling between the intrinsic 
states and are coupled linearly to the lattice distortion and 
contribute the strongest to the phonon scattering. The "ver- 
tical" transitions, denoted by the dashed line, are coupled to 
the higher order strain (see Appendix[XJ and contribute only 
to Rayleigh type scattering, which is much lower in strength 
than that due to the resonant transitions. 
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combined internal (configurational) and ripplonic density 
of states, as depicted in Fig^] has the following moti- 
vation behind it. We include the possibility of distorting 
the domain wall during the tunneling transition by pro- 
viding a set of vibrational states on top of the alterna- 
tive internal state. The arrangement of the energy states 
as depicted in Fig ll4l insures that only thermally active 
tunneling centers have mobile (and thus vibrationally ex- 
citable boundaries). The atomic motions at the inactive 
defects' sites (i.e. that cannot tunnel or cross the bar- 
rier) would be indistinguishable from the regular elastic 
lattice vibrations. Importantly, direct transitions to the 
ripplonic state can occur from one of the two lowest - 
"TLS" - energy states of a tunneling center. This in- 
herent assymetry between the two structural states of 
a tunneling center actually reflects the thermodynami- 
cal inequivalence of the two states at the glass transition 
temperature. While one of the states represents the lo- 
cal structure in (meta-stable) equilibrium with the cur- 
rent liquid arrangement around it, the other state is a 
configuration that must only be regarded as one of the 
structures along the many escape routes from the cur- 
rent equilibrium local state. At T gi most of those escape 
routes become too costly energetically. 

This is a good place to remind the reader that ex- 
isting explanations of the large density of states at the 
BP energies have to do either with purely harmonic ex- 
citations of disordered, but perfectly stable lattice (see 
Introduction), or by generalizing the low energy inelas- 
tic, two-level degrees of freedom to multilevel systems, 
as was done e.g . by the soft potential m odel (SPM) 
l|Buchenau et all Il992t iKarpov et all . Il983|) . Such gen- 
eralizations imply a connection between the anomalies 
seen in the TLS regime and at these higher energies. Such 
a connection is strongly suggested by experiment, most 
prominently by the strength of phonon scattering. The 
latter is inelastic at the BP energies, as it was at the 
TLS energies. We stress, the rate of increase of the rip- 
plonic density of states is much much higher than that 
empirically assumed in the purely empirical SPM. Again, 
there is virtually no freedom to adjust the numbers in our 
theory. 

In order to compute the heat capacity of the ripplons 
on top of the structural transitions we will need to con- 
sider the (classical) density of the inelastic states in more 
detail than in the previous section. The density of states 
n(e) = ijre t l Ta was derived earlier taking as the reference 
state the generic global liquid state corresponding to the 
(high-energy) configuration frozen- in at T g . It turns out 



that only transitions to states with e < (relative to 
the liquid state!) will contribute to the TLS density of 
states. Indeed, as we have shown, the size of the re- 
gion that permits a low-barrier rearrangement must be 
slightly (by ~18 molecules) larger than the generic co- 
operative size at T g . On the other hand, we know from 
the RFOT theory that larger cooperative regions corre- 
spond to lower energy liquid structures. Therefore one of 
the two alternative states must be lower in energy than 
the generic liquid state at T g . As a result, the negative 
e's correspond to some of the very numerous but mostly 
unavailable lower lying energy states, now accessed by 
tunneling. Now, if each of those true local ground states 
is taken as the reference one, the spectral density will be 
now n(e) = -Le~ c / Ts (e > 0). We consequently can let 
e from Fig^Jtake both positive and negative values by 
writing 

„( e ) = i- e -M/r.. (34) 
1 g 

We can now calculate each domain's partition function 
by including all possible ways to excite the system: 

Z e = 1 + ]T e-^+^lm ni m u, lm ) = 1 + e -0e Z 2l+1 } 
{n, m } I 

(35) 

where Zi = 1/(1 — e /3uJl ) is the partition function of an 
Ith order ripplon mode and we used m = —I.. I. Here we 
assume each ripplon is a harmonic oscillator. Note that 
since the "harmonic" excitations of frequency lji are on 
top of another (structural) excitation, we must consider 
the issue of the zero-point energy of these "harmonic" 
excitations, that is no longer a matter of simply choosing 
a convenient reference energy. Note that this zero-point 
energy is actually several orders of magnitude higher than 
the subKelvin energies that are sufficient to excite some 
of the local structural transitions. And indeed, the en- 
ergy that comprises the ripplons' ground state energy 
is not extracted from the thermal fluctuations of the 
medium, but, one may say, is simply "converted" from 
the zero-point energy of local elastic vibrations of the 
lattice. At the site of a "slow" (or, thermally inactive) 
structural transition, domain wall vibrations are indis- 
tinguishable from the regular lattice phonons, as already 
mentioned. 

The specific heat corresponding to the partition func- 
tion in Eq. i|35|) is found by computing c e — 2 d i° 9 g' : 
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Expression i|36|) clearly becomes the TLS specific heat 

°TLS = ( 2cos h( £ ;3 e /2) ) forT< W; . 

In order to obtain the amorphous heat capacity per 
domain, we (numerically) average c e with respect to n(e); 
the result is shown in figure IT51 with the thin solid line. 




FIG. 15 The bump in the amorphous heat capacity, di- 
vided by T 3 , as follows from the derived TLS + ripplon 
density of states. The thin curve corresponds to Eg. 1361 . 
The thick s olid line is experimental data for a-Si02 from 
llPohl ll98lT) . The experimental curve, originally given in 
J/gK 4 , was brought to our scale by being multiplied by 
h 3 pc 3 s {6n 2 )(^/a) 3 /k%, where we used lo d = (c s /a){6-K 2 ) 1/3 , 
e/a) 3 = 200, p = 2.2g/cm 3 ,Cs = 4100m/sec and T D = 342K 
Freeman and Anderson! H986). The other curves take into 
account effects of friction and frequency shift in the ripplon 
frequencies. They will be explained later in subsection IIV.FI 
The Debye contribution was included in our estimate of the 
total specific heat; it was calculated according to (per par- 
ticle) cd = 9(r/Tp) 3 f T/T ° dx ^f^ llKittelll956ft . This 

equals to 234(T/T ) 3 at the low T. When multiplied by 
(£/a) 3 ~ 200, this gives a value 4.6 ■ 10 4 in good ag reement 
with the experiment (see Fig. 3. 10 from iPohll Il9 81)); note, 
however that the amorphous Tb is lower than the correspond- 
ing crystalline one, still it seems T D morph > \T c D ryst . We 
remind the reader that no adjustable parameters have been 
used so far. 



D. The Density of Scatterers and the Plateau 

In order to estimate the phonon scattering strength 
and thus the heat conductivity, we need to know the effec- 
tive scattering density of states, the transition amplitudes 
and the coupling of these transitions to the phonons. 

Any transition in the domain accompanied by a change 
in its internal state is coupled to the gradient of the elas- 
tic field with energy g ~ pc 2 s J dsd(r), where d(r) is the 
molecular displacement at the droplet edge due to the 
transition (see section UlLBfl . An additional modulation 
in the domain wall shape due to the current vibrational 
state cancels out due to the high symmetry (I > 2), as 
easily seen when computing the angular part of the sur- 
face integral. We therefore conclude that any transitions 



between groups marked with solid lines in Fig ll4l are cou- 
pled to the phonons with the same strength as the under- 
lying (TLS-like) transition. (Notice this also implies in- 
elastic scattering off those transitions!) Incidentally, no 
selection rules apply for the change in the ripplon quan- 
tum numbers, being essentially a consequence of strong 
anharmonicity of the total transition. 

We do not possess detailed information on the transi- 
tion amplitudes, however they should be on the order of 
the transition frequencies themselves, just as is the case 
for those two level systems that are primarily responsi- 
ble for the phonon absorption at the lower T which also 
have their transition amplitudes comparable to the total 
energy splitting. The argument is thus essentially the 
same as proposed earlier in section 1111. Al It should be 
noted, however, that the Hilbert spaces corresponding to 
the quantum in nature ripplons and the classical inelastic 
states are quite distinct (although overlapping); it thus 
should not be surprising that the matrix element beween 
superpositions of these spaces is on the order of the energy 
differences themselves. In what follows, we circumvent to 
an extent the question of what the precise distribution of 
the tunneling amplitudes of the TLS+ripplon transitions 
is and simply calculate the enhancement of the bare TLS 
induced scattering due to the presence of the ripplons. 
This is suggested by an earlier notion that the struc- 
tural transitions in glasses couple to the phonons with 
the same strength even if accompanied by exciting vibra- 
tional modes of the mosaic. 

We now calculate the density of the phonon scattering 
states. Since we have effectively isolated the transition 
amplitide issue, the fact of equally strong coupling of all 
transitions to the lattice means that the scattering den- 
sity should directly follow from the partition function of 
a domain via the inverse Laplace transform. We will not 
proceed this way for purely technical reasons. In addi- 
tion, we will separate the cases of positive and negative 
e (see Fig ll4[l . corresponding to absorption from ground 
and excited states respectively. 

The phonon-ripplon interaction exhibits itself most ex- 
plicitly through the phonon scattering, which becomes 
so strong by the end of the plateau as to cause complete 
phonon localization. This interaction also results in other 
observable consequences, such as dispersion (or frequency 
shift) of the ripplon frequencies, as well as rendering the 
resonances finite width. Furthermore, we will argue, this 
interaction suffices to account for the non-universality of 
the plateau. First, however, we consider a simpler sit- 
uation, where we assume the ripplon spectrum itself is 
unaffected by coupling to the phonons. 



E. Phonon scattering off frictionless ripplons 

If e > 0, the phonon absorbing transition occurs from 
the ground state. The total number of ways to admit 
energy oj into the system is 
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/•OC 

p(w) = / tfenfc) ^ 5(« - [e + ^ n; ro w im ]) = 1/T g £ 0(u - £ »imWim)e- /s « (w - i: »» (37) 

^ {til™} im Zm 

where we sum over all occupation numbers of the ripplons with quantum numbers I, m (m = —I.. I). Using an integral 
representation of the step function 9, this can be rewritten as 

1 f°° dk 

T g /^ e 1 ^0+J_ oo 27r(ik + e 1 ) v ' 

" {n tm } 



The integral in Ea. l|55|l will be taken by the steepest de- 
scent method (SDM) . The reason why we do not apply an 
analogous technique directly to the 5-function in Eq. I|37|) 
is not only because we want to get rid of the e integra- 
tion, but also because that the SDM proved more forgiv- 
ing in terms of accuracy when used to approximate the 
^-function, rather than the 5-function. 

For each k on the real axis, the sum over the occupation 
numbers n; m diverges, so each integral should be taken 
before the summation. However, in the vicinity of the 
point that will turn out to be the saddle point k (S/cq < 
— P g ) all the sums are finite, so we reverse the order of 
summation and integration. The integration contour is 
shifted as shown in see Figll(j| 



FIG. 16 The integration contour from Ea. l)M8|l is distorted 
as explained in text. The pole k = iei is shown by a cross. 
Note the contour does not cross the pole when being shifted 
off the real axis. 
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Hence, the saddle-point approximation yields for the 
value of the integral in Eq.|j38J) (k = ik): 
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(39) 
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where the saddle point is determined from 

c = V^-^ - + - (40) 

Im 

and the curvature at the saddle point is equal to 

\f"U )\ = Y 5 ^ + \. (41) 

As is clear from (|40H , the approximation amounts to find- 
ing the effective temperature so as to populate the rip- 



plonic states to match the excitation energy to. The ex- 
pression for the curvature (|41|) appropriately involves the 
corresponding heat capacity of the excitations. 

The lu — > and the barely relevant lu — > oo asymp- 
totics are easily found. As luck has it, the lu — ► limit 
of Ea.l|39|l. apart from 1/T g factor, gives exp(l)/v27r ~ 
1.08, only 8% away from the correct 1. The ui — > oo 
yields, on the other hand, p(cu) oc Y\i m (^ / <^im) oc to 96 , 
as expected (£f =2 (2^ + 1) = 96). The SDM is thus rea- 
sonably accurate in this case, which could be at least 



31 



somewhat evaluated by computing the value of the fourth 
order term under the exponent at "one sigma" distance 
from the extremal action point. This turns out to be sat- 
isfactorily small, as demonstrated in Fig ll7l along with 
the density of states itself as a function of u. 
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FIG. 17 The solid line shows p{uj)T g from Ea.p^. The dash- 
dotted line shows the value of the fourth order term. The 
third order term, being purely imaginary, contributes only in 
the sixth order; it is shown as the dashed line, however there 
are other contributions to the 6-th order. 

When estimating absorption from the ground state, 
we totally ignore the depletion of ground state popula- 



tion at finite temperatures, when the system spends some 
time in an excited state. This is fine because by the rel- 
evant temperatures, the excited state absorption domi- 
nates anyway (see Fig^Jand note that uji — \ e \ < cj; + 1 e|). 
This case, i.e. e < 0, is somewhat less straightforward. 
Let us calculate 



N e (lu)= / den(e) ^ S(oj - n lm u lm - e]). (42) 

{";„,} Im 

This expression gives the cumulative density of absorbing 
states between energies and E (note the change of sign 
in front of e). This expression can be used to estimate 
the total excited state absorption by computing 

Pexc (u,,T) = J°° dEf{E,T)^^L, (43) 

where f(E, T) = 2/(e l3E + 1) gives the appropriate Boltz- 
mann weights. The factor of 2 is used in order to calibrate 
the excited state absorption relative to the ground state 
case: /(0, T) = 1. We now have two 6- functions and con- 
sequently two integrations. The SDM value for Ne(ui) is 
given by 



Ne(u) = , tl 1 „ |1/2 , 1 exp{(/? g + A - Mo )w + X E - ]T Iog[l - e 



T g 27r|"Det"| 1/2 AoAio 
The corresponding saddle points are determined from 

u + E = Y^ 

and 
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(44) 



(45) 



(46) 



(47) 



r 



is the determinant of the curvature tensor in the direc- 
tion (i.e. 2D subset) of the fastest descent in the 4- 
dimensional (complex) A, fi space. The steepest descent 
approximation turns out to perform well, except at very 
low frequencies (oj < lO -2 ^). However, even though it 
overestimates the answer, it is still very small compared 



to the p(u>) calculated earlier at these frequencies, much 
as the complete result would be. The appropriate graph 
is shown in Fig ll8l 

An accurate calculation of the heat conductivity re- 
quires solving a kinetic equation for the phonons cou- 
pled with the multilevel systems, which would account 
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FIG. 18 The value of the fourth order correction in the ex- 
ponent at the distance of one a from the saddle point, used 
to evaluate the adequacy of the SDM, is shown here. Each 
axes is shown in the base 10 logarithm scale. The error at 
the low energies (< 10 _2 cj_d) is extremely large, which is a 
consequence of the fact that the expression in Eq. l44H gives 
an incorrect asymptotics as u),E — > 0. However it still gives 
a density of states which is less than 1, which is all that mat- 
ters to us. The performance at the plateau energies, which is 
essential here, is good. Actually, in order to save space, we 
have presented here the data for the case with the phonon 
coupling effects on the ripplon spectrum taken into account 
{Tg/uiD = 5). This case is more interesting anyway, and the 
error here is slightly larger (but still tolerable!), we thus have 
covered all the cases. 



for thermal saturation effects etc. We encountered one 
example of such saturation in the expression l|19l) for the 
scattering strength by a two-level system, where the fac- 
tor of tanh(/?a>/2) reflected the difference between ther- 
mal populations of the two states. Neglecting these 
effects should lead to an error of the order unity for 
the thermal frequencies. Within this single relaxation 
time approximation for each phonon frequency, the Fermi 
golden rule yields for the scattering rate of a phonon with 



(48) 



The heat conductivity then equals re = 
I ^mfp(w)C w c s . The mean free path cannot be 
less than the phonon's wave-length A (which occurs at 
the Ioffe-Riegel condition). Since our theory does not 
cover the phonon localization regime we account for mul- 
tiple scattering effects by simply putting i m f p = c s t u + A. 
At high T, the heat is not carried by "ballistic" phonons, 
but rather is transfered by a random walk from site 
to sit e, as originally anticipated by Einstein l|Einsteint 
11911ft for homogeneous isotropic solids. The resultant 
heat conductivity is shown in Fig^| We postpone 

further discussion of the results above until we include 
the effects of coupling of the resonant transitions to the 
phonons on the transitions' spectrum. 
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FIG. 19 The predicted low T heat conductivity. The 
"no coupling" case neglects phonon coupling effects on 
the ripplon spec trum. The (scaled) experimental data 
are taken from JSmithl I1974T) fo r a-S iQ-2 (fcsTg/fttJD ~ 
4.4) and iFreeman and Andersonl Il986t) for poly butadiene 
(kBTg/hitiD — 2.5). The empirical univers al lower T ratio 
i m f p /Z ~ 150 llFreeman and Andersonl I1986T) . used explicitly 
here to superimpose our results on the experiment, was pre- 
dicted by the present theory earlier within a factor of order 
unity, as explained in Section IIII.BI The effects of "non- 
universality" due to the phonon coupling are explained in 
Section HVTI 



F. The effects of friction and dispersion 

A transition linearly coupled to the phonon field gra- 
dient will experience, from the perturbation theory per- 
spective, a frequency shift and a drag force owing to 
phonon emission/absorption. Here we resort to the sim- 
plest way to model these effects by assuming that our 
degree of freedom behaves like a localized boson with 
frequency u>i. The corresponding Hamiltonian reads: 



H = LUia'a 



(gk) 



-bla). 



(49) 



The ensuing equations of motion are 

(gk) 



E 



iv k b 



(gk) 



(50) 



We next introduce the following (retarded) Green's 
functions A{t) = -i9{t) ([a(t), d(0)]) and B(t) = 
— iO(t) (\b(t), al(0)]Y The fourier transforms of these 
Green's functions will consequently obey 



(u-wt)A=—+J2 

k 



(gk) 

V^kVp' 

(gk) 



--B v 



A. 



(51) 



33 



From Eans. (|51|l . one determines the real and imaginary 
parts of the Green's functions self-consistcntly. We how- 
ever can disregard the phonons' dispersion and damping 
which introduces an error in a higher order, in so far as 
the shifted frequencies w;'s are concerned. This yields 

A = — ( lo — loi % — - lim / — ; | 

2n \ 3 4:TT z pcj. £i-+o+ J co/c s + ie\ — k J 

(52) 

where k c is the cut-off wave-vector whose value will 
be discussed shortly (we have also replaced ^ k — > 

V J (2n-) 9 )' Ean. (|52*|l gives immediately for the inverse 
life-time of the internal resonance 



Uj\ 377 ■ / UJ 
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and its frequency shift 



LO < L0 C 

(53) 



LO/ 



cLuj'(uj' /lUcY 



lo' — co 



^ duj'{uj' /uj c f 



ujd J Jo uj' — lo 



where the factor of 1 /3 has disappeared because we have 
accounted for the three phonon polarizations and also ig- 
nored the distinction between the longitudinal and trans- 
verse branches. The singularity in Ea. H54|) at lo — ► ui c is 
completely artificial, as the cut-off is not supposed to 
be sharp. In our numerical estimates, we use a cut-off 
smeared by Suj c — uj c /\f~D, where D is the glass' fragility 
(see Appendix^); this is however totally unimportant as 
the divergence is only logarithmic. According to Eq. (|5Hl . 
the frequency shift scales roughly with w 3 and is thus 
rather sensitive to its value. Due to the dispersion, the 
resonance in Eq. (|52J) is effectively broadened because the 
value of the integral in Eq. (|54|l is positive for sufficiently 
small lo, but turns negative at a frequency which is a 
multiple of uo c . 

We approximate the phonon coupling effects by replac- 
ing in our spectral sums in Eas. l|39l41|) . (|44I47|) the dis- 
crete summation over different ripplon harmonics by in- 
tegration over "lorentzian" profiles: 



/ du) 5{uj — uji) — ► / dui 



(55) 



where 7w = r" 1 is a (frequency dependent) friction co- 
efficient and loi(lo) is the ripplon frequency shifted due 
to the dispersion effects. This approximation amounts to 
having the total inverse life-time of a transition involv- 
ing more than one mode being the sum of the inverse 
life-times of the participating modes. This would be cor- 
rect in the case of a frequency independent 7, but should 
be still adequate at the low T end of the plateau, where 
the absorption is mostly due to single ripplon mode pro- 
cesses. 



The value of the cut-off frequency lo c is close to but 
larger than (a/£)oj£> (see Appendix IB|l . as the phonons 
whose wave-length is shorter than £ cause an increasingly 
smaller effective gradient of the phonon field as sensed by 
a region of size £. These shorter wave-length phonons will 
still strongly interact with the droplets, however at this 
point we could only emulate that to some extent by in- 
creasing lo c . This also brings us back to the radiation 
life-time's frequency dependence. It is now clear that 
for loi(uj) > lo c , 7^ will not follow the simple cubic de- 
pendence cited above, the latter being probably still a 
safe lower estimate. We will thus use the above expres- 
sion as it makes little difference computationally in the 
region of such intense damping. At the corresponding 
temperatures, the scattering is probab ly better formally 
descr ibed by the stochastic resonance (jGamaitoni et all 
1998) methodology anyway. 

We are now ready to discuss the non-universality of 
the plateau. It is evident from Eas. (|53[) - I|54ll that even 
though the absorbers' frequencies are determined by the 
quantum energy scale lo d , the overall effective frequency 
shifts scale with T g . The ratio Tg/uoo seems to vary 
within the range of between 2 and 5 among different 
glasses, and the non-universality in this number could 
have a substantial effect subject to the value of ui c . As 
argued in Appendix iBl a value of ui c < 2.5(a/ £,)ujd is 
justified. ui c = 1.8(o/£)wd seems to yield the experi- 
mentally observed spread in the plateau's position. Our 
results for three values of T g /uJo are shown in FiglI9l 
Since lo c should be regarded as an adjustable parameter 
we can claim to possess only circumstantial evidence that 
the plateau's non-universality is caused by the spread 
in the value of the ratio of the two main energy scales 
in the problem: the classical T g and the quantum lojj. 
On a speculative note, this phenomenon may be a sign 
of strong mixing (and thus level repulsion) between the 
ripplons and the phonons, as implicitly confirmed by a 
phonon localization transition at frequencies just above 
those at the plateau. Indeed, the self-energy of an in- 
ternal resonance of dimensions £ coupled with strength 
g to an elastic medium scales (within perturbation the- 
ory) as g 2 1 pc^ 3 (x T g . This can be viewed as lowering 
of an impurity band edge due to the interaction with 
the phonons, yet another way to express the existence of 
mixing between the resonant transitions and the elastic 
waves. Within our theory, the non-universality of the 
plateau is an internally consistent proof that the degrees 
of freedom causing the Boson Peak are inelastic ones, 
whose coupling with the phonons then must be equal to 
g related to the value of T g through the stability require- 
ment explained in Section IlIII 

We now comment on the plateau slopes in Fig lI9l being 
noticeably more negative than the experimental value. 
The explaination is, we did not solve the full kinetic equa- 
tion for the interacting system, but used a simplistic sin- 
gle life-time approximation. We demonstrate this issue 
by briefly presenting a slightly different way to estimate 
P cxc (lo,T) from Eq. l|43|l . Here, we imagine we do not 
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exactly know the thermal weight function f(E,T) due 
to the lack of knowledge of the life-times in the multi- 
level system. On general grounds, however, this func- 
tion should decrease rapidly for lo > aT, where the a is 
of order unity. This yields p exc {uj,T) ~ N a T(uj) (where 
Ne (<jj) was defined in Eq. (|42|) ) . We show the result of this 
approximation for reasonable a = 1 and uj c = 2(a / '£)u; £> . 
Even though these curves resemble the experimental data 
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FIG. 20 The predicted low T heat conductivity for three 
diferent values of T 3 /ujd = 5,3.5,2 for a simpler model of 
the scattering density as explained in text. The a-SiC>2 is the 
same as in Fig |19l 



better than in the previous figure, they do not really pro- 
vide more material support for the theory than the earlier 
method. This discussion simply demonstrates that the 
basic estimates are robust enough to "survive" different 
levels of treatment. Also, curiously, these curves reflect 
the experimental tendency that the higher T plateaux 
seem to have a more negative slope as compared to the 
low T ones (see Fig^l, which was less obvious in Figll9l 

Finally in this subsection, we return to the specific 
heat. The effects of the phonon coupling on the rip- 
plon spectrum can be taken into account in the same 
fashion as in the conductivity case. Here, we replace 
the discrete summation in Eq.JJSJ) by integration over 
the broadened resonances, as prescribed by Eq.i|55fl. The 
bump, as shown in Fig ll5l is also predicted to be non- 
universal depending on T g /uiD- The predicted bump for 
Tg/ujo — 2 seems to match the best the available data 
for a-Si02, whereas the more appropriate Tg/ojo ~ 4 is 
about a factor of 3 lower in temperature. It is somewhat 
unsatisfying that the plateau's and the bump's position 
can not be thus both made to exactly match the exper- 
iment at the same time say by adjusting u>i, which is 
certainly allowed given the qualitative character of some 
of the estimates. However, since we had to employ an 
approximation when calculating the scattering density of 
states, the discrepancy does not warrant too much con- 
cern, in our opinion. 



G. The Relaxational Absorption 

In addition to the resonant absorption, an internal res- 
onance will also provide a so called "relaxational" scat- 
tering mechanism. Since a cross-over to the multilevel 
behavior of the tunneling centers leads to an increased 
resonant scattering, we must check whether the relax- 
ation mechanism is enhanced as well. This latter mech- 
anism arises because a passing phonon modifies the en- 
ergy bias of a particular pair of internal states. This 
causes irreversible thermal equilibration pr ocesses within 
each pair, resulting in energy dissipation ijjackle et all 
Il976t iMavnardl Il975j) . This phen omenon is sometimes 
referr ed to as the bulk viscosity ijLandau and Lifshit A 
Il987|h One important difference between the relaxational 
and resonant absorption is that the former does not sat- 
urate and can easily exceed the latter at low enough 
temperature and high enough sound intensity, which is 
what is usually observe d in ultrasonic experiments un- 
less s pecial care is taken l)Hunklineer and Ravchaudhuril 
1986) (this saturation is not an issue in heat conduc- 
tance, owing to the rather low sound intensities in these 
experiments). Applying the notion of the relaxational 
absorption to the two-level systems explained well the 
shape of the maximum in the temperature dependence 
of the sound speed at very low freq uencies at ~ IK 
ijHunklineer and Ravchaudhuril 11986^1 . which is one of 
the impressive achievements of the TLS model. In 
(Hun klinger and Ravchaudhuril Il986f) . the relation be- 
tween the slopes of the logarithmic temperature profiles 
around the maximum was explained. At higher T, the 
logarithmic decrease in c s is followed by what has been 
viewe d by others as a mysterious linear law l)Belessa|. 
1978). At higher frequencies still, the logarithmic de- 
crease is outweighed by the just mentioned linear T de- 
pendence. We have argued that the increase in the den- 
sity of the scattering states is due to thermal activation 
of the vibrational states of the domain walls, or match- 
ing of the thermal phonon frequency with that of a rip- 
plon on a mobile domain wall. Does the existence of the 
vibrational modes modify the relaxational scattering as 
compared to a bare underlying two-level system? The 
answer is: not significantly, for the following reason. The 
magnitude of the dissipation due to the bulk viscosity de- 
pends on the number of local distinct molecular configu- 
rations, populated according to the Boltzmann statistics. 
A shift in this population results in relaxational dissipa- 
tion. While having a domain wall excited may modify the 
energy scale in the Boltzmann distribution, which may 
produce some effect, it does not change the number of 
the intrinsic ("inelastic") glassy states, and thus will not 
on average enhance the relaxational scattering. This is to 
be compared to the resonant scattering, which depends 
on the degeneracy of the ripplon states and will thus in- 
tensify at higher T, subject to the degree of the ripplon's 
linearity. While the relaxational mechanism thus seems 
to play only a minor role in the phonon absorption at the 
plateau temperatures, its effects are observable and can 
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explain, as we will argue below, the temperature indepen- 
den t log u part in the sound speed v ariation as m e asure d 
in i|Belessai Il978l) . According to i|,Iackle et all \W7(j) . 
the variation in the speed of sound due to a collection of 
two-level systems is 
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is the radiative life-time of the ith TLS ll.Tacklel Il97l 
(see also Ea. (|53|l '). and the double angular brackets de- 
note averaging with respect to Ei, Ai and T{. While it 
would seem that detailed information on the relevant pa- 
rameters' distribution is necessary to use Eq. l|56fl . some 
qualitative conclusions can be made on general grounds. 
First, for small lj the average is dominated by the long 
life-time systems, i.e. those with A <C E and thus 
e ~ E. As a result, the averaging over these systems 
is not very sensitive to the possible correlation between 
Ei and r^, and thus the summation over the two-level 
system (nearly flat!) spectral density ~ f G?e(l/T 9 £ 3 ) 
introduces, within order unity, only a numerical factor 
proportional to T/T g (and eliminates the explicit tem- 
perature dependence). As just argued, the (e/E) 2 fac- 
tor should only give a correction factor of order unity, 
and we are left with averaging expression 1/(1 + lo 2 t 2 ) 
with respect to the life-time distribution. At low fre- 
quencies u), this averaging will be dominated by the TLS 
with the long life-times. Quite generally, for large r, 
P(r)dr tx dr/r because r _1 scales algebraically with A, 
and the distribution of log A is flat (at least for small 
A), or, almost flat, up to a weak power law, as argued 
earlier. More specifically, for a two-level system coupled 
linearly to the elastic strain, t -1 cx A 2 E fEalSTI). There- 
fore at each E (which is incidentally only weakly depen- 
dent on A in the relevant long life-time case A/E <C 1), 
obviously d(log A) = const =>■ G?(log r) = const. Thus 
the averaging w.r.t. r will produce a term of the order 
~ (g 2 / pc 2 )(l/T g £ 3 ) logw, which is of the right order of 
magnitude (and sign!). Since the dimensionless factor in 
front of the log lu term has been shown to be universal 
(cx (a/£) 3 ), the present theory predicts that it should not 
vary significantly among the insulating glasses; in fact, 
according to our argument, it is proportional to the co- 
efficient a at the logarithmic temperature dependence of 
the sound speed varianc e in the TLS reg ime, a rather uni- 
versal quantity indeed ijLeggettl Il99l[) . We stress how- 
ever that the just predicted TLS-like property should be 
observed in the plateau regime. A deviation would be a 
sign of more than two inelastic states playing a role in 
the transition. We finally mention that the lower limit 
in the integral over the life-time distribution should pro- 
duce a logT term, which would be however masked by 
the stronger linear dependence. 



V. QUANTUM EFFECTS BEYOND THE STRICT 
SEMI-CLASSICAL PICTURE 

A. Quantum Mixing of a Tunneling Center and the 
Black-Halperin Paradox 

The preceding sections have shown that structural 
transitions, accompanied at high enough temperatures 
by vibrational excitations of the mosaic, account for the 
most conspicuous departures of the low temperature be- 
havior of glasses from the prescriptions of a standard 
harmonic lattice theory - namely the existence of mul- 
tilevel intrinsic resonances in a amorphous sample made 
by quenching a supecooled liquid. At the lowest temper- 
atures these resonances behave for the most part as if 
they were two-level systems, while at higher T the den- 
sity of states of these intrinsic excitations grows consider- 
ably and leads to the Boson Peak phenomena. While we 
have computed the density of states accessible by tunnel- 
ing even at the lowest temperatures, we have assumed, 
within a semi-classical approach, that having a small tun- 
neling barrier between alternative local structural states 
does not affect significantly the corresponding spectrum 
n(e) of the lowest energy transitions from its classically 
defined value. Likewise, we have assumed that the vi- 
brational spectrum of moving domain walls is unaffected 
by the presence of tunneling, that would in principle 
mix those vibrations quantum mechanically. Clearly, the 
transitions that are active at low T must have some signif- 
icant (even if small) overlap between the wave-functions 
corresponding to the alternative structural states. This 
overlap would lead to the familiar effects of repulsion be- 
tween the semi-classically determined energy levels. This 
could be described as partial quantum melting of some 
tunneling centers, but it is probably better to use term 
"quantum mixing" . 

In this section we estimate the magnitude of these 
quantum mixing effects. Even though the strictly semi- 
classical theory agrees well with experiment as is, making 
such estimates that go beyond it is useful for two distinct 
reasons. First, we must check to what extent the semi- 
classical picture, tacitly assumed earlier, is a consistent 
zeroth order approximation to a more complete treat- 
ment. Second, it is important to ask whether the ex- 
pected corrections to the strict semiclassical theory lead 
to observable consequences. In what follows, we provide 
approximate arguments that indeed such corrections are 
discernible and may even potentially answer some long- 
standing puzzles in this field. 

Quantum Mixing. As the starting point in the discu- 
sion, we consider a simplified version of the diagram of a 
tunneling center's energy states from Fig^Jwith e < 0, 
as shown on the left hand side of Fig|2J We remind 
the reader that the e < situation, explicitly depicted in 
Fig l21l implies lower transition energies than when the 
semiclassical energy difference e > and thus dominates 
the low temperature onset of the Boson Peak and the 
plateau. 
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FIG. 21 A low energy portion of the energy level structure 
of a tunneling center is shown. Here, e < 0, which means that 
the reference, liquid, state structure is higher in energy than 
the alternative configuration available to this local region. A 
transition to the latter configuration may be accompanied by 
a distortion of the domain wall, as reflected by the band of 
higher energy states, denoted as "ripplon" states. 



Accounting for tunneling, the low energy portion of the 
Hamiltonian that corresponds to Figl^is as follows: 
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where the semi-classical values of the ripplonic energies 
are denoted as fuvi and the transition amplitudes to those 
levels are Aj respectively (only two lowest of those rip- 
plonic states are shown in Eq.JSSJ)). As argued in detail 
in Section llVl only one of the lowest two energy levels in 
a tunneling center (the top one in this case) is directly 
coupled to the higher, ripplonic, energy states. Obvi- 
ously, virtual transitions to those high energy states will 
result in lowering the energy of the higher level. There 
are no direct transitions from the bottom state of energy 
0, as explained in the previous Section, and therefore its 
position is unaffected by the presence of the ripplons. 
Consequently, the effective energy splitting of the two- 
level system (with e < 0) will be lower than the classical 
value obtained earlier, and the smaller the original value 
of e was, the more pronounced the effect will be. In what 
follows we estimate the consquences of this effect on the 
apparent energy spectrum of the lower excitations, i.e. 
the empirical two-level systems. In the limit of inifinitely 
small tunneling amplitude A, the decrease in e could be 
estimated using a perturbative expansion: 



E 



(59) 



Here, e is the new value of the energy splitting, Wj's are 
the ripplon frequencies and A, 's are tunneling amplitudes 
of transitions that excite the corresponding vibrational 
mode of the domain wall. Those amplitudes will be dis- 
cussed in due time; for now, we repeat, the expression 
above will be correct in the limit Aj/fiw.j — > 0. Fi- 
nally, the renormalized value e was used in the denomina- 
tor. While, according to Feenberg's expansion l)Feenberd . 



1948), including e in the resolvent is actually more accu- 
rate, we do it here mostly for convenience. 

Given that the semi-classical values of eigen- values hwi 
are known, the low energy portion of the energy level 
structure of the tunneling center, as shown in Fig l21l 
gives a quantitative idea of the eigen-energies of the full 
Hamiltonian only in the limit of a very small tunnel- 
ing splitting A. In a complete treatment, all transition 
amplitudes must be included and the Hamiltonian diag- 
onalized. In general, such diagonalization (and, in our 
case, the system's "quantization") is difficult, however 
could still be conducted approximately in some cases of 
interest. Consider, for the sake of argument, the follow- 
ing situation, where A is not necessarily smaller than e 
but Af/huji is. In this arrangement, the energy shift 
due to the higher lying states can be computed using 
perturbation theory and yields a "renormalized" value of 
the classical energy difference that we have called e. This 
procedure also modifies the tunneling amplitude A of the 
underlying TLS by a multiplicative factor according to 



A = A 1 
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Following this, the full energy splitting of the TLS tun- 
neling transition is computed using E = V A 2 + e 2 . The 
important feature of the argument is that A (or A) is 
allowed to take arbitrarily large values relative to e and 
the ratio of the two parameters is not treated perturba- 
tively. The lowering of e due to virtual transitions among 
the higher energy states changes somewhat the effective 
density of transition energies E that directly enters into 
the heat capacity and conductivity calculations. While 
Eq.(|Hnj) is perturbative, it should accurately give finite ef- 
fects in the mean-field limit of infinitely many transitions 
huji coupled infinitely weakly to one of the two bottom 
states of the tunneling center. We will analyze the phys- 
ical consequences assuming the accuracy of Ea. (|59f) . We 
must, of course, bear in mind that while a tunneling cen- 
ter is nearly a meanfield entity, owing to the strong cor- 
relations, it is actually of finite, albeit molecularly large 
size. Let us plot |e| as a function of |e| (see Fig l22T a)). 
Clearly, the smallest allowed value of the effective classi- 
cal splitting |e| of zero corresponds to a finite value of |e| 
equal to 
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Therefore, smaller values of |e| do not correspond to phys- 
ically realizable systems, according to Ea. (|59(l . The over- 
all excitation spectrum of the structural transitions with 
those (small) values of e < E is strongly affected by the 
ease of tunneling. As a consequence, the e and A dis- 
tributions become correlated. The quantity ^Af/hoji 
is of central importance for this section of the article, 
therefore we will discuss it now in some detail. The 
wave-functions of the highly quantum tunneling centers 
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FIG. 22 Shown in panel (a) is the relation between the bare 
energy difference e between frozen-in structural states in a 
glass and the effective splitting e that is smaller due the level 
repulsion in the tunneling center. Panel (b) depicts schemat- 
ically the derivative of e with respect to e, which is used to 
compute the new effective distribution P(e) of the transition 
energies. 



are heavily mixed combinations of the classical states 
corresponding to potential energy minima. Transitions 
between such states are strongly coupled to lattice vi- 
brations and, among other things, would strongly scat- 
ter phonons. This result is expected on rather general 
grounds and was exploited earlier when we noted that the 
eigen-states of the classical potential energy are largely 
unrelated to the eigenstates of the vibrational modes of 
the domain wall, hence one expects that transition am- 
plitudes of exciting a ripplon (which, loosely speaking, is 
a highly anharmonic combination of both structural and 
vibrational modes of the lattice) are expected to be com- 
parable to the ripplon energy itself. If this is the case, 
then number £ is actually very large relative to e and 
one, strictly speaking, should not use a perturbative ex- 
pression, such as in Eq. lf5^|) in order to assess the lowering 
of e due to quantum effects. We note, however, it is more 
likely that the seeming discrepancy simply stems from 
our "quantization" procedure being so far rather naive, 
so let us slow down a bit and attempt to outline briefly a 
more careful way to quantize the tunneling centers 's dy- 
namics. First of all, recall that we actually know the val- 
ues uji in a strongly quantum regime, because they were 
computed assuming a freely moving membrane. On the 
other hand, we know that in the classical limit domain 
wall vibrations arc indistinguishiblc from the lattice vi- 
brations. Remarkably, the vibrational eigen-frequencies 
of a "box" of dimensions £ - wn times a number from 
a/£ to 1 - span roughly the same range of energies that 
the u>iS do. Therefore, even though the quantization 
procedure will "reshuffle" all the ripplonic states, it will 
not significantly shift their position as a whole. Next, 
since static lattice inhomogeneities scatter phonons only 
elastically, the coupling of the ripplonic excitations to 
the phonons must scale with a positive power of A so 
as to vanish in the classical limit. Consequently, in the 
limit of small tunneling matrix element A, the transi- 
tion amplitudes Aj's must scale with a positive power 



of A. On the other hand, as mentioned many times, in 
the quantum regime, the Aj's are of the order hiOi and 
are not directly related to A. Therefore, a careful quan- 
tization procedure of introducing quantum tunneling in 
the system must combine and rationalize both of those 
seemingly conflicting notions. It will turn out that in the 
quantum regime, the Aj's are related to the strength of 
the underlying tunneling transition A only in a certain 
renormalized sense. We will now outline this renormal- 
ization/quantization procedure. This procedure imposes 
certain restrictions on how adding the possibility of tun- 
neling to a local structural transition can be performed, 
so that the structure of the energy levels of the transi- 
tion, that we know from general arguments, is preserved. 
Recall that "switching on" tunneling to the higher energy 
states e+hu)i not only lowered e, but also made A smaller 
by a factor of [1 — J2i / (So^) 2 ] . This expression is only 
valid if the sum is a small number, so that the whole 
correction factor is necessarily positive and only changes 
the effective magnitude of the matrix tunneling element, 
but not its sign. We must require that A not change 
its sign in the course of the "quantization" , but only its 
absolute value, because the (ordinarily small) value A 
only reflects the (ordinarily small) configurational over- 
lap between two local structural states, while the sign 
(or complex phase, in general) bears no special meaning 
here because no particular spatial symmetry is involved 
in the problem. (Such spatial symmetry is important, for 
instance, when computing overlaps between eigenstates, 
or near-eigenstates of orbital momentum centered around 
close locations in space.) Thus, the final answer should 
depend only on |A| 2 . This becomes especially evident af- 
ter the following realization. Note that the expression in 
the brackets in Ea. (|6U|l is a small coupling limit of what 
can be considered a Franck-Condon factor. The appear- 
ance of such a factor after the introduction of non-zero 
transition amplitudes is natural: the degrees of freedom 
that used to be classical and static, can now follow to 
some extent a selected motion in the system. The Franck- 
Condon (FC) factor is the overlap between the the initial 
and final wave- functions of these other ( "ripplonic" ) de- 
grees of freedom, corresponding to the initial and final 
configuration of that selected motion. (A well known ex- 
ample of a FC factor arising in an analogous, dynamical 
fashion is the tunneling m atrix element renorm alization 
in the spin-boson problem ijLeeeett et al\ . \l98l\\ .) Let us 
suppose, in a simplified manner, that the effective renor- 
malization of all of the newly introduced tunneling am- 
plitudes occurs in a similar fashion. This allows one to 
sclf-consistently close Ea. (|6U|l and rewrite it for a repre- 
sentative amplitude Aq: 
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where replacing Aq by Aq inside the brackets preserves 
the approximation's order. B is a numerical constant, 
reflecting the sum over ripplon states with their vibra- 
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tional frequencies LOi, and we have replaced LOi by the De- 
bye frequency u>r>. The two must be related at the end 
of the renormalization, as we already know. Identifying 
the expression in brackets with a Franck-Condon factor 
reminds us that Aq (or Aq) is not only a (generic) tun- 
neling amplitude but also can be considered a coupling 
constant, therefore only its absolute value is physically 
relevant in the present context, not its sign. The small- 
ness of Aq and Aq lets us recast Ea. (|62|l as 
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where the reference state is Aq = Aq = 0, 
limA Q -+o(AQ/AQ) = 1 (that is "no tunneling" =>■ "no 
renormalization"). Notice that the r.h.s. of Eq.JBSJ) 
depends explicitly only on the effective tunneling ele- 
ment Aq, but not on the original (tunable) perturbation 
strength Aq. We therefore can use the differential rela- 
tion in Ea. H63(l to extend the perturbative construction 
from Eq. (|62l) into the region of arbitrarily large values of 
the bare coupling Aq by using the outcome of the pre- 
vious (infinitesimal) change in Aq as the initial input in 
the subsequent increment g?Aq . Each of these increments 
is a small perturbation around a new, sclf-consistently 
determined, value of Aq. One gets the following self- 
consistent equation for the effective tunneling amplitude 
as a result: 



Aq = A Q e" 



(64) 



Our renormalization procedure is internally consistent in 
that the physical value of the tunneling amplitude de- 
pends on the scaling variable - the bare coupling Aq - 
only logarithmically. This bare coupling must scale with 
the only quantum scale in the problem - the Debye fre- 
quency, as pointed out yet in the first section. 

In a more complete treatment, the A renormalization 
would not be characterized by a single "Franck-Condon" 
parameter, but by a distribution of Franck-Condon fac- 
tors. Therefore, the exponential form in Ea. (|64|l might 
be possibly replaced by a different, perhaps a polynomial 
expression. In fact, one may think minimally of Eq.JfJlJ 
as of one of the possible Pade extensions of the perturba- 
tive formula At any rate, such a Pade approximant 
will retain the main feature of Eq. (|64|l in that the value of 
the observable tunneling matrix clement Aq is bounded 
from above and depends strongly on the (semi-)classical 
energies tkoi. According to the discussion above, this re- 
striction stems from a self-consistency condition, namely 
that the leading term in the exponent in Eq. l|64[) must 
scale with A 2 ,, if this same number Aj~ is on the l.h.s. 
in that equation. An important corollary of this is that 
the perturbative term inside the brackets of Eq. (|60|l must 
scale with A 2 itself, hence the perturbation correction £ 
will scale with A 2 too (from now on, we will drop tildes 



from the symbols denoting the physical tunneling ampli- 
tudes, but retain them for the effective e's). That is, with 
our definition of B, it is roughly true that 
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We have thus demonstrated explicitly that the magnitude 
of quantum effects on the classical energy splitting e on 
a particular site should depend on the facility of tunnel- 
ing at that same site. We have therefore established that 
the fact of quantum A 4 being close in value to a rather 
large energy scale Ui is consistent with a relatively small 
value of the correction in Ea. (|59|) and its scaling with 
A 2 . As another dividend from the argument, we obtain 
a ballpark estimate of the constant B. A structural tran- 
sition that is thermally active at a plateau temeperature 
ksT ~ hbji and being an efficient phonon scatterer, will 
have A ~ hiOi. Therefore, B = ^iLOjj/LOi, will be a 
number on the order of several hundreds, since the total 
number of the ripplonic modes (at the laboratory glass 
transition, at which £/a ~ 5.. 6) is approximately a hun- 
dred, and LOi is proportional to, but somewhat smaller 
than the Debye frequency lob- 

Eq. H65fl implies that while the distributions of the cla- 
sical energy splittings e and the bare semiclassical tun- 
neling amplitude A may be uncorrelated, quantum cor- 
rections require that e and A be correlated for systems 
with a sufficiently low barriers and which simultaneously 
have small energy difference between the initial and final 
structural state. Conversely, the independence approxi- 
mation is valid when, roughly, |e| > B-^-^. Since A is 
proportional to Hlod, this criterion is a formal restate- 
ment of an earlier comment that the theory is strictly 
valid in the classical limit. (Note that there is also a 
(much stronger) h dependence in the exponent of the tun- 
neling element A (see Eq. (12111 ). While, obviously, only 
a negligible fraction of the total number of structural re- 
arrangements in the liquid at T g would not satisfy the 
classicality criterion, these particularly facile transitions 
do actually comprise a significant portion of those transi- 
tions that are thermally active at cryogenic temperatures. 
We will now indicate what the observable consequences 
of this deviation for the strict semiclassical limit are. In 
order to do this, let us discuss first the difference between 
the strongly quantum and the bare "classical" structural 
transitions. 

According to Eq.j55J), for all transitions, whose diag- 
onal energy difference would be |e| < S in the classical 
limit, the effective diagonal splitting e is actually zero, 
meaning that the full energy splitting E is entirely com- 
prised of the originally off-diagonal energy scale A. This 
implies that the energy eigen-states of such highly mixed 
tunneling centers are heavy superpositions of the original 
classical structural states and would not be easily inter- 
preted in terms of the atomic coordinates of the potential 
minima alone, but must include the kinetic energy term 
as well. This is directly related to the well known ambigu- 
ity in separating the energy of such systems into potential 
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and kinetic components even at conditions that are en- 
tirely classical, such as at a high temperature. Of course, 
in such cases free energy formulations must be employed 
that allow one to count the number of configurational 
states unambiguously, while using "inherent" structures 
based on potential energy stationary points alone is of 
limited utility. The strongly quantum case can be loosely 
understood by transcribing the complex multiparticle re- 
arrangements onto a single collective "reaction" coor- 
dinat e (as in the soft potential model ijGalnerin et all 
Il99l|) h In fact, this analogy to a single coordinate soft 
potential model is quite loose because of the much higher 
density of states of the ripplons (that give rise to the Bo- 
son Peak and correspond to the vibrations of the mem- 
brane) compared with the density of states of the soft 
potential model, which is one dimensional so that only 
one coordinate is vibrationally excited. Nevertheless, fol- 
lowing this analogy, consider a two-well potential (with 
very steep outer walls) with a barrier high enough so that 
the physical coordinate eigen-states corresponding to the 
particle being in the left or the right well are unambigu- 
ously definable and the diagonal component of the tran- 
sition's energy is equal to the difference in the potential 
energy of the two well with high accuracy. Imagine next 
lowering the barrier. In the limit of zero barrier the sys- 
tem is simply a particle in a square box, whose energy 
scale is determined by the quantum energy scale in the 
problem - that is the particle's kinetic energy alone. This 
analogy reminds us that just like the transition from a 
largely classical to quantum behavior in a double well 
potential, the transition at |e| = S is not sharp (note, 
however, that unlike in a one dimensional soft-potential 
model, the density of excited states of a tunneling center 
is very high thus possibly leading to a sharper cross-over). 
Put another way, this "phase transition" clearly corre- 
sponds to term-crossing and therefore would be gradual 
in a finite system. From a mean-field perspective, the 
transition at e = £ resembles a de-localizatio n phase 
transition (see e.g. ijAbou-Chacra et all Il973|) ') which 
we may think of as quantum depinning of the domain 
wall. Alternatively, one could say that the local struc- 
ture of classical energy levels melts out locally in that 
the energy variation on the mostly classical landscape 
(determined by T g ) happens locally to be smaller than 
the confinement kinetic energy of the domain wall mo- 
tion. Of course, this is occuring for only small parts of 
an otherwise rigid matrix. Again, since the system is fi- 
nite, one expects a soft cross-over rather than a sharp 
transition when such "melting" occurs. Both ways of in- 
terpreting the quantum mixing/melting described above 
are consistent with our view of the tunneling process lead- 
ing to the expression (|21fl for the tunneling amplitude A. 
The action exponent in Ea. l|21[l scales as the height of 
the barrier relative to the under-barrier frequency. The 
former quantity, while distributed, scales with the classi- 
cal energy scale in the problem - T gi while the latter is 
proportional to the Debye frequency (and, most likely, is 
somewhat distributed too). The quantum limit of large 



hwD corresponds to a narrow barrier and a short tunnel- 
ing path. This would imply the relative unimportance 
of the classical energy landscape modulation during the 
tunneling process. Finally, in order to avoid ambiguity, 
we stress that the structural transitions of both types of 
tunneling centers, that we have called "classical" (in that 
the wave functions are well localized near minima and 
are well defined structurally, i.e. in position space) and 
"quantum" (i.e. in a superposition of structural states), 
at low temperature occur in a purely quantum mechani- 
cal fashion, that is by tunneling. 

We now show that the presence of a somewhat distinct 
class of such low barrier, or "fast" , two-level systems, 
whose effective diagonal splitting is zero, leads to addi- 
tional phonon scattering in comparison with the strictly 
semiclassical analysis, which neglects the renormalization 
from quantum mixing effects. This additional scatter- 
ing at low energy is consistent with the apparent sub- 
quadratic temperature dependence of the heat conductiv- 
ity in the TLS regime. The mixing also leads to a super- 
linear addition to the heat capacity at subKelvin tem- 
peratures. These highly quantum tunneling centers in 
strongly mixed superpositions of structural states, there- 
fore, give a mechanism to resolve a quantitative deviation 
from the standard tunnelin g model, which was brou ght 
up by Black and Halperin ijBlack and Halperinl Il977|) in 
1977. They noted that the short time heat capacity of a- 
Si02 is larger than would be predicted by the logarithmic 
dependence obtained in the STM, if one uses the TLS pa- 
rameters extracted from ultra-sonic measurements. The 
quantitative mismatch appears to be as if there were two 
kinds of two-level systems: one set obeying the distribu- 
tion postulated in STM, and another set of "fast" tun- 
neling centers responsible for the short time value of the 
heat capacity. We can see our analysis of mode mixing 
leading to the existence of a finite number of two-level 
systems with e very nearly 0, as suggested by Ea . (|59fl is 
quite consistent with this empirical notion 6 . 

To see this more explicitly we note that Eq. Ij59(l allows 
one to formulate the effects of quantum mode mixing as 
a change in the apparent distribution of the diagonal en- 
ergy splitting. Whatever the old distribution of classical 
energy difference n(e), the new distribution of the effec- 
tive classical component of the transition energy can be 
found using n(e)|de| = n(e)|de|. For e's not too close to 
huji (case e <~ Hu>i will be discussed later), which is ap- 
propriate in the TLS regime, the function \de/de\ that 
describes the relative probability distribution of the two 



We must stress however that the Black-Halperin analysis has 
been conducted only for a single substance, namely amorphous 
silica, and systematic studies on other materials should be done. 
The discovered numerical inconsistency may well turn out to 
be within the deviations of the heat capacity and conductivity 
from the strict linear and quadratic laws repsectively. Finally, 
a controllable kinetic treatment of a time-dependent experiment 
would be necessary. 
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quantities, is given by 
de 



(66) 



where the <5-function is positioned to the right of the ori- 
gin: f deS(e) — 1 (see also FigE2b)- Consequently, the 
distribution of the effective diagonal splitting is: 
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The coefficient of the <5-function reflects the "pile-up" 
of the two-level systems that would have had a value 
of |e| < S were it not for quantum effects. These fast 
two level systems will contribute to short time value of 
the heat capacity in glasses. The precise distribution in 
Eq. (I67|l was only derived within perturbation theory and 
so is expected to provide only a crude description of the 
interplay of clasical and quantum effects in forming low 
barrier TLS. Quantitative discrepancies from the simple 
perturbative distribution may be expected owing to the 
finite size of a tunneling mosaic cell, as mentioned ear- 
lier, and the finite life-times of each energy state due to 
phonon emission. These effects would also smoothen the 
local quantum melting transition as e — > 0. While vari- 
ous improvements of the functional form of n(e) might be 
suggested, it seems unwarranted, at present, to use any 
more complicated expressions for this function. Thus, 
to see the main consequences of the quantum mixing ef- 
fect, we will proceed with the perturbative expression. 
Assuming a particular value of the coefficient B allows 
one to derive the contribution of the fast two-level sys- 
tems to the heat capacity and scattering of the ther- 
mal phonons. Before we start, let us note that since 
we now have to deal with a specific coupled distribution 
of e and A, the generic two- level system model that only 
specifies the distribution of the total splitting E is not 
sufficient. We must use the full tunneling model where 
the tunneling elements A's are distributed according to 
Eq. (1221 . The exact value of constant A in equation 1221) 
depends (weakly!) on the (possibly e-dependent) cut-off 
value of the P{A) distribution. Both the heat capacity 
and the phonon scattering strength depend on the co- 
efficient A, therefore it is possible to check the relative 
contribution of the "quantum" centers to both of those 
quantities, regardless of A's value. The n(e, A) distribu- 
tion obtained in this way is now a product of the -P(A) 
distribution from Eq.()22p and the density of states from 
Eq.lJHJJ. The new normalization coefficient A% is found 
from the requirement that f dedAn(e, A) = l/£ 3 . This 
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der to compute the life-time of a phonon of energy E, one 

2 a 2 io rp 

averages the Golden Rule scattering rate 9 ^ E tanh ^- 
with respect to n(e, A), subject to the resonance condi- 
tion E = \J e 2 + A 2 l)Anderson et al l ll972H.Iacklelll972t 
Phillips, 1981). This yields two contributions to the de- 
cay rate: 
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The first term in the square brackets is the contribution 
owing to the fast, or highly quantum, two-level systems. 
Note that this term scales faster with E than the other 
term. Provided the magnitude of this first term is com- 
parable to the other term, the fast modes will somewhat 
modify the overall scaling of the heat conductivity n. 
Without the first term, k scales superquadratically ac- 
cording to T 2+c (recall that the heat conductivity is in- 
versely proportional to the scattering rate from Eq. (|68(l ). 
If we use a numerical value of B of the order 100, this 
leads to a subquadratic T dependence of k: Experimen- 
tally, k(T) scales like T 19±1 as extracted from a decade 
and a half of data (see Fig|I}. Without the fast TLS, 
one, again, would have k oc T 2+c . Using the theoreti- 
cal approximation for c, this differs from the empirically 
observed value at least by a factor of (IO 1 " 5 ) " 1 "' 1 ~ 2 at 
T ~ 10~ 2 Td. Obviously, this is a very crude estimate be- 
cause, first, we do not know how far down in temperature 
the power law scaling of K goes; second, our correction, 
while going in the right direction, summed with the older 
result, is not strictly a power law. Since the integral in 
the square brackets of Ea. H68(l varies between 1 and tt/2 
for < c < 1 (A min /E < 1, surely at E ~ 10~ 2 T D ), 
we conclude that the first term must be between 10° 
and 10 1 in order to make a sizable contribution to the 
phonon scattering and modify its functional form. Since 
E ~ I0~ 2 Td, this shows that B indeed must be of the 
order of several hundreds, consistent with our expecta- 
tions based on the number of vibrational modes in the 
Boson Peak. 

Does this mixing induced correction to the density of 
states with the value of B around a hundred make an 
appreciable contribution to the time-dependent heat ca- 
pacity? Following the calculation from subsection llll.CI 
but now using the new distribution n(e, A), one finds: 
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where 0(t) is the usual step-function and r m i„ is the 
fastest possible relaxation time of a TLS with the to- 
tal energy splitting E, defined in Ea. 1)25(1 . Again, the 
first term in the square brackets gives the contribution 
of the "fast" TLS. Using the same numbers as given in 
subsection IIII.CI it is straightforward to show that the 
second, regular, term is of the order a hundred at tem- 
peratures T ~ 10 _2 Td when measured on the time scale 
of minutes. At the same time, the first term is at most of 
order ten. Note that at the shortest times t ~ T m i n , when 
the regular two-level system only begin to contribute to 
the heat capacity, the theory with quantum corrections 
says the actual heat capacity is finite and is at the most 
one tenth of the long-time value. At the same time, the 
fast tunneling centers do not seem to contribute signifi- 
cantly to the long-time heat capacity. We note however 
that the result obtained c(T) oc T 1+c / 2 with c = .1 gives 
a somewhat slower rise with temperature than seen in 
experiment. The quantum correction again goes in the 
right direction of increasing the rate of the heat capacity 
growth with temperature relative to the T 1+c / 2 law. 

We have established that effects beyond the strict semi- 
classical analysis give rise to a subset of tunneling centers 
that undergo faster tunneling than the rest. Neverthe- 
less, there are some quantitative issues in the heat capac- 
ity magnitude that remain to be understood, namely that 
the computed contribution of the "fast" centers seems 
somewhat lower than what is necessary to explain the 
deviation of the experimental T dependence from the 
supelinear dependence T 1+c / 2 predicted by the present 
(approximate) argument. It is posssible that ultimately 
a broader view of the time-dependence of the heat ca- 
pacity needs to be taken. Since, in fact, the system will 
clearly be aging by tunneling at those low temperatures, 
the notion of fixed frozen-in "defects" may no longer be 
adequate - essentially interactions between defects play 
a role. "Aging" by definition implies irreversible struc- 
tural changes. More work on understanding the long time 
evolution of the tunneling centers is necessary. 

We have concentrated on the quantum corrections to 
the low lying tunneling states with low barriers. Quan- 
tum mixing applies to the higher energy states too. En- 
ergy shifts and quantum melting occur within sub-bands 
of the ripplonic states of order I and respective degener- 
acy (2Z + 1), thus mixing these states. As tunneling can 
take place on a given time scale and the vibrationally 
excited levels become observable, their apparent ener- 
gies can not be degenerate because the levels are coupled 
through those same tunneling transitions. The magni- 
tude of energy level repulsion from the quantum mixing 
can be assessed qualitatively. In the limit of weak cou- 
pling, the deviation of a ripplonic frequency from its clas- 
sical value scales J2i Af/fojj. The width of the ripplonic 
band of order I is probably limited from above by the 
tunneling amplitude Aj itself. Does this band broaden- 
ing affect our previous results on the Boson Peak phe- 
nomena? Not very much. Since the observables depend 
mostly on the number of new excitations and the number 



of the ripplonic modes is not changed by these mixing ef- 
fects, the essential core of our conclusions from Section 
IIVI remains intact. Nevertheless, some quantitative mod- 
ifications are to be expected. For example, the lowest 
ripplonic energies may be lowered to the extent so as to 
cause a cross-over to a multi-level behavior in some of the 
internal resonances, thus possibly modifying the derived 
magnitude of the heat capacity and phonon scattering at 
sub-plateau temperatures. This effect will further con- 
tribute to the phonon interaction induced broadening of 
the ripplonic transitions, as estimated in Section llVl 



B. Mosaic Stiffening and Temperature Evolution of the 
Boson Peak 

Ea. l(59|l raises another interesting point. According to 
that equation, the values of both the bare and the ef- 
fective classical energy bias of a transition - e and e re- 
spectively - are limited from above by the lowest ripplon 
frequency (0*2)- (Note that this is only realized in the 
e < case, discussed in this section.) This is unimpor- 
tant at low temperatures. But what happens at higher 
T, near this limit? Unlike in the low energy situation 
just discussed, one simply cannot ignore here that all the 
energy states have a rather short life-time. Therefore 
the singularity in Eq.JSHJ does not occur, but will be 
rounded. This observation does not completely answer 
the question that one should have asked in the first place 
on general grounds alone: what happens to the structure 
of the energy spectrum of a tunneling center, when the 
energy of the transition becomes comparable to a vibra- 
tional eigen- frequency of the domain wall 7 ? 

When attempting to answer this question, a general 
multi-level perspective on each tunneling center is some- 
what easier to use than the very mechanical view of the 
wall's excitations that we have mostly employed so far, in 
which the ripplonic energy states are obtained by quan- 
tizing vibrations of a freely moving classical mambrane. 
The "singularity" at |e| ~ hwi is actually a term-crossing 
phenomenon that, again, would not take place in the 
strict classical limit. Let us go back to our argument on 
the density of states, but consider a case when e is larger 
than a ripplonic frequency. As mentioned many times 
already, vibrational excitations of a domain wall can be 
defined meaningfully only when a structural transition 
takes place in a given region of the material. The energy 
of the transition must be the lowest excited state of a 
mosaic cell. On the other hand, the values of the ripplon 
frequencies are determined by a (fixed) surface tension 
coefficient and the wall's mass density. They have fixed 
values. The necessary conclusion from this is that the 



We remind the reader that the tunneling transition energy could 
be also thought of as an eigen-energy of the wall's motion, but 
of a lower, 1 = 1 order, associated with the translational motion 
of the shell's center of mass 
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tunneling centers will not have ripplons whose frequency 
is lower than the transition frequency. We provide a car- 
toon illustrating this idea in FigEU We see the quantum 
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FIG. 23 This caricature demonstrates the predicted phenom- 
ena of energy level crossing in domains whose energy bias is 
comparable or larger than the vibronic frequency of the do- 
main wall distortions. The vertical axis is the energy mea- 
sured from the bottom state; the horisontal axis denotes tem- 
perature. The diagonal dashed line denotes roughly the ther- 
mal energies. A tunneling center, that would become ther- 
mally active at some temperature To will not possess ripplons 
whose frequency is less than To- 



mixing reduces the number of the lower frequency vibra- 
tional modes. The mosaic appears stiffer than expected. 
This effect may contribute to the temperature evolution 
of the Boson Peak as observed in in elastic scattering ex- 
perim ents. Wischnewski et al. find (Wis chnewski et all 
1998) that at temperatures between 51 K (numerically 
close to silica's hu>2) and above the glass transition, the 
left hand side of the Boson Peak decreases in size as the 
temperature was raised. At the same time, the high fre- 
quency side remained relatively unchanged. Note that, as 
temper ature is raised, the total a rea of the peak in Fig. 2 
of Ref . ijWischnewski et all . Il998|) does not increase. In 
this temperature range mosaic cell motion loses oscillat- 
ing character and becomes a rather featureless activated 
relaxation process. 

To summarize this section, we have seen that the pos- 
sibility of quantum tunneling between structurally close 
states in glass does have a predictable effect on the spec- 
trum and must be taken into account when computing 
the density of low (and not so low) energy structural 
excitations in these materials. At the same time, the 
main conclusions of the original semi-classical argument 
remain valid: each structural transition may be thought 
of as a rearrangement of about 200 molecules accompa- 
nied by distortion of the domain wall that separates the 
two alternative local atomic arragements. 



VI. THE NEGATIVE GRUNEISEN PARAMETER: AN 
ELASTIC CASIMIR EFFECT? 

With the exception of the plateau's position and the 
quantum mixing effects, we have so far dealt with those 
anomalies in low temperature glasses that are more or 
less universal. These universal patterns are of particular 
interest because they cannot be easily blamed on chemi- 
cal peculiarities of each substance. Indeed, given the flat- 
ness of the low energy excitation spectrum in glasses, the 
apparent universal ratio l m f p /\ ~ 150 is the dimension- 
less quantity that seems to express the general, intrin- 
sic character of those low energy excitations, as arising 
from the non-equilibrium nature of the glass transition. 
The number 150 reflects the size of nearly independent 
fragments into which a supecooled liquid is broken up 
at the laboratory glass transition. Yet, there is another 
dimensionless quantity, namely the Gruneisen parame- 
ter 7, that also reflects the necessity of going beyond a 
harmonic picture for amorphous solids. This parameter 
is always a positive number of order one for simple cu- 
bic crystals (at low e nough T), but varies wi ldly among 
amorphous m aterials llAckerman et all . Il984|) (see also a 
discussion in l|Leggettj ll99ll) V -y in glasses has been re- 
ported to be as large as several tens and often negative 
in sign! A negative 7 implies a negative thermal expan- 
sion coefficient y (the linear expansion coefficient 

a = £ (fy) p = Ip 5 1S a commonly used quan- 

tity chracterizing anharmonicity too). Contraction with 
heating is observed in some crystals at not too low tem- 
peratures, owing to the details of the anharmonic cou- 
plings in a specific substance that may result in the neg- 
ativity of the Gruneisen para meter of a, lattice mode of 
a finite frequency (see e.g. ijWei et all Il994j) 1. Ther- 
mal contraction along a single direction in anisotropic 
materials is even more common. Nonetheless, as the 
temperature is lowered, the thermal expansion coefficient 
in an insulating crystal eventually becomes positive and 
approaches the cubic T dependence predicted by stan- 
dard thermodynamics. In contrast, an isotropic negative 
thermal expansivity is observed in many amorphous sub- 
stances even at the lowest temperatures. In addition, the 
expansivity is not cubic in T. The most widely known 
example of a substance with a negative a is rubber. Rub- 
ber owes this property to the largely entropic nature of 
its elasticity. Here, we will see that a distinct mechanism 
of thermal contraction in glasses in the TLS temperature 
range arises, which is a direct consequence of the ex- 
istence of the spatially extended tunneling centers that 
give rise to the universal phenomena considered earlier. 

As shown above, the excitation spectrum of the tunnel- 
ing centers may be represented as a combination of the 
two lowest energy levels, corresponding to the structural 
transition and a set of higher energy states involving vi- 
brations of the moving domain wall. By the exchange of 
phonons, these local (quantum) fluctuations in the elastic 
stress will be attracted to each other much like in the Van 
der Waals interaction between neutral molecules. The 
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elastic Casimir effect seems a more appropriate name for 
this phenomenon, since the moving domain walls are not 
point-like but, instead, resemble fluctuating membranes. 
While we do not claim this attraction is solely responsi- 
ble for the negative expansion coefficient, it turns out to 
provide a large contribution to the thermal contraction 
in glasses. We will see how this effect arising from in- 
teraction of amorphous state excitations depends on the 
material constants and the preparation speed of the glass 
is derived and, therefore, is not universal! 

We note first that not all amorphous substances ac- 
tually exhibit a negative a in the experimentally probed 
temperature range. In such cases, it is likely that the con- 
traction coming from those interactions in these materi- 
als is simply weaker than the regular, anharmonic lattice 
thermal expansion. Other contributions to the Griineisen 
parameter will be discussed below as well. 

Coupling the motion of the mosaic cell (TLS and Bo- 
son Peak) to phonons is necesssary to explain thermal 
conductivity, therefore the interaction effects discussed 
below follow from our identification of the origin of amor- 
phous state excitations. The emission of a phonon fol- 
lowed by its absorption by another cell will give an effec- 
tive interaction, in the same way that photon exchange 
leads to inter-particle interactions in QED. The longest 
range coupling between local degrees of freedom coupled 
linearly to the elastic stress has the form of a dipole- 
dipole interaction. Since the structural transitions are 
of finite size, the dipole assumption is only approximate 
for the closer centers. For the time being, we take for 
granted that there is no first order, static, interaction 
between the vibrating domain walls, which, if non-zero, 
could be a priori of either sign. The next, second order 
interaction is always negative in sign and is proportional 

to - Ey ^ « " (1 - IT 1 )' - -1 + 2 TT- Thi s favors a 
sample's contraction (V is the volume). This attractive 
force, which will be temperature dependent, is balanced 
by the regular temperature independent elastic energy of 
the lattice: F e i ast /V = y (^) . Calculating the equi- 
librium volume from this balance allows us to estimate 
the thermal expansion coefficient a. More specifically, 
the simplest Hamiltonian describing two local resonances 
that interact off-diagonally is H — ^a^ + ^-al + Jijalai, 
where u>i and ujj would be the frequencies of ripplons on 
sites i and j and 

, 3 (K,S,i 3lK,i\;KS,i V 

13 - tope* r% ( W) 

is the dipole-dipole interaction following from Eas. (|15|) 
and (|16|) . (Having the interaction be off-diagonal auto- 
matically removes the first order term in J,j .) The factor 
3 accounts in our usual simplistic way for all three acous- 
tic phonon branches. This ignores a distinction between 
the longitudinal and transverse speed of sound. This sim- 
plification is however accurate enough for our purposes. 
Since g ~ yf pcja 3 kBT g , the Jy's turn out to scale in 
a very simple way with the glass transition temperature 



and the molecular size a, giving ~ ksT g (^) . 

Since only mobile domain walls give rise to local dy- 
namic heterogeneities, one may conclude intuitively that 
only the sites of thermally active structural transitions 
can contribute to a. Therefore one expects that as tem- 
perature in increased, more tunneling centers will con- 
tribute to the Van der Waals attraction thus leading to 
negative expansivity. As already mentioned, the excita- 
tions of a tunneling centers are conveniently subdivided 
into a low energy TLS-like pair of states, and higher 
energy, "ripplonic" excitations corresponding to distor- 
tions of an active center's domain wall. Hence we may 
view the total Van der Waals attraction as having three 
somewhat distinct contributions: "TLS-TLS", "ripplon- 
ripplon" and "TLS-ripplon" attractions. In this section, 
we focus on the relatively low, subplateau temperature 
regime, for reasons that will be explained later. At these 
low temperatures, transitions to the ripplonic states are 
only virtual, whereas the TLS structural may well be 
thermally active. This, in addition to the differences in 
the respective spectra of these excitations, will lead to 
some difference in the dependence of the mutual inter- 
actions between those excitations on temperature and 
other parameters. In order to assess the magnitude of 
those interactions let us consider the following, very sim- 
ple, three-level Hamiltonian that is designed to model a 
transition of energy between two different structures 
that may also be accompanied by a wall vibration of fre- 
quency LOf. 

/o \ 
Ht = ei (71) 
\0 e l +uj, J 

Note that, even though, for simplicity's sake, we use the 
semiclassical energy e in the Hamiltonian above, the lat- 
ter is meant as (the lowest energy portion of) the full, 
diagonalizcd Hamiltonian with quantum corrections in- 
cluded. This corresponds to the plain two-level system 
formalism that does not specify a distribution of the tun- 
neling matrix element A. Also, in comparison with the 
general case of Eq. I|35|l , we only include an excitation by 
a single quantum of a single ripplon. The latter simplifi- 
cation is obviously justified in the lowest perturbation or- 
der, where all pairs of excitations contribute to the total 
in an additive fashion. Considering only single-quantum 
excitations is a low temperature approximation, made 
mostly to avoid adopting extra modelling assumptions 
necessary to embody the mixed spin/boson statistics on 
each site. This simplification is nevertheless adequate, as 
will become clear later in the discussion. 

Since the contributions of the three constituents of the 
Van der Waals attraction are additive, one can consider 
each contribution separately. This indeed proves to be 
convenient not only because all the contributions exhibit 
distinct scaling with the parameters, but each contribu- 
tion comes to dominate the expansivity at somewhat dis- 
tinct temperatures. We consider first the ripplon-ripplon 
attraction. This contribution appears to dominate the 
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most studied region around 1 K. The off-diagonal (flip- 
flop) interaction between the ripplons has the form: 



H$* = J ij \2 i 3 j )(3 i 2 j \ + H.C.. 



(72) 



where the rows and columns in the unperturbed Hamilto- 
nian from Eq. I|7f(l are numbered in the conventional way 
from the upper left corner. The "ripplon-ripplon" case 
appears the simplest of the three because here, the issue 
of how many tunneling centers contribute to the effect is 
more or less separate from the strength of the interaction. 
The former is (qualitatively) determined by the number 
of thermally active two-level systems, that scales roughly 
with the heat capacity. The latter is nothing but the 
ground state lowering of a pair of resonances after inter- 
action is switched on, which scales as —Jfj/ (ui + Uj) and 
is T-independent at these low temperatures. This con- 
tribution to the negative thermal expansion is therefore 
expected to be roughly quadratic in temperature (this 
corresponds to linear expansivity), which is similar, if 
not somewhat slower than observed in amorphous silica 
around 1 degree K. 

Calculating the correction to the system's free energy 
in the lowest order in J^-, that corresponds to the inter- 
action term from Eq.JT^J) is entirely straightforward and 
yields: 



SFrr = -J24 



-/»(et+<y)(l + e -^')(l + e"^) 
ZiZj 



LUi tanh(/3ciJi/2) — tUj tanh(/3wj/2) 

X 7i o , \'<Jj 



Here, Zi = 1 + e _/3ei + e~^ ei+u}i ^ is the unperturbed on- 
site partition function, corresponding to Ea. (|7f fl . Here, 
subscript "rr" signifies the "ripplon-ripplon" contribu- 
tion. 

At low - subplateau - temperatures T < a;,, that we are 
primarily interested in here, the expression above reduces 
to the following Van der Waals energy: 



J 2 



1 



1*2 L l l 2 



(74) 



where we have explicitly written out summation over dis- 
tinct ripplon harmonics h and h at sites i and j. 

A few intermediate calculations are needed to com- 
pute the sum in Ea.Q74[l. First, averaging of Jfj with 
respect to different mutual orientations of gi, g^ and 
Vij yields an effective isotropic attractive interaction 

I (^) 2 Tg (^) 6 - Second, the sum over all harmonics 
amounts to vj m " x (2Zi+i)(2Z 2 +i) w j lere w , j s f ounc i us _ 

ing the dispersion relation from Eq. (|33[1 . Here we assume 
that u>iS are not correlated with Jij and e 2 ;. As we al- 
ready know from the discussion in the previous section, 
the latter assumption is adequate for values e smaller 
than the Boson Peak frequency. Now, recall that l max 
actually depends on the droplet's perimeter, thus intro- 
ducing an additional (cubic!) scaling with £/a. In the 



end, the sum over the l's is equal, within sufficient accu- 
racy, to 1.5w~ 1 (3/47r)7r 3 (^/a) 5 / 4 (^/a) 3 . Finally, assum- 
ing Jij's and e's to be uncorrelated enables one to present 
the double sum over ej as a product of two identical sums: 
+ /3 £ i)~ 1 ) ■ Each sum is the effective concentra- 
tion of thermally active tunneling centers: /cs(ln2)y^y 
as computed by integrating 1/(1 + e@ e ) with the density 



of states from Eq. 

C3 



Note that here we use the sim- 



ple 1 /Tgt; expression for density of the tunneling transi- 
tions, in keeping with the assumption E ~ ej of the plain 
two-level system model adopted in this section. This is 
reasonable, given the qualitative character of this calcu- 
lation. Finally, the summation over the ripplon sites can 
now be reduced to an integration with the lower limit 
equal to £(3 /An) 1 / 3 . 

As a result of the previous discussion, one recov- 
ers the following expression for the energy gain (per 
volume) due to a volume change SV: SF rr /V ~ 

,7/4 

'tjH. This works against the 

regular elastic energy 5F e i ast /V = -j (ir) 2 ; introduced 
earlier. The equilibrium relative change 5V/V as a func- 
tion of T is obtained by setting dF/dV = 0. Differen- 
tiating the equilibrium value of 8V with respect to tem- 
perature yields the following estimate for the thermal 
(volume) expansion coefficient: 



1.5(ln2)Vfjg(f)' 



V\dTj v y ' K?T D \£ 



7/4 



(75) 



This can already be used to estimate the magnitude of 
the ripplon-ripplon contribution to the "Casimir" effect 
numerically. One can do it in several ways. The simplest 
thing to do that does not require knowing K , is simply 
to use Ea. (|75fl to calculate the Griin eisen parame ter 7 
itself according to 7 = (dp/dT)v/cy l)Kitte]tll956|) . also 
using (d P /dT) v = -(dp/dV) T /(dT/dV) p . This yields a 
temperature independent Griineisen parameter: 



7rr c-3.0(ln2) 2 7 r 2 ^(| 



7/4 



(76) 



Using (£/a) 3 ~ 200 and silica's T g /T D ~ 1500/350 
one obtains 7 ~ —3., within an order of magnitude of 
what is observed in amorphous silica at low tempera- 
tures (that experimental number varies between —5 and 
—20 among different kinds of silica at 1 K and seems to 
grow larger with lowerin g the temperature, see Fig. 3 from 
l)Ackerman et alj\l9 84)). We will argue shortly that this 
growth may be explained by other contributions to the 
attraction between local resonances. 

We can also directly compare the conribution in 
Eq. (|75|) to the linear thermal expansio n coefficient a = 
W (f rOn f° r smca as measured in ijAckerman et all . 
19^). According to the Fig. 2 from l|Ackerman et all 
1984^) . the a of silica is linear (possibly slightly sub- 
linear) in temperature and equals —1.0 ■ lQ^K^ 1 at IK. 
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The c ompressibility K was obtained in ijAckerman et all 
Il984|) from measured speed of sound and density. For in- 
ternal consistency, we use the scalar elasticity to estimate 
K in this way. Summing up three single polarization 
phonon Hamiltonians from Eq. (|15fl yields K ~ pc 2 s /3 (re- 
member, AV/V = 3A0). Using silica's constants, given 
in Fig^l and the earlier obtained £ = 20 A and recall- 
ing that 81/1 — SV/3V, Eq.J75J) gives linear expansion 
coefficient a ~ —0.4 ■ lO -9 ^ -1 at IK, indeed strongly 
suggesting that attraction between the tunneling centers 
is a significant contributor to the negativity of the ex- 
pansion coefficient. The numbers just obtained are also 
a convenient benchmark in assessing other contributions 
to the negative thermal expansivity. 

Next, we estimate the magnitude of the attraction be- 
tween virtual transition and the direct, lowest energy 
transitions on different sites. The corresponding coupling 
term - Jy |2i2 3 -)(3jlj| +H.C. - leads to the following con- 
tribution to the free energy in the lowest order: 



5R 



rT 



-"«(! + . 



-0OJ, 



oji tanh(/3wj/2) — ej tanh(/?ej/2) 

At subplateau temperatures, when /3wj S> 1, tanh(/3wi/2) 
can replaced by unity. Furthemore, the summation 
with respect to ej is no longer cut off by the temper- 
ature and the respective integral (weighted by n(e) = 
JL e — l«l/r s ) pjcjjg U p m ost of its value at e ^> T. There- 

fore tanh(/3ej/2) may be replaced by unity as well. (Ac- 
tually, both of those replacements must be made simul- 
taneously lest the sum becomes potentially ill-behaved 
when uji ~ e.) As a result, the expression in Ea. l|77|) 
simplifies: 

SF rT = -J24 ,\ e . 1 . ( 7 «) 

The ej integral is related to an exponential integral E\ 
and yields in the two lowest orders: {hi{T g /uii) — 7^), 
where 7^ = 0.577... is the Euler constant. As in the 
previous calculation, we regard u>i and Jy as uncor- 
related. The summation over uii can be approximately 
represented as a continuous integral between and l ma x 
and leads to a quantity that scales as the area of the do- 
main wall with a logarithmic correction. The final result 



is 5F rT /V = 0.5X(a/£) 4 ln 



2.0 



jd \a J 



1/4' 



(-1 + 2^) 



Up to a logarithmic correction, the expression is inde- 
pendent of the energy parameters in the problem and 
thus must scale linearly in T. Note that we have written 
out the full expression of 8F t t /V that includes the big- 
ger, SV independent term "—1", for the following reason: 
This larger negative term is linear in temperature, which 
apparently would lead to a non-zero (positive) entropy 
at T = 0. This observation signals a breakdown of a per- 
turbative picture of largely non-interacting two-level sys- 
tems. For the sake of argument, let us estimate at what 



temperature this breakdown occurs we compare the mag- 
nitude of the SFj-t/V term, assuming it is correct, to the 
free energy of non-interacting two-level systems per unit 
volume: J 7j^e~ t l T <> [— Tln(l + e _,9£ )], where we have 
appropriately chosen E — as the reference energy. The 
latter expression is equal to (ir 2 /12)T 2 /Tg^ 3 and becomes 
smaller (in absolute value) than the 5F r x /V term at tem- 
peratures below 10~ 3 T g . This temperature is actually 
less, but still within an order of magnitude from the lower 
end of the plateau, which is well within the empirical 
validity of the non-interacting two-level systems regime. 
Let us recall, however, that a perturbative expansion is 
an asymptotic one and therefore always overestimates the 
magnitude of a correction (we suspect that most of the 
error comes from the low e two-level systems). There- 
fore, a more accurate estimate would probably yield a 
break-down temperature lower in value than the estimate 
above. There is a reason to believe the "break-down" 
temperature is just at the edge of the lowest temperatures 
routinely accessed in the experiments. This is suggested 
by several experiments such as on internal friction where 
deviations from the standard non-interacting two-level 
system picture have bee n seen (see, for exa mple, a recent 
review by Pohl et al. llPohl et all 12002^ . In general, 
the effect of interaction between two-level systems could 
exhibit itself under several guises. One of those is an 
apparent gap in the excitation spectrum of the effective 
individ ual TLS. Such effects may have in fact been ob - 
served l|Lasiaunas et all Il978t iThompson et all l2000() . 
The estimates above show thise effects are more likely 
to be observed in substances with a higher glass tran- 
sition temperature, such amorphous silica, or, germania 
(GeC>2). Note, however, that the effects of interaction 
on the apparent TLS spectrum must be separated from 
quantum effects of level repulsion on each sites, that we 
have considered in Section[V] At any rate, the volume ex- 
pansion coefficient, corresponding to the computed value 
of the ripplon-TLS term, is approximately equal to 
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T 
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1/4 



(79) 

Substituting the numerical values for a-SiC>2 in Eqs.JJHJ 
and (|75|l shows that at 1 K, the ratio of the ripplon- 
TLS contribution to the ripplon-ripplon term is about 
1.2 - that is they contribute comparably to the "contrac- 
tion" free energy at this temperature. However, since 
the ripplon-TLS a is temperature independent, it will 
dominate at sw&Kelvin temperatures. The Griineisen pa- 
rameter's value corresponding to Eq.J75J) is 
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1-0^ (I 1 in 
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2.0^ 
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1/4' 



(80) 



The ripplon-TLS term, as estimated here, therefore 
seems somewhat larger relative to the ripplon-ripplon 
term than seen in experiment, consistent with our earlier 
notion that it is somewhat overestimated. Still, qual- 
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itatively our estimates are consistent with the observed 
tendency of 7 to increase in magnitude, when the temper- 
ature is lowered. We point out that the results obtained 
above disregard possible effects of a specific distribution 
of A that will influence the precise value of the coupling 
between phonons and tunneling centers. 

Note that the heat capacity like expression reflecting 
the number of thermally active sites i enters into the 
expressions from Eqs.JTSJ) and (|80() in a linear fashion. 
Therefore, in contrast to Ea. l|7()l) . the temperature de- 
pendence of expression 18Ufl is expected to be largely 
independent of the exact T-scaling of the heat capac- 
ity. Therefore, according to Eq.JSUJ), the Griineisen pa- 
rameter should eventually scale as 1/T at low enough 
temperatures in all substances (however, unrcalistically 
long observation times may be required to verify this pre- 
diction; see the discussion at the end of this section). 
And again, the apparent density of states of the tun- 
neling centers may be modified at those low tempera- 
tures due to interaction e ffects (such as the Burin-Kagan 
(|Burin and Kaganl Il996j) effect). 

According to Eas. (|76|l and (|80|l . the dimensionless con- 
tribution of the attractive forces between the tunneling 
centers can be expressed in a simple manner through the 
Tg/To and T g /T ratios, as well as the relative size of 
the mosaic. Note that the effects of varying the quench- 
ing speed of the liquid on the number in Eas. (|76l) and 
(|80|l add up. For instance, making the quenching faster 
will increase T g and decrease £. The ripplon-TLS term is 
especially convenient with regard to testing our results, 
because it is nearly insensitive to changes in the Debye 
temperature potentially induced by altering the speed of 
glass preparation. 

Finally, we show that the second order coupling be- 
tween direct tunneling transitions is subdominant to the 
already computed quantities. Consider an interaction of 
the form Jij\lj2j) + H.C.. If one repeats simple- 
mindedly the steps leading to Ea. (|73|l . one obtains the 
following simple expression for the free energy correc- 
tion due to interaction between the underlying structural 
transitions: 



TT 



E 2 ej tanh(/3ej/2) 



e i tanh(/3e i /2) 



(81) 



Assuming, again, that Jy and e,-'s are uncorrelated, the 
e summation can be performed via averaging with re- 
spect to the distribution from Eq. l|34l) . One can show 
that the low temperature expansion of the expression 
above yields, within two leading terms, 8Ftt/V = 

-(2T,/3£ 3 ) (f) 6 [1 + (nT/T g )Hn(T g /T)/3}. The T- 

independent term in itself is curious in that it is a con- 
tribution to the "vacuum energy" of the lattice that is of 
purely glassy origin and is entirely due to the locality of 
the free energy landscape of a liquid. Indeed, as attested 
by its scaling with T g /£ 3 , this "vacuum energy" contribu- 
tion would disappear at the ideal glass transition at which 
the whole space is occupied by a non- extensive number 



of distinct aperiodic solutions of the free energy func- 
tional. However, this constant term will have no effect 
on the thermal expansion. The lowest order T-dependent 
term - T 2 In T - actually has a slightly stronger temper- 
ature dependence than the ripplon-ripplon contribution, 
however the latter is larger by at least three orders of 
magnitude, mostly owing to the large number of ripplon 
modes. Apropos, we would like to stress again that the 
presence of vibrational modes of the (extended) mosaic 
walls is essential to the existence of the negative ther- 
mal expansivity effect that we just estimated. There- 
fore, while the present theory predicts that many (and 
most conspicuous) effects that distinguish amorphous lat- 
tices from crystals should be described well by a set of 
non-interacting two-level-like entities at cryogenic tem- 
peratures, the intrinsic multilevel character of the struc- 
tural transitions, that follows from the present theory, in 
glasses exhibits itself even at these low energies in higher 
order perturbation theory. 

To complete the discussion of the second order inter- 
action between tunneling centers we note that the corre- 
sponding contribution to the heat capacity in the lead- 
ing low T term comes from the "ripplon-TLS" term and 
scales as T 1+2q , where a is the anomalous exponent 
of the specific law. Within the approximation adopted 
in this section, a = 0. However it is easily seen that 
the magnitude of the interaction induced specific heat 
is down from the two-level system value by a factor of 
10.(a/£) 5 (<ii/a) 2 ~ 10~ 4 ' 5 and therefore may be safely 
neglected. 

We have so far considered the second order part of the 
induced interactions (square in J^, but forth order in 
g). There could be also, a priori, lower order contribu- 
tions - first order in g, and first order in Jy. First, let 
us consider the term linear in , which also has to do 
with interaction, mediated by the phonons. If non-zero, 
it could be of either sign. In our case, it is identical l y zero 
for the following reas on. It is known ijNeu et all Il997t 
lYu and Leggetttll988j) . that the apparent TLS's are only 
weakly interacting (one could also infer this implicitly 
from the smallness of the second order term that we have 
already estimated. The first order term, if non-zero, is 
comparable to the second order in a mean-field disor- 
dered system. The dipole-dipole 1/r 3 interaction is long 
range and is indeed well described by the mean field). 
But we are dealing here with a non-polarized state, for 
which the first order term, linear in the average on-site 
magnetization, vanishes. In any case, even if the sys- 
tem were in a "ferromagnetic" state, the first order term 
would still be only very weakly temperature dependent 
and thus would not contribute to the thermal contrac- 
tion. Whether to consider such first order term non-zero 
or not is, to some degree, a matter of choice. If non-zero, 
it must be simply thought of as the effective Weiss-like 
field that is part of molecular field at each site. That field 
implies a hard gap of the order T g and indeed is negli- 
gible at low T . Yet, at low enou gh temperatures - mi- 
croKelvins or so l)Neu et a/.lll997j) . the phonon-mediated 
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first order interaction between the tunneling centers may 
become important and one can no longer use the bare 
frozen-in values of the on-site TLS energies, but those 
determined by the interaction. In this regime an inde- 
pendent two-level system picture breaks down and more 
compli cated renormalized excit ations may begin to play 
a role ijBurin and Kagarlll996fl . 

On the other hand, the other possible contribution to 
a, a term linear in g does not have to do with interac- 
tions between the anharmonic amorphous solid excita- 
tions but is due to the direct coupling of the tunneling 
centers with the phonons. This direct T LS-phonon inter- 
action has so far been the main susp ect l|Galperin et all 
11983 iPapoulait Il972t iPhillmsl [l973) behind the anoma- 
lous thermal expansion properties of the glassses. This 
mechanism req uires however the existence of a correlation 
iPhillmsllT97l ( in our notation) between the on-site val- 
ues g and e, or else between A and dA/d(j>a. In other 
words, the value of either classical or quantum splitting 
of a two-level system must be correlated with the way its 
energy changes when elastic stress is applied locally. The 
A with dA/dipa correlation has been argued to make 
a small contribution relative to the g versus e correla- 
tion because of the smallness of the value of the A's 
for th e majority of the thermally active TLS ijPhillipj , 
Il973|) . On the other hand, a correlation between g and 
e could produce, in principle, both a negative or pos- 
itive Griineisen parameter and therefore could explain, 
by itself, the observed variety of expansion anomalies 
in the low T glasses. However, the degree of correla- 
tion between g and e and its temperature dependence 
is not really known and has to be parametrized. The 
soft-potential model offers enough richness in behavior 
to accomodate two possible contributions - one dilating 
and the other contra cting - to the sample 's volume. In 
fact, Galperin et al. l|Galperin et aZJ . fl985^ suggest that 
those two types of the TLS may well be the two types 
of the tunneling cent ers that were postulated ea rly on by 
Black and Halperin ijBlack and HalperinL Il977|) in order 
to resolve the apparent discrepancy between the value of 
the TLS density P as deduced from the phonon scatter- 
ing experiments and the equilibrium and time dependent 
heat capacity measurements. This, of course, could be 
checked experimentally by comparing the degree of the 
discrepancy in P and the sign of the thermal expansion 
coefficient in different substances. (We have shown in 
the previous section how the Black-Halperin paradox is, 
at least partially, explained by quantum corrections to 
the semi-classical landscape picture of structural transi- 
tions in glass.) With regard to the linear in g effect, we 
suggest here a modification to the original argument of 
Phillips from l|Phillipslll973|) . According to Phillips (note 
some notational differences), I7I = 12gaoln2/ir 2 kBT (he 
also assumed a linear heat capacity linear in T). Here, 
|ao| < 1 is an (unknown) coefficient that reflects the de- 
gree of correlation between g and e: {eigi{ei)) = aotig- 
cxq = ±1 means complete correlation and ao — means 
no correlation. Now, due to symmetry, ao must be odd 



power in e, the dominant term being therefore linear (see 
the form of ao{e), somewhat cryptically mentioned as a 
remark of B. Halperin, at the end of Phillips' article). We 
must note, that although we have pretended, within our 
one-polarization phonon theory, that g is a vector quan- 
tity, it is in reality a tensor, if the phonons are treated 
properly. The off-diagonal terms, corresponding to in- 
teraction with shear, will indeed be uncorrelated with 
e due to symmetry. However, the trace of the tensor, 
corresponding to coupling of the TLS to a uniform vol- 
ume change could be, in principle, correlated with the 
energy of the transition. For example, it may happen 
that when the sample is locally dilated, the structural 
transitions in that region will require less energy to oc- 
cur. At present, we do not have an argument in favor of 
or against such a correlation. Note, however that at the 
glass transition temperature, when the current arrange- 
ment of the defects freezes in, most structural transitions 
involve a thermal energy around T g . On the other hand, 
the energy spliting e of the tunneling centers relevant at 
the cryogenic temperatures is significantly smaller. In- 
formally speaking, relative to the thermal energy scale 
at T g all two-level systems with low splitting will feel the 
same to the phonons. Therefore, qualitatively, the corre- 
lation factor ao should be at least a factor of e/T g down 
from the largest value of one. Note, this coincides with 
the form ao(e) oc e, suggested by Halperin. Therefore, 
the contribution of TLS-phonon coupling to the thermal 
expansivity of Phillips (who left the issue of the degree 
of correlation open at the time) should be multiplied by 
a factor of T/T g . This takes into account, in a very naive 
way, both the symmetry and our knowledge of the energy 
scales relevant at the moment of the tunneling centers 
formation. This modifies Phillips' result to yield 
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The temperature independence of this contribution to 
the Griineisen constant is the main difference between 
Ea. (|82|) and the original calculation of Phillips. The nu- 
merical value of the expression should be nearly the same 
for all substances and is about 8. This suggests that the 
direct coupling to phonons is a potential contributor to 
the elastic Casimir effect at temperatures around 1 K. Re- 
member, however, the sign of the expression in Ea. H82|) 
is unknown and its numerical value of 10 1 only provides 
an estimate from the above. 

From the qualitative analysis in this section, we con- 
clude, tentatively, that there are several contributions of 
comparable magnitude to the thermal expansion at low 
temperatures. Higher order effects may also be present. 
In this case, it may be more straightforward to estimate 
the interaction between ripplons as extended membranes 
without using a multipole expansion, as indeed is done 
when computing the regular Casimir force between ex- 
tended plates. 

The qualitative treatment above of the second-order 
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interaction between the ripplons on different sites can be 
extended to higher temperatures as well. It is easily seen 
from Eq.JJsU that an excitation of energy u>i will con- 
tribute only f3Jfj at temperatures comparable to cjj and 
above. Therefore one might expect that at the tempera- 
tures near the end of the plateau the ripplonic transitions 
become thermally saturated and this attractive mecha- 
nism becomes increasingly less important. The expres- 
sion in Ea. (ff7jl . in contrast, is subject to thermal sat- 
uration to a lesser degree. Still, we have seen that its 
scaling with temperature is subdominant to the ripplon- 
ripplon term at temperatures above 1 K. Finally, we re- 
mind the reader about the effect of mosaic stiffening ex- 
plained in the previous sections. This should also dimin- 
ish the attraction between the tunneling centers, owing 
to a smaller number of resonant modes at the sites of cen- 
ters thermally active at these higher temperatures. On 
the other hand, the usual anharmonic effects also become 
more significant at a higher T leading to a turn-over in 
the temperatrure dependence of a, as circumstantially 
sup ported by the old data on se veral materials cited in 
Ref. ljKrause and KurkiianL ri968) . However, in order to 
assess this "cross-over" temperature, one needs to know 
the magnitude of the regular thermal expansion due to 
the non-linearities of the lattice. This is something that 
would be extremely difficult to measure independently, 
because even a crystal with the same stoichiometry as 
the respective glass, is not guaranteed to have the same 
non-linearity. Direct computer simulation estimates of 
the Griineisen parameter, on the other hand, may be 
problematic due to the current difficulty of generating 
amorphous structures corresponding to realistic quench- 
ing rates. This is the main reason we have confined our- 
selves here to sub-plateau temperatures. 

Finally, we note again that even at the low tempera- 
tures we have been discussing, not all glasses have been 
shown to exhibit a negative a. According to our the- 
ory, however, the "Casimir" contribution to a is nega- 
tive and sub-linear in T, whereas the regular non-linear 
expansion coefficient is positive but only cubic in tem- 
perature. Therefore, there should be a (perhaps very 
low) temperature at which the Casimir force should dom- 
inate. Data for many substances, although still positive 
at the achieved degree of cooling, do extrapolate to neg- 
ative values of a at finite tempe ratures. This is not th e 
case, however, for all substances l)Ackerman et aZl ll9841. 
Even excluding the possibility of error in these difficult 
experiments, this is not necessarily inconsistent with our 
theory for the following reasons. As the temperature is 
lowered, it takes a long time (proportional to T~ 3 ) for the 
tunneling transitions to occur and appear thermally ac- 
tivated. For these same reasons, like the amorphous heat 
capacity, the direct interaction effect is time dependent 
at low temperatures. It may therefore take an exces- 
sively long time to actually observe the effects, discussed 
in this section, at very low temperatures, thus making it 
difficult to see a sign change in a for lattices with rel- 
atively large anharmonicity. Incidentally, this analysis 



predicts that the response of the length of an amorphous 
sample to a temperature change at sub-plateau temper- 
atures must be time dependent (such time-dependence, 
acompanied by h eat release, has been observed in poly- 
crystalline NbTi ijEscher et a?! I200C|) 1. Since the inter- 
action effect is quadratic in concentration, one expects 
qualitatively that the relative rate of the expansion's time 
dependence should be twice that of the specific heat. 



VII. CONCLUSIONS 

In summary, this work elucidates the origin of the ther- 
mal phenomena observed in the amorphous materials at 
temperatures ~ Td/3 and below, down to the so far 
reached milliKelvins. The nature of these phenomena 
can be boiled down to the existence of excitations other 
than elastic strains of a stable lattice. The peculiarity of 
these excitations is exhibited most conspicuously in the 
following phenomena: The specific heat obeys a nearly 
linear dependence on the temperature at the lowest T, 
greatly exceeding the Debye contribution. At the same 
temperatures, the heat conductivity is nearly quadratic 
in T and is universal if scaled in terms of the elastic con- 
stants. At higher temperatures (~ Td/30), the density of 
these mysterious excitations grows considerably leading 
to enhanced phonon scattering and thus a plateau in the 
temperature dependence of the heat conductance. This 
increase in the density of states is also directly observed 
as the so called Boson Peak in the heat capacity data, as 
well as inelastic scattering experiments. 

We have argued that the origin of these excitations is 
a necessary consequence of the non-equilibrium nature 
of the structural glass transition. This transition, not 
strictly being a phase transition at all in a regular equi- 
librium sense, occurs if the barriers for molecular mo- 
tions in a supercooled liquid become so high as to pro- 
hibit any macroscopic shape c hanges in the material o n 
the scales of hours and longer IjXia and Wolvnest (2000) . 
The origin of these high barriers lies in a cooperative 
character of the molecular motions, which involve around 
200 molecules at the glass transition temperature. Un- 
like regular crystals, where the correlation between the 
molecular motions is rather long range, thus leading to 
the emergence of translational symmetry below solidifi- 
cation, the motions within the cooperative regions in a 
supercooled liquid, or entropic droplets, are only weakly 
correlated with their surrounding. In the language of the 
energy landscape paradigm, a crystal is a (possibly non- 
unique) ground state of the sample (thus the long-range 
correlation!), whereas a glass is caught in a high energy 
state, not being able to reach the true ground state for 
kinetic reasons. The respective dense energy spectrum at 
these energies exhibits itself in the existence of alterna- 
tive mutually accessible conformational states of regions, 
or domains, of about 200 molecules in size. It was ar- 
gued that quantum transitions between these alternative 
states are the additional excitations observed in glasses 
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at low temperatures. The knowledge of the spectral and 
spatial density of these excitations allowed us to estimate 
from first principles the magnitude of the observed linear 
specific heat. The relevant energy scale here is the glass 
transition temperature T g itself. 

Stability requirements for the existence of these alter- 
native conformational states at T g allowed us also to esti- 
mate the strength of their coupling to the regular lattice 
vibrations, which is determined by T g , the material mass 
density and the speed of sound. This enabled us to un- 
derstand the universality of the phonon scattering at the 
low temperatures. 

The novelty of this picture is that we have established 
rather generally a multiparticle character of the tunnel- 
ing events. This is counter-intuitive because, naively, the 
larger the number of particles involved in a tunneling 
event, the larger the tunneling mass is, and the harder 
the tunneling becomes. This is indeed the case for sys- 
tems like disordered crystals or crystals with substitu- 
tional impurities, where the tunneling mass is that of 
an atom, and the barrier heights are determined by the 
energy of stretching a chemical bond by a molecular dis- 
tance; this virtually excludes the possibility of tunneling. 
The existence of structural rearrangements in a macro- 
scopically rigid system is a sign of the system being in 
a high energy state in which the available phase space 
is potentially macroscopically large. However, a decrease 
in this density of states for glass transitions occuring at a 
slower pace of quenching would result in the necessity to 
engage a larger number of atoms in these structural rear- 
rangements. Transitions between the internal states of a 
domain involve only a very minor length change of each 
individual bond and atomic displacements not exceeding 
the Lindemann length, which is of the order one-tenth of 
the atomic length scale. It is not particularly beneficial 
to picture the tunneling events as individual atomic mo- 
tions but rather as the motion of an interface between the 
alternative states of the domain. This domain wall is a 
quasi-particle of a sort, which has a low mass indeed: per 
molecule in the domain, it is only about one-hundredth 
of the atomic mass. The contributes to the ease of the 
tunneling events that are thermally relevant at cryogenic 
temperatures: These events are subject to only very mild 
potential variations and are possible, again, because the 
lattice is frozen-in in a high energy state. 

The spatially extended character of the domain wall 
excitations along with their strongly anharmonic nature 
explains also higher temperature phenomena, such as the 
Boson peak and the plateau in the heat conductivity. 
By using our knowledge of the surface tension and the 
mass density of the domain wall we were able to calcu- 
late the energy spectrum of vibrational excitations of the 
active domain walls, or ripplons. This spectrum is in 
good agreement with the observed frequency of the Bo- 
son peak. The ripplonic excitations accompany the tran- 
sitions between the domain's internal states and thus are 
strongly coupled to the phonons. This has enabled us 
to understand the experimentally observed rapid drop in 



the phonon mean free path at the plateau temperatures. 
In addition, we have investigated the effects of phonon 
coupling on the spectrum of the ripplons. These spec- 
tral shifts scale with T g and seem to be the cause of the 
non-universal position of the plateau. 

We have carrried out an analysis of the multi-level 
structure of the tunneling centers that goes beyond a 
semi-classical picture of the formation of those centers at 
the glass transition, that was primarily employed in this 
work. These effects exhibit themselves in a deviation of 
the heat capacity and conductivity from the nearly lin- 
ear and quadratic laws respectively, that are predicted 
by the semi-classical theory. 

A Van der Waals attraction between the domain walls 
undergoing tunneling motions was argued to contribute 
to the puzzling negative expansivity, observed in a num- 
ber of low T glasses. 

Finally, we note that the conclusions of this work 
strictly apply only to glasses made by quenching a su- 
percooled liquid. One may ask, nevertherless, to what 
extent the present results are pertinent to other types 
of disordered solids, such as "amorphous" films made 
chemically or by vapor depositions, or, say, disordered 
crystals. Indeed, phenomena, reminiscent of real glasses, 
such as an excess density of states, are observed in many 
types of disordered materials, although they do not ap- 
pear t o be as universal as in true glasses (see, for ex- 
ample (|Pohl et all 12002^ . In this regard, we note that 
most of the phenomena discussed in the present work 
should indeed take place in other types of aperiodic struc- 
tures. What makes quenched glasses special is the intrin- 
sic character of their additional degrees of freedom that 
stems from the non-equilibrium nature of the glass tran- 
sition. Since the characteristics of this transition (while 
not being a transition in a strict thermodynamic sense!) 
are nearly universal from substance to substance, many 
low (and not so low) temperature properties of all those 
substances can be understood within a unified approach. 
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APPENDIX A: Rayleigh Scattering of the Phonons due to 
the Elastic Component of Ripplon-Phonon Interaction 

In this Appendix, we present an argument on the 
strength of the phonon scattering due to the direct cou- 
pling with the ripplons via lattice distortions, but not due 
to the inelastic momentum absorbing transition in which 
the internal state of the domain changes. We thus con- 
sider phonon scattering processes which do obey selection 
rules and couple to the lattice strain only in the second 
and higher order. This scattering is of the Rayleigh type 



50 



(and higher order) and occurs off the domain walls as 
localized modes. Importantly, we will use only derived 
quantities and no adjustable parameters in this estimate. 
We show here that, indeed, this absorption mechanism is 
not significant compared to the resonant scattering by 
the inelastic transitions between the internal states of a 
thermally active domain. 

First, it proves handy to rederive the ripplon spectrum 
from Eq. jS^I in the less general case p g = (but non-zero 
pressure!). As argued in Section HV1 the droplet wall is 

O o 1/2 

at equilibrium pressure p = ^ — | ^ 3/2 . If the surface 
is distorted locally by fl, this results in an extra force 
on this portion of the wall due to a changed curvature 
ijMorse and Feshbachl Il953). The second Newton's law 
(as applied per unit area) yields then: 
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where 9 and (f> are the usual polar and asimuth angular 
coordinates on the surface and we took into account the 
r dependence of pressure. The equation above can be 
solved by a linear combination of the eigen-functions of 
angular momentum in 3D: 



(A2) 



Yi m (9, </>) are the spherical Laplace functions (m = — 
Substituting a harmonic of Z-th order in Ea. l|Al(l yields 
the equation for u>i derived in text as Ea. (|32|) . We will 
absorb the 9/8 factor into the definition of a in the rest 
of the Appendix. 

A (fake) potential energy, yielding the equation of mo- 
tion (|A1|) . is (c.f. the dis cussion of surface waves on a 
spherical liquid droplet in ijLandau and Lifshitd . Il98^) ): 
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Although varying Eci . (|A3|I w.r.t. fl does produce the 
Ea. (|Al|l . note that it differs (by a factor of 9/8!) from 
the original surface energy a4irr 2 . The resulting error is 
sufficiently small for our purposes, however this subtlety 
may be worth thinking about as this could reveal an extra 
friction mechanism due to the wetting phenomenon and 
surface tension renormalization mentioned in our discus- 
sion of the random first order transition in Section ILTI 
While the domain wall positions are not strictly tied to 



the atomic locations, they are tied to the lattice as a con- 
tinuum and follow the lattice distortions. Let us employ 
our usual "scalar" phonons descibed by Hamiltonian 



pc 2 s {V^f 
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where [ip(ri), 7r(r2)] = iH6(ri — r2). The surface energy 
due to the presence of both and ip is: 



J 



H, 



surf 



d{cos9)\(R+[^~*jj(n)} + n) 2 



d(4> + n) 

89 



i fd{ij + n) 



sin 2 9 V 



(A5) 



where ip is taken on the sphere of radius R with the cen- 
ter located at r^. The potential energy in Ea. (jA5|) thus 
provides an explicit form of phonon-ripplon interaction 
due to the liquid free energy functional solutions being 
imbedded in the real space. 

If we expand the value of the displacement field <p 
in terms of spherical harmonics according to ipi m = 
J d(j)d{cos 9)4>(r = R)Y l * m ((j),9), it is then possible 
to write down equations of motion for the (l,m)- 



components of both ripplon and phonon displacements: 



d 2 n 



I in 



dt 2 



+ Ujf (Cll m + i>l m ) = 0. 
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The equation of motion for the phonon field can be 
obtained e.g. from ip — i[Hph + H sur f, Tt/p] to yield: 
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The terms with -0 on the r.h.s. serve only to modify 
the local elastic constants, and therefore give rise to the 
regular Rayleigh scattering, so we will ignore them from 
now on. 

Equations l|A6IA7fl can be used to write down equa- 
tions of motion for the retarded Green's functions, which 
are prefera ble due to their convenient analytical prop- 
erties (see ijZubarevl ^960) for our conventions). We 
are interested in the system's response to "plucking" 
the latice at site r = at time zero, hence the choice 
of the Green's function corresponding to an operator 
X: -i0(t - t') {[X{t), i/>(r = 0, t' = 0)]). Eas. (|A6IA7jl . 



J 



if rewritten for the corresponding Green's functions, 
will preserve except there will be an additional term 
— - 5(t)S 3 (r), corresponding to the "plucking" event, in 
the r.h.s. of Eci. (|A7() (note also a change in units). Thus 
obtained equations are possible to rewrite in the Fourier 
space: 
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where ipl 
and 



nm = Jd s k^e^(4ny 3l (kR)Y* n (k/k) 
we used the expansion of a plane wave 
in terms of the spherical harmonics: e lkr = 

4^ EZo TL=-l il h ( kr ) Y i* m WWhn (r/r) . Here, 
jl(x) = y // n/2x Ji+i/2{ x ) is the spherical Bessel function, 
which scales as x l for small x, hence we see that the 
ripplons' coupling with the phonons is quadratic or 
higher order in k as the second harmonic is the lowest 
order term allowed. Modes I = and I — 1 have 
the meaning of the droplet 's growth and translation 
respectively, as was discussed in Section IIV.CI These 
modes are not covered by this Section's formalism. Even 
though the theory as a whole could be thought of as a 
multipole expansion of a molecular cluster interacting 
with the rest of the lattice, the modes of different orders 



J 



end up being described by different theories. 



The system of Eqs. (|A8ll and (|A9|) can now be used 
to determine the sound dissipation due to the interac- 
tion with the ripplons. Since the system is infinite and 
has a continuous spectrum, all excitations will have fi- 
nite life-times, which can be, in principle, obtained self- 
consiste ntly by using e.g. the Feenberg's perturbativ e ex- 
pansion jAb ou-Chacra et OR 119731 iFeenberd. ll94Sft (one 
in the end arrives at Green's functions that are well be- 
haved at infinity, as implied in the thus greatly simplified 
derivation). We do not have to do this self-consistent 
self-energy determination as long as we are interested in 
the lowest order estimate, as justified in the end by the 
smallness of the obtained value of the perturbation. Sub- 
stituting Eq.JXSJl into Eq.l|A9|l yields 
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Since the spatial locations of active droplets are not 
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correlated 8 , we can replace the summation over the 
droplets by a continuous integral, assuming at the same 
time that the ripplon frequency corresponding to ioi 
varies from droplet to droplet within a (normalized) dis- 
tribution Vi(lu) centered around tui and having a charac- 

8 This is not strictly true - they, of course, can not be on top of 
each other. 

I 



where n is the concentration of the active domain 
walls to be estimated shortly and we have displaced 
u by e into the upper half-plane because we are look- 
ing for the retarded Green's function. Also, in or- 
der to derive Ea. ljAll|) . we have used the summa- 
tion theorem for the spherical functions P;(nn') = 

2TfTEL=-i^( n ')^™(n), as well as fl(-l) = (-1)', 
where Pi is the Legendre polynomial. If we ignore the 
real part of the r.h.s. of Ea. (|All|l . responsible only for 
the dispersion, the poles of the resultant phonon Green's 
function are found by solving ui 2 — c 2 k 2 + i2u)T~ 1 = 0, 
where r^ 1 clearly has the meaning of the inverse life-time 
of a phonon of frequency lo and is given by 

9 

r- 1 =n-^ 7 r 2 [2 + ?(/ + l)](2/ + l) J f(fc J R)^H, (A12) 

where we have ignored the contribution of the peaks cen- 
tered around (— u>i). We remind the reader that l m ax — 9 
is dictated by the finite size of a droplet. 

One can find the valu e of 5loi from an argumen t iden- 
tical to the one used in (|Xia and WolvnesL l2001a|) to ob- 
tain the width of the distribution of the barriers for the 
droplet growth free energy profile. At the glass tran- 
sition, a liquid breaks up into dynamically cooperative 
regions, so that a translation of one atom involves mov- 
ing about 200 atoms around it, which involves overcom- 
ing a large (on average) barrier. This barrier's height 
is determined, together with the domain surface ten- 
sion coefficient, by the configurational entropy density, 
which in its turn reflects the number of metastable states 
available to a particular volume of liquid at this tem- 
perature. Even though a good description of freezing 
is achieved by assuming that this number of available 
states does not strongly depend on where exactly on 
the free energy surface a particular molecular cluster is 
ijXia and Wolvnesl l2000|) , it should vary from domain to 
domain. The size of the variation can be estimated from 
the known magnitude of the entropy fluctuations at con- 
stant energy, so that the ratio of the variance to the mean 
is related to the jump in the heat capacity at T g and 
subsequently turns out to be 1/2^/T) ijXia and Wolvnesl 



teristic width 5loi, whose value will be discussed shortly. 
There is no reason to believe that the frequency and lo- 
cation of the tunneling centers are correlated, therefore 
one obtains 



(All) 

I 

l2001a|) . where D is the liquid's fragility, entering the 
Vogel-Fulcher law for relaxation times in a supercooled 

liquid T re i axat ion oc e T ~ T K . We conclude then that the 
lower bound on the fluctuations of the ripplon frequency 
uji is given by Sloi ~ loi/2\/~D. 

Lastly, in order to use Eq. (|A12(1 to compute the phonon 
absorption due to this particular mechanism, we need 
to estimate the density of the active domain walls. It 
will suffice for our purposes here to consider as active 
the defects that contribute to the specific heat, that is, 
roughly, n ~ -^T ' jT g . A more accurate estimate would 
be similar to the one we made when calculating the bump 
in the heat capacity in Section II V. 01 

We are now ready to give a numerical estimate of the 
expression in Eq. (|A12|I . We will compute here the contri- 
bution of the I = 2 term in the plateau region. It is conve- 
nient to represent kR from Eq. (|A12(1 as kR ~ Q 4 ( a /^ D ■ 
For the reference, (a/£)Trj ~ 0.2Tq is at the high tem- 
perature end of the plateau, whereas its middle is about 
an order of magnitude lower depending on the substance 
(see k vs. T/Trj plot in FigQ. We can now use our usual 
expressions connecting <j,T g ,LOD,c s , p,a etc to obtain a 
numerical estimate of Ea. (|A12|l at the plateau frequen- 
cies UJplateau 

~ 10 ujd- Even if one favorably assumes 
that Cl>2 ~ ^plateau (it is somewhat larger according to 
Section llV.Dfl . one still gets Z m f p /A >^ 10 4 at the plateau 
frequency, whereas the resonant absorption by the TLS 
would give Z m f p /A ~ 10 2 . The amplitude of this type of 
absorption is small due to the weakness of direct coupling 
to the ripplons for the processes not accompanied by a 
change in the domain's internal state. 



APPENDIX B: Frequency Cutoff in the Interaction 
Between the Tunneling Centers and the Linear Strain 

As argued in Section IIII.BI the coupling of the tun- 
neling transition to a phonon can be found from an ad- 
ditional energy cost of moving the molecules within the 
domain in the presence of a strain and is given by an 
integral over the droplet 's volume (we consider only lon- 
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gitudinal strain for simplicity): 

g = pc 2 s f d 3 r(V<?)(Vd), (Bl) 
Jv 

where pc 2 s is basically the elastic modulus, (j) arl d d are 
elastic and inelastic components of the atomic displace- 
ments respectively. If the phonon's wave-length is much 
larger than £, the elastic component is constant through- 
out the integration region and the integral reduces to 
one over the droplet 's surface and thus the g estimate 
obtained in text. Otherwise, one obtains: 

g = pc 2 s j dS d - d 3 r (dV)(V<?) j . (B2) 

The volume integral will give a higher order term 
in k, so for now, we focus on the surface integral. 
The displacement due to the phonon is conveniently 
expanded it terms of the spherical waves: e lkr = 
4»rESoli=-i*'ji(fcr)^(V*)^m(r/r). Since it is 
the first derivative with respect to r that we are inter- 
ested in, we only need the / = 1 term from this expansion. 
The angular part contributes only to the overall constant, 
but it is the spherical function ji(fcr) that sets the cut-off 
value of the wave- vector, above which the phonons do not 
produce significant linear uniform stress on the domain. 
In Fig[^ we plot the derivative dji(x)/dx (or, rather, 
we plot the square of it, which enters into all the final 
expressions) . 
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FIG. 24 Shown is the derivative of the 1st order spheri- 
cal Bessel function determining the effective decrease in the 
elastic field gradient produced by a phonon of wave-length k 
(x = kR). 

We see that it is not unreasonable to assume that only 
the phonons with kR <~ 6 will exert an appreciable 
linear strain on the domain. kR = 6 translates into oj c ~ 
2.5(a/pw D . 

While we are at it, we estimate the interaction of the 
domain with the higher order strain, at least due to the 
term I|B1|) . in the frequency region of interest. The next 
order term in the k expansion in the surface integral from 
Ea. (|B2l) has the same structure, but is scaled down from 
the linear term by a factor of kR. At the plateau frequen- 
cies ~ w.d/30, kR < 0.5 as immediately follows from the 
previous paragraph. While this is not a large number, it 



is not very small either. Therefore, this interaction term 
is of potential importance. 

The volume integral in Eq. ljB2fl produces a quadratic 
term, which is roughly equal to 

(V(£) J v d 3 r (d k). We then proceed in a completely iden- 
tical fashion to our earlier estimate of g. Assuming the 
diplacements within the droplet are random, one gets 
for the integral \ \J N* a 3 d^k, where factor of 1/4 comes 
about because the displacement is assumed do decrease 
from dr, in the center of the dro plet to zero at the edge 
l|Lubchenko and Wolvnesl l200l|) . This yields that this 
term becomes compa rable to the linear one at frequen- 
cies lo ~ uu D ^(a/£) 4/ (67T 2 ) 1 / 2 ~ 0.4 uo D - well beyond 
the high T end of the plateu. 

We must note, there are other sources of non-linearity 
in the system, such as the intrinsic anharmonicity of the 
molecular interactions present also in the corresponding 
crystals. While these issues are of potential importance 
to other problems, such as the Griineisen parameter, ex- 
pression (|B1J| only considers the lowest order harmonic 
interactions and thus does not account for this non-linear 
effect. We must note that if this non-linearity is signif- 
icant, it could contribute to the non-univrsality of the 
plateau, in addition to the variation in T g /ujD ratio. It 
would be thus helpful to conduct an experiment com- 
paring the thermal expansion of different glasses and see 
whether there is any correlation with the plateau's loca- 
tion. 
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